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Resume 


Ce sujet de these s’inscrit dans le cadre des recherches visant a obtenir des ve- 
hicules aeriens miniatures ayant a la fois les performances des avions pour le vol 
horizontal et la manoeuvrabilite des helicopteres pour le vol stationnaire. L’objectif 
est de developper, modeliser et commander un drone combinant la manoeuvrabilite 
des vehicules a voilure tournante (helicopteres) telle que l’avance lente, le decollage 
et l’atterrissage vertical, et les performances d’un vehicule a voilure hxe (avions) 
telle que l’avance rapide, la longue portee et une endurance superieure. Ce type 
de recherche a suscite beaucoup d’interet dans le passe car ce type d’apparcil ne 
necessite pas de piste de decollage et sa capacite de vol stationnaire le rend tres 
utile pour des missions de surveillance aerienne. L’objectif du sujet est en effet de 
concevoir et de realiser des prototypes d’un systeme de drone pouvant effectuer 
des decollages/atterrissages verticaux de maniere autonome, puis de realiser une 
transition autonome vers un vol d’avancement rapide. 

Deux prototypes experimental^ ont ete developpes au laboratoire : l’avion conver¬ 
tible et 1’avion basculant tri-rotors. Les modeles dynamiques de chaque vehicule 
ont ete obtenus en utilisant la methodologie de Newton-Euler prenant en compte 
les forces et les couples aerodynamiques. Une loi de commande non lineaire pour 
la stabilisation en attitude et en position d’un corps rigide a ete proposee. Elle 
est basee sur l’approche de commande bornee qui garantit la convergence du corps 
rigide vers une attitude et position desirees. Sa performance est aussi renforcee par 
1’algorithme de generation de trajectoire et les reseaux de neurones. La loi de com¬ 
mande est ensuite appliquee pour commander un quadrirotor, l’avion convertible et 
Lavion basculant tri-rotor. Pour l’avion convertible, la commande d’attitude a ete 
adaptee en utilisant le formalisme des quaternions. Par allelement, un simulateur 
a ete developpe, il permet de developper, d’executer et de tester le programme du 
prototype sur un ordinateur grace a un rnodele dynamique. Du cote experimental, 
l’electronique embarquee a ete congu pour avoir un systeme embarque de com¬ 
mande plus puissant et adaptable. La validation des lois de commande proposees 
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a ete realisee sur la plate-forme experimental de l’avion convertible qui execute 
les algorithmes en temps reel avec une bonne performance. 
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Chapitre 1 


Introduction 


Ces dernieres annees, l’interet pour les Vehicules Aeriens Autonomes ( Unman¬ 
ned Aerial Vehicle, UAV) s’est constamment accru. Aussi appeles drones, ils ont 
d’abord ete employes, principalement dans le domaine militaire, pour des missions 
de reconnaissance, de surveillance et aussi dans un role offensif. De nos jours, nous 
pouvons trouver des applications civiles telles que par exemple la surveillance du 
trafic routier, la surveillance et la protection de l’environnement, la recherche et 
le secours de blesses, ou la gestion de grandes infrastructures telles que les lignes 
haute tension, les barrages et les ponts. Les progres considerables de la technologie 
ont permis d’augmenter la variete et le nombre d’applications de ces vehicules au- 
tonomes. Les chercheurs et les ingenieurs ont travaille de faqon intense au cours des 
deux dernieres decennies pour developper des machines volantes efhcaces qui effec- 
tuent des nombreuses fonctionnalites. Le developpement de ces machines volantes 
implique les domaines de la mecanique, de l’electronique, de l’aerodynamique, de 
Lautomatique et des systemes informatiques embarques. 


1.1 Vehicules aeriens autonomes 

Un vehicule aerien autonome (UAV) est un terme auto-descriptif utilise par l’ar- 
mee americaine pour decrire les dernieres generations d’avions sans pilote. En fait, 
ce terme decrit un objet volant avec un certain degre d’autonomie. Les conhgura- 
tions des vehicules peuvent etre tres variees : helicopteres ou avions partiellement 
radio-commandes, cerfs-volants motorises, missiles de croisiere. Les drones sont de 
plus en plus importants dans le domaine de hespace, au niveau militaire, pour 
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des operations commerciales et civiles. Les progres technologique ont rendu le vol 
sans pilote non seulement possible, mais aussi pratique et economique. Puisque 
les drones ne portent pas les humains, ces vehicules peuvent etre congus pour etre 
avec une petite taille, faciles a stocker et transportes a moindre cout. Les drones 
peuvent egalement accomplir des missions qui pourraient placer les humains en 
danger : par exemple lors de guerres, les Etats-Unis ont developpe des avions sans 
pilote capables de faire de la reconnaissance sans engager de vies humaines sur 
le terrain ennemi. Tous ces avantages font des drones une plate-forme plus eco¬ 
nomique et plus appropriee que les avions pilotes pour de nombreuses missions. 
Finalement, les drones avec leur capacite de voler de fagon autonome, seul ou en 
formation et la capacite de vol stationnaire et de la manoeuvre, represented une 
nouvelle classe de vehicules ayant le potenticl de renforcer la securite nationale et 
d’augmenter le soutien des fonctions publiques. 

De nos jours, les drones sont devenus beaucoup plus accessibles. En effet, la tech- 
nologie evoluee permet la miniaturisation du systeme de drone et la baisse du cout. 
Des nombreuses equipes de recherche travaillent maintenant sur des drones minia¬ 
tures et la production scientihque sur le sujet est abondante. Les applications ne 
sont plus reservees aux domaines militaires mais aussi aux domaines civils. Pour- 
tant, elles restent dans le domaine de l’observation, le transport de marchandises 
ou de personnes par exemple n’etant pas d’actualite. 


Autonomie 

Au debut, certains UAVs ont ete appeles drones, parce qu’ils n’etaient plus simple- 
ment controle par un pilote humain (parfois appele l’operateur) via une radiocom- 
mande. En effet, ils ont integre un certain degre d’autonomie grace a un systeme 
de controle integre et/ou un systeme de guidage pour effectuer des fonctions de 
bas niveau tels que la stabilisation et des fonctions simples de navigation comme la 
poursuite de trajectoires predehnies. L’autonomie du vehicule aerien est un champ 
recemment emergent, dont le support est largement finance par les militaires pour 
developper une technologie utilisee en cas de guerre. En comparant revolution de 
la fabrication du materiel des drones et le marche de la technologie pour l’autono¬ 
mie, ce dernier est assez immature et peu developpe. Dans une certaine mesure, le 
but ultime dans le developpement technologique de l’autonomie est de remplacer 
le pilote humain. Ainsi, l’autonomie est couramment dehnie comme la capacite a 
prendre des decisions sans l’intervention humaine. Le developpement de l’autono¬ 
mie du drone est associe au developpement du domaine de 1’intelligence artihcielle 
tres populaire dans les annees 1980 et 1990, comme les systemes experts, les re- 
seaux de neurones, l’apprentissage automatique, le traitement de la vision. La 
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technologie de l’autonomie importante pour le developpement des futurs drones 
concerne les categories suivantes : 

• La fusion des capteurs : La combinaison de l’information provenant des dif¬ 
ferent s capteurs utilises a bord du vehicule. 

• Communication : Tenir compte de la communication et de la coordination 
entre plusieurs agents en presence d’informations incompletes et imparfaites. 

• Planification de mouvement : Determination d’un chemin optimal pour le 
vehicule aerien en repondant a certains objectifs et contraintes, tels que les 
obstacles, l’economie de l’energie... 

• Generation de la trajectoire : Determination d’une manoeuvre de controle 
optimale a prendre pour suivre un chemin determine ou pour aller d’un 
endroit a un autre. 

• Repartition des taches et de Vordonnancement : Determination de la repar¬ 
tition optimale des taches avec des contraintes de temps reel et de l’equipe- 
ment. 

• Cooperation tactique : Realisation d’une sequence optimale et la distribution 
des activites entre les drones afin de maximiser les chances de succes dans 
tout scenario de mission. 


1.2 Applications 

II existe differents types de drones qui effectuent de nombreuses taches pour des 
applications civiles et militaires. 


1.2.1 Applications civiles 

En ce qui concerne les applications civiles, les UAVs ont la capacite de realiser une 
variete de mission. Les utilisations potentielles des drones dans l’industrie civile 
sont : 

• Dans le secteur agricole : les drones equipes d’engrais et de materiel de disper¬ 
sion des pesticides peuvent etre utilises pour vaporiser ces substances dans 
les champs. 

• Applications industrielles : Ces applications concernent, par exemple, la sur¬ 
veillance d’usines nucleaires, la surveillance des pipelines, etc.. 
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• Controle de l’environnement et la recherche meteorologique : La recherche 
scientifique sur la nature (environnement, atmosphere, pollution, etc..), sur 
les sites archeologiques peut etre effectuee par des drones equipes d’outils 
appropries. 

• Detection des incendies : Les vehicules aeriens equipes des capteurs infra¬ 
rouges peuvent detecter les incendies dans le foret et ainsi les pompiers 
peuvent etre prevenus plus rapidement. 

• Exploration miniere : Les drones sont utilises dans des etudes aeriennes et 
terrestres pour trouver des mineraux sur les regions isolees et difhciles a 
atteindre. 

• Relais de communication : Les drones a haute altitude et longue endurance 
peuvent etre utilises comme des satellites pour la transmission des signaux. 

• Surveillance des cotes : Les drones sont utilises comme garde cote pour les 
missions de defense. 

• Gestion des catastrophes et des accidents : Les plates-formes aeriennes avec 
des cameras de surveillance peuvent fournir des informations en temps reel 
dans des situations dangereuses telles que des seismes, des tsunamis... 

• Recherche et sauvetage : La recherche des survivants dans le naufrage de 
navire, d’accident d’avion, etc.. 

• Defense frontaliere : La patrouille des frontieres grace a des plates-formes 
aeriennes 

• Soutien de la police : Les UAVs peuvent remplacer les hclicopteres de la 
police d’une maniere plus rentable. 


1.2.2 Applications militaires 

Les drones ont la capacite d’executer des missions a des fins militaires. La liste 
ci-dessous presente quelques applications militaires : 

• Cible et leurre : Fournir aux artilleries au sol et aeriennes une cible simulant 
un avion ennemi ou un missile. 

• Reconnaissance : Offrir des renseignements sur le champ de bataille. 

• Combat : Fournir une capacite d’attaque pour les missions a haut risque. 
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1.3 Classifications des drones 

1.3.1 Les drones a voilure fixe 

Les drones de cette categoric sont principalement des avions. Les avions neces- 
sitent en principe une piste pour decoder ou atterrir. Cependant les petits modeles 
peuvent etre lances a la main, en mettant prealablement un niveau de gaz suffi- 
sant. La propulsion de l’avion est assuree par un ou plusieurs moteurs et l’energie 
depensee ne sert qu’a le faire avancer, et non pas a le porter. Par consequent, la 
vitesse d’avance est tres rapide. L’avion est done une configuration a privilegier 
pour des missions de longue duree ou de grande distance. Pourtant, l’utilisation 
de ce type de drone necessite presque toujours un pilote experiments a cause de la 
grande vitesse. Le deplacement des differentes gouvernes permet de controler les 
couples de tangage, roulis et lacet (figure 1.1). 



Figure 1.1 — Le drone Solar-Storm de l’ENAC [72] avec alimentation hybride (solaire et 
batteries) 


1.3.2 Les drones a voilure tournante 

L’exemple general des drones a voilure tournante est l’helicoptere classique : une 
machine volante ayant deux rotors tournants, avec un rotor principal donnant la 
poussee et un rotor anticouple a la queue (voir figure 1.2). Le rotor arriere peut 
etre remplace par un rotor tournant dans le sens inverse au rotor principal (voir 
figure 1.3). Cette configuration est plus compacte et la poussee utile est plus grande 
que celle de l’helicoptere classique. Pourtant, la stabilisation de Tangle de lacet est 
plus compliquee car elle est connnandee par la difference de la vitesse de rotation 
des deux rotors contrarotatifs. Le CoaX (figure 1.3) est un lielicoptere miniature 



1.3 Classifications des drones 


6 


robotique congu pour la recherche et l’education. Le systeme pese environ 300g, 
a une envergure de 340mm et une hauteur de 274mm. Le pilotage est assure par 
deux servomoteurs, qui controlent le plateau cyclique du rotor inferieur. Le rotor 
superieur est attache a la barre de stabilisation agissant cornme un regulateur 
passive de roulis et de tangage de Lhelicoptere. 



Figure 1.2 - Helicoptere classique 



Figure 1.3 — Helicoptere CoaX a deux rotors contrarotatifs [73] 


Les configurations utilisant de multiple rotors permet de simplifier la structure 
mecanique de la creation des forces et des couples. Dans cette categorie de drone, 
un nombre pair de rotor est souvent utilise, par exemple quatre, six ou huit ro¬ 
tors autour de l’axe principal vertical du drone. La raison principale d’utiliser des 
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paires de rotor est que le couple resistant de chaque rotor s’annule en forgant les 
deux rotors a tourner dans des sens differents (contra-rotatifs). En consequence, 
ce type de l’helicoptere peut voler sans tourner antonr de son axe. Cette configu¬ 
ration permet d’eliminer le plateau cyclique ( swash-plate ), qui est une structure 
mecanique tres complexe et couteuse. En effet, la complexity mecanique dn drone 
est diminuee et elle entraine non senlement une baisse dn cout pour fabriquer le 
drone mais aussi une augmentation de la fiabilite dn systeme. Les helicopteres 
a quatre rotors (figure 1.4), a six rotors (figure 1.5) on huit rotors ont le meme 
principe. Le nombre de paires de rotors est ajonte par rapport an besoin de la 
poussee on de la charge utile. Pourtant, en ajontant des moteurs supplementaires, 
la batterie se decharge plus vite et l’autonomie du drone diminue. Le qnadrirotor 
est tres souvent utilise pour la simplicity de sa structure mecanique. Les couples de 
roulis et tangage sont obtenus en variant les vitesses de rotation de denx moteurs 
opposes. Les travaux de recherche sur le quadrirotor se concentrent generalement 
sur la stabilisation d’attitude et de la position. En effet, une fois que le quadrirotor 
est stabilise, on peut ajouter des fonctionnalites technologiques supplementaires en 
embarque pour augmenter son autonomie : par exemple, traitement d’images pour 
la perception de l’environnement et la navigation. 

Certains helicopteres utilisent un nombre impair de rotor, notamment trois rotors 
(voir figure 1.6). En fait, un des trois rotors est inclinable, ce qui va creer une petite 
force pour generer un couple contre la rotation a cause du couple resistant de ce 
rotor. L’inclination est controlee par un moteur servo. Cette structure mecanique 
est similaire au plateau cyclique dans la configuration d’helicoptere classique mais 
seul un degre de liberte est controle. L’inconvenient de ce type de configuration 
est son asymetrie qui peut poser le probleme pour l’equilibrage. L’avantage de la 
configuration a trois rotors par rapport a l’helicoptere classique est la simplicity 
de cette configuration par rapport a la structure du plateau cyclique et les helices 
utilisees a pas fixe. Puisque la configuration n’utilise que trois moteurs, le trirotor 
est plus economique que le quadrirotor. Ainsi, V autonomie est augmentee mais au 
detriment de la charge utile. 

Les travaux de l’helicoptere ayant deux rotors places en tandem se trouvent dans 
[35] et [36]. Sa configuration est montree dans la figure 1.7. Chaque rotor est 
capable de s’incliner suivant deux degres de liberte grace a un systeme electrome- 
canique complexe, des servomoteurs et des bras de transmission, ce qui permet 
de produire les couples de commande d’attitude. Le deplacement de translation 
horizontale se fait en inclinant les rotors vers l’avant. La position des rotors en- 
gendre une composante de poussee qui fait avancer le drone en translation. Les 
mouvement de tangage et lacet sont obtenus en basculant deux moteurs dans des 
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Figure 1.5 - Hexarotor 



Figure 1.6 — Trirotor et son troisieme rotor 


directions inverses. Le mouvement de roulis est du a la difference de poussee des 
moteurs. 


En general, l’helicoptere est un systeme sous-actionne, car il manque deux action- 
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Figure 1.7 — Prototype de birotors pivotants et son systeme de basculant des moteurs 

neurs pour le deplacement lateral. Dans [37] et [38], une configuration d’helicoptere 
a huit rotors est presentee. Cette configuration possede quatre rotors verticaux 
comme le quadrirotor classique mais il possede aussi quatre rotors supplementaires 
places a l’horizontal pour controler le deplacement lateral. Cet octarotor permet de 
decoupler la dynamique de rotation et la dynamique de translation. L’inconvenient 
principal de cette plateforme est son poids, du a l’ajout des moteurs lateraux, il 
reduit en fait significativement son autonomic et sa charge utile. Le prototype ne 
peut en effet voler que quelques minutes; la plupart des essais ont d’ailleurs ete 
faits avec une alimentation de secteur et non une batterie. L’autre inconvenient est 
le mauvais placement des moteurs lateraux. En effet, ils sont places a proximite des 
rotors principaux, ce qui genere une perturbation du flux d’air. Ainsi, le flux d’air 
genere par un rotor lateral augmente la portance du rotor principal, provoquant 
un desequilibre au niveau des couples de tangage ou de roulis, les deux dynamiques 
ne sont done pas tout a fait decoup lees. 

En conclusion, en comparant les configurations presentees ci-dessus, nous pouvons 
constater que le meilleur prototype du type a voilure tournante est le quadrirotor. 
En fait, cette configuration a le meilleur compromis entre le poids, la charge utile, 
le cout, le niveau de complexity de la structure mecanique et de la commande. 
C’est la raison pour laquelle le quadrirotor est le plus utilise non seulement dans 
la communaute scientihque mais aussi pour les loisirs. 


1.3.3 Les drones hybrides 

Ce sujet de these s’inscrit dans le cadre des recherches visant a obtenir des ve- 
hicules aeriens miniatures ayant a la fois les performances des drones a voilure 
hxe (avions) pour le vol horizontal telle que l’avance rapide, la longue portee et 
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Figure 1.8 - Octarotor 


une endurance superieure, et la manoeuvrability des drones a voilure tournante 
(helicopteres) pour le vol stationnaire telle que l’avance lente, le decollage et l’at- 
terrissage vertical (voir un exemple dans la figure 1.9). Ce type de recherche a 
suscite beaucoup d’interet dans le passe car ce type d’appareil ne necessite pas 
de piste de decollage et sa capacity de vol stationnaire le rend tres utile pour des 
missions de surveillance aerienne. Le travail de these s’insere aussi dans le cadre 
d’un contrat de recherche exploratoire innovante (REI) avec la DGA sur 1’etude 
et la conception d’un avion miniature autonome a decollage vertical. L’objectif du 
projet est en effet de concevoir et de realiser un prototype d’un systeme de drone 
pouvant effectuer des decollages/atterrissages verticaux de maniere autonome, puis 
de realiser une transition autonome vers un vol d’avancement rapide. Ce type de 
vehicule va etre aborde en details dans la section 1.4. 



Figure 1.9 — Le drone hybride qui combine les avantages de l’helicoptere et de l’avion 
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1.4 Aeronef a decollage et atterrissage verticaux 

Un aeronef a decollage et atterrissage verticaux (dont I’acronyme est ADAV, en 
anglais Vertical Take-off and Landing aircraft ou VTOL) est un aeronef a voilure 
fixe congu pour s’affranchir des pistes qui lui sont normalement necessaires pour le 
decollage et 1’atterrissage [75]. Les ADAV ont ete initialement developpes par les 
militaires, qui craignaient que leurs avions ne soient cloues au sol apres un bombar- 
dement de la piste. Ce type d’avion peut en effet etre utilise avec une infrastructure 
reduite au sol, tout comme un helicoptere. Dans le cadre naval, il est possible de 
construire des navires d’emport qui n’auront pas besoin d’une piste d’appontage 
aussi longue que sur un porte-avions classique (cas des porte-aeronefs). Dans le 
domaine civil, l’interet est surtout de combiner les avantages d’un helicoptere a 
ceux d’un avion, permettant ainsi le decollage depuis une plate-forme reduite et le 
vol stationnaire, avec une vitesse de croisiere et une charge utile superieure. 

Un defaut majeur des ADAV est d’avoir une consonnnation de carburant tres 
elevee dans les phases de decollage et d’atterrissage, d’ou une autonomie et une 
charge utile limitees, ce qui les rend non viables comme aeronefs civils et limite leur 
interet comme aeronefs militaires. Ces defauts ont ete partiellement resolus avec 
le concept de STOVL (Short Take Off Vertical Landing ), ou le decollage vertical 
est remplace par un decollage court. 


Techniques utilisees 

Plusieurs voies ont ete explorees pour permettre de raccourcir, voire de supprimer 
la distance parcourue horizontalement [75] : 

• L’aeronef decode a la verticale comme une fusee, bascule a l’horizontale pour 
le vol normal, puis rebascule a la verticale pour l’atterrissage (systeme designe 
"tail-sitter", par exemple : Snecma C450 Coleoptere, Convair XFY-1 Pongo, 
X-13 Vertijet) 

• L’aeronef reste a l’horizontale pendant toutes les phases de vol, ce sont uni- 
quement ses moteurs qui basculent de la verticale pour le decollage/atterrissage 
a l’horizontale pour le vol normal (systeme designe "tilt-rotor" dans le cas 
d’un moteur a helice et "tilt-jet" avec un moteur a reaction ou encore "tilt¬ 
wing" si c’est l’aile entiere qui bascule, par exemple : Boeing V-22 Osprey 
(figure 1.12), EWR VJ 101) 

• L’aeronef utilise la combinaison de rotors pour la sustentation et pour la 
translation, les moteurs entrainant l’un ou l’autre grace a un systeme de 
crabotage progressif, par exemple : Kamov Ka-22 (figure 1.13). 



1.4 Aeronef a decollage et atterrissage verticaux 
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Figure 1.10 - Le X-13 Vertijet en vol vertical 



1 tJiKH 


Figure 1.11 - Le X-13 Vertijet en vol horizontal 


• L’aeronef utilise deux groupes de moteurs differents suivant la phase de vol : 
des moteurs dit de sustentation places verticalement pour le decollage et 
1’atterrissage, et des moteurs de propulsion classiques places horizontalement 
pour le vol normal (systeme designe "flat-riser", par exemple : Dassault Mi¬ 
rage IIIV) 

• L’aeronef reste a l’horizontale, il dispose d’un senl reactenr fixe dont l’orien- 
tation du flux est variee (systeme de deflexion de la poussee, par exemple : 
Hawker Siddeley Harrier) 

• Une combinaison des denx dernieres solutions, avec des moteurs de sustenta- 
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Figure 1.12 - Boeing V-22 Osprey 



Figure 1.13 - Kamov Ka-22 


tion a l’avant et un moteur de propulsion avec deflexion de poussee a l’arriere 
(par exemple : Yakovlev Yak-38 Forger, Lockheed F-35 Lightning II) 



Figure 1.14 - STOVL Lockheed F-35B 


On appelle "transition" le moment pendant lequel l’ADAV passe du vol vertical 
au vol horizontal a Tissue du decollage, ou du vol horizontal au vol vertical avant 
l’atterrissage. La premiere transition est tres delicate car la portance doit etre creee 
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au moment ou la vitesse horizontale de l’avion est suffisante. Le basculement se 
fait done de fagon progressive. 


1.5 L’etat de Part 

Cette section donne une breve presentation de drones aerodynamiques qui sont 
capables de realiser les transitions du vol vertical au vol horizontal et inversement. 


1.5.1 T-wing 

Le T-Wing a ete construit par PUniversite de Sydney, en Australie. Ce vehicule 
est un drone a decollage et atterrissage vertical qui est capable d’effectuer le vol 
horizontal et le vol vertical (voir figure 1.15). Sa masse au decollage (avec le plein 
d’essence pour le moteur thermique) est de 30 kg, sa plus grande dimension, la 
hauteur, est de 1,50 m. II possede deux rotors contrarotatifs positionnes en tandem. 
Pendant la phase du vol vertical a une vitesse reduite, le T-Wing utilise la surface 
immergee dans le flux d’air cree par Phelice pour effectuer le controle. Au niveau 
de la commande, le vehicule utilise une combinaison entre la commande classique 
LQR [39]. II utilise deux systemes de coordonnees pour presenter l’attitude : un 
pour le vol vertical et un pour le vol horizontal. En effet, il faut changer le systeme 
de coordonnees lorsqu’il effectue la transition entre les modes de vol horizontal et 
vertical. 



Figure 1.15 - T-wing UAV 
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1.5.2 V-bat 

La societe MLB congoit et developpe une grande variete de projets sur les drones. 
Un de ceux-ci concerne un drone a decollage et atterrissage vertical appele V-Bat. 
Le 1 7-Bat est un type de "tail-sitter" qui peut decoder et atterrir verticalement ainsi 
qu’effectuer un vol stationnaire autonome, ce qui permet d’eliminer la necessity 
d’une piste ou d’une catapulte. En plus, il peut faire une transition doucement a 
l’horizontale pour une vitesse de croisiere jusqu’a 100 m/h (voir figure 1.16) 



Figure 1.16 - V-bat UAV 


1.5.3 Monoplane 

Le monoplane [14] est un aeronef a voilure fixe autonome disponible dans le com¬ 
merce avec la capacite de voler stationnairement, de faire la transition avec le vol 
de croisiere et de se percher sur une plate-forme d’atterrissage vertical. Ces ma¬ 
noeuvres ont ete demontrees dans l’environnement d’un espace tres limite : Real¬ 
time indoor Autonomous Vehicle test ENvironment (RAVEN) a MIT. Le RAVEN 
permet d’estimer tres precisement la position et la vitesse de drone dans l’espace 
grace a un systeme de capture de mouvement par cameras ( VICON ). L’attitude et 
les vitesses de rotation sont estimees par une centrale inertielle a bord. Pourtant, 
ce type de l’avion ne permet pas de decoder et d’atterrir verticalement, il faut uti- 
liser une plate-forme d’atterrissage vertical sur laquelle le drone se perche (figure 
1.17). Les auteurs dans [14] utilisent des lois de commande par retour d’etat pour 
la stabilisation et le vol horizontal a partir de systemes lineaires invariants dans 
le temps. Notons qu’ils ont acces a l’etat complet du systeme grace au systeme 
(VICON). Pour augmenter la precision du maintien en position, des integrateurs 
sont ajoutes dans le processus de commande. Les gains de controle sont optimises 
par des techniques Lineaires Quadratiques (LQ) afin de s’assurer du maintien en 
altitude du vehicule et de realiser la phase de transition.Les transitions entre les 
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modes stationnaire et horizontal sont realisees en boucle ouverte excepte pour les 
ailerons. 



Figure 1.17 - Le Monoplane dans RAVEN de MIT 


1.5.4 Vertigo 

Le Vertigo est un petit drone qui est dessine et congu par le Centre d’aeronau- 
tique et d’espace de PISAE (Institut superieur de l’aeronautique et de l’espace, 
Toulouse, France). Ce vehicule a la capacite de realiser des vols stationnaires et 
liorizontaux. Dans la configuration du Vertigo, deux helices coaxiales sont utilisees 
connne systeme de propulsion. Les surfaces de controle (ailerons et gouvernail) 
sont utilisees pour controler Pangle de roulis, de tangage et de lacet. Le poids de 
Vertigo est environ 1600g et son autonomie peut atteindre environ 12min (figure 
1.18). L’observation du comportement dynamique du drone a conduit au choix de 
deux structures de pilotage pour couvrir toute Penveloppe de vol; une structure 
de pilotage a gains statiques auto-sequences de type retour d’etats longitudinaux 
et lateraux pour le vol d’avancement; et une structure hierarchique axe par axe 
pour le vol stationnaire et pseudo stationnaire, proche de ce qui peut etre utilise 
pour le pilotage des helicopteres. Des outils d’analyse de robustesse (tel que la p- 
analyse) sont utilises pour prouver la stability de la boucle fermee tout au long des 
trajectoires en vol d’avancement [65]. 
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Figure 1.18 - Vertigo UAV 


1.5.5 Vision’Air 

Le Centre d’aeronautique et d’espace de 1’ISAE a aussi construit un autre drone 
a decollage et atterrissage appele Vision’Air. Le Vision’Air est capable de realiser 
la manoeuvre de transition entre le vol stationnaire et le vol horizontal. Le niveau 
d’autonomie permet a ce vehicule de realiser des missions telles que l’observation 
de l’environnement urbain (voir figure 1.19). 



Figure 1.19 - Vision’air UAV 
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1.5.6 C-Plane 

Le C-Plane [63] (voir figure 1.20) est la version precedente de l’avion convertible 
developpe dans la these. Le C-Plane a ete congu totalement au laboratoire par Oc¬ 
tavio Garcia. C’est un avion birotors contrarotatifs qui fonctionne bien en mode 
de vol horizontal comme un avion conventionnel. Pourtant en mode vertical, le 
flux d’air etroit de deux moteurs ne permet pas de creer la force aerodynamique 
necessaire aux ailerons pour stabiliser le mouvement de tangage. Pour resoudre le 
probleme de la stabilisation de Pangle de tangage, les moteurs sont places en tan¬ 
dem dans la nouvelle configuration (figure 2.1). Une loi de commande en utilisant 
des saturations separees des etats a ete proposee pour la stabilisation d’attitude 
et une approche de suivi de trajectoire pour le systeme linearise a ete developpee 
pour la transition. 



Figure 1.20 - Le C-Plane 


1.5.7 Drone tri rotors Panther 

Le Panther de PIAI (Israel Aerospace Industries ) (voir figure 1.21) possede une 
configuration a trois rotors electriques : deux rotors contrarotatifs situes sur les 
ailes basculant a l’avant pour le vol d’avancement, tandis que le rotor vertical 
fixe a l’arriere, uniquement actionne en mode helicoptere permet d’augmenter la 
stability et la poussee de l’appareil en phase verticale. Avec une masse d’environ 
65kg et une envergure de 6m, le drone affiche un rayon d’action de 60km, et peut 
voler a 3km d’altitude et mener une mission de surveillance de 6h au-dessus d’une 
zone d’interet. Un systeme de commande automatique controle la transition entre 
le decollage vertical et le vol d’avancement; le decollage et Patterrissage se font en 
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automatique, sur un simple "click" de l’operateur. Des prototypes du Panther ont 
deja accompli des essais et devaient etre operationnels en 2011. 



Figure 1.21 - Le Panther de 1’IAI [74] 


1.5.8 Drone quadri rotor 

Ce type de drone (figure 1.22) utilise quatre rotors places par paires devant et 
derriere l’aile. Chaque paire est montee sur un axe pivotant parallele a l’axe de 
tangage du drone. Un mecanisme de basculement base sur un actionneur de type 
servo moteur et un renvoi par biellette per met a I’ensemble des rotors de basculer 
simultanement afin de passer du mode helicoptere au mode avion. Ce type de sys- 
teme a deja ete mis en oeuvre par des modelistes. La commande etant entierement 
manuelle, ils utilisent deux gyroscopes places a 90° afin d’assurer la stability en 
mode helicoptere. 



MUfl 


Figure 1.22 - Drone quadri rotor 
























1.6 Organisation de la these 


20 


1.6 Organisation de la these 

Cette these est organisee de la fagon suivante. D’abord un etat de Part sur les 
diflerentes configurations de drones a ete decrite dans ce chapitre. Puis, les proto¬ 
types experimentaux choisis et developpes, la structure electronique et le systeme 
informatique embarque sont presentes dans le chapitre 2. Dans le chapitre 3, les 
modeles dynamique, aerodynamique et mathematique des drones sont analyses en 
details. La chapitre 4 traite des lois de commande. En se basant sur un modele 
du corps rigide dans l’espace tridimensionnel, des lois de commande non-lineaire 
sont etudiees. Une loi de commande extremement simple stabilisant la position 
et l’attitude d’un corps rigide est proposee. Cette loi de commande est ensuite 
appliquee pour commander le quadrirotor, l’avion convertible et havion basculant 
tri-rotor. Le chapitre 5 presente les resultats experimentaux. Les tests sont realises 
a l’interieur et aussi a l’exterieur. Enfin, les conclusions sur ces travaux de these 
sont donnees dans le chapitre 6. Les perspectives ainsi que les travaux futurs y 
sont discutes. 



Chapitre 2 


Presentation des plateformes 


Ce chapitre presente les plateformes developpees dans le cadre de cette these. II 
s’agit de denx prototypes : un avion convertible et un avion basculant tri-rotor. 
L’architecture electronique embarquee est decrite. Elle est composee des capteurs 
utilises pour obtenir l’attitude et l’altitude, des actionneurs du drone (moteurs 
sans balais et servomoteurs), du microcontroleur et de la station sol. Notons que 
ces elements peuvent etre utilises pour d’autres configurations de drones et/on de 
robots mobiles. 


2.1 Les prototypes 


L’objectif de la conception est d’avoir des drones legers, petits et qui combinent les 
deux modes de fonctionnement. Par le passe, plusieurs configurations ont ete de¬ 
veloppees au laboratoire (monorotor basculant et contraroratif, birotor) et testees 
[62], [63]. Les analyses de ces travaux ont permis de proposer des solutions ame- 
liorees par rapport aux configurations anciennes. L’avion convertible est inspire de 
havion basculant T-Wing (figure 1.15) et le drone Vision’Air (figure 1.19) pour 
leurs avantages et performances respectives. En effet, le couple de lacet est mieux 
controle avec une solution birotor en tandem (T-Wing) et nous nous sonnnes inspi¬ 
res de la structure mecanique et de l’emplacement de l’electronique embarquee du 
Vision’Air. L’avion basculant tri-rotor est inspire du drone Panther (figure 1.21). 



2.1 Les prototypes 
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2.1.1 Avion Convertible 

L’avion convertible est un avion du type "Tail-sitter" (voir figure 2.1). II a ete 
congu et fabrique completement au laboratoire. Dans la configuration, deux rotors 
contrarotatifs sont places cote a cote dans la partie superieure du vehicule et dans 
1’aile pour diminuer les risques lors des accidents imprevus. L’avion a ete bien ca¬ 
libre pour que le centre de gravite coincide avec le centre aerodynamique de 1’aile 
et ils se situent au quart du bord d’attaque de l’aile. Son poids est environ 900g et 
il utilise une batterie LiPo 2200mAh. La configuration choisie est tres performante 
en vol stationnaire puisqu’elle atteint environ 7 minutes en vol vertical. 

Dans le vol vertical, l’altitude est controlee par la poussee collective des deux 
moteurs electriques DC sans balais places en tandem et tournant chacun dans un 
sens different. Le mouvement vertical de roulis est realise par la difference de pous¬ 
see entre les moteurs. Le mouvement de tangage est controle par le deplacement 
des ailerons dans le meme sens, le mouvement vertical de lacet est controle par le 
deplacement des ailerons differentiellement. Les ailerons sont immerges dans le flux 
d’air des deux helices pour augmenter l’efficacite dans les phases du vol vertical. En 
vol horizontal, le mouvement de roulis est controle en deplagant differentiellement 
les ailerons et le mouvement horizontal de lacet est controle par le deplacement du 
gouvernail. Les vitesses de rotation des deux moteurs sont egales dans les phases 
de transition et de vol horizontal. 



Figure 2.1 - Avion Convertible 
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2.1.2 Avion basculant tri-rotor 

L’avion basculant tri-rotor a aussi ete construit et entierement developpe au labo- 
ratoire (voir figure 2.2). II a un poids total d’environ 1250g et utilise une batterie 
LiPo 2200mAh. Dans le vol vertical, 1’altitude est controlee par la poussee col¬ 
lective des trois moteurs electriques. Le mouvement de roulis est obtenu par la 
difference de poussee entre les deux moteurs avant : la rotation du moteur gauche 
dans le sens des aiguilles d’une montre alors que le moteur droit dans le sens 
contraire des aiguilles d’une montre. Le moteur arriere compense la force generee 
par les deux moteurs avant pour stabiliser Tangle de tangage quand le drone vole 
en mode vertical (decollage, atterrissage et vol stationnaire). 

Apres le decollage, les deux moteurs avant basculent vers l’avant grace a un sys- 
teme electromecanique de basculement. Le moteur arriere diminue sa vitesse pro- 
portionnellement a Tangle de basculement des deux moteurs avant. Quand Tangle 
de basculement atteint zero degre, c.-a-d. Taxe des deux moteurs avant coincide 
avec Taxe du fuselage, le moteur arriere s’arrete totalement. L’avion vole alors lio- 
rizontalement en croisiere comme un avion conventionnel avec une grande vitesse 
vers la cible. Lorsqu’il arrive a la cible, les deux moteurs avant basculent vers la 
verticale et le drone revient en mode lielicoptere, puis atterrit verticalement. Les 
vitesses de rotation des deux moteurs avant sont egales dans les phases de transi¬ 
tion et de vol horizontal. 



Figure 2.2 - Avion basculant tri-rotor 


Le tri rotor offre le meilleur compromis poids/stabilite, pourtant il rencontre un 
probleme pour stabiliser Tangle de lacet a cause du couple resistant du moteur 
arriere. Deux solutions ont ete proposees : soit le moteur est incline et le couple 
de reaction est annule par Tinclinaison de celui-ci, ce qui engendrera un vecteur de 
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poussee permettant de reequilibrer le lacet; soit deux helices contrarotatives sont 
montees sur ce moteur arriere par le biais d’une boite de transmission. Cependant, 
ces methodes augmenteront le poids ou la complexity de la structure mecanique. Si 
le probleme de la connnande de l’angle de lacet devient un critere important dans 
les phases de decollage et d’atterrissage, une des deux solutions sera appliquee. 


2.2 Capteurs utilises 

2.2.1 Capteurs utilises pour l’estimation de l’attitude 

Aujourd’hui, l’avancee technologique dans la construction et le developpement de 
capteurs a permis d’aborder les problemes d’estimation d’attitude dans differents 
domaines applicatifs. Les nouveaux composants integres sur silicium, les MEMS 
(Micro Electro Mecanical System ), permettent de realiser en grande quantity des 
capteurs performants a des dimensions micrometriques. Parmi la grande diversity 
de capteurs disponibles en technologie MEMS, nous pouvons citer les capteurs 
inertiels constitues par les accelerometres et les gyrometres. Mais aussi nous pou¬ 
vons trouver des capteurs qui permettent la mesure du champ magnetique ter- 
restre, les magnetometres. Le point fort de ce type de capteurs est leur taille, leur 
faible consonnnation energetique et leur cout. Par consequent, ils sont parfaitement 
adaptes pour les applications embarquees grand public. 


Capteurs de vitesse angulaire - Gyrometres 

Un gyrometre ou gyroscope vibrant est un appareil permettant d’effectuer une 
mesure de la vitesse de rotation. Les gyrometres utilisent plusieurs elements vi- 
brants qui produisent des signaux sinusoidaux identiques et d’amplitude constante 
lorsque le capteur est au repos. Lorsque celui-ci est soumis a un mouvement de ro¬ 
tation, les forces de Coriolis induisent une variation d’amplitude de ces signaux. Le 
fonctionnement du gyrometre vibrant a effet Coriolis est le suivant : il s’agit d’un 
diapason dont on entretient la resonance du mode de flexion a l’aide de 1’effet pie- 
zoelectrique indirect (utilisation d’un champ electrique pour creer des contraintes). 
Lorsque l’ensemble du capteur est soumis a une rotation, 1’effet Coriolis cree un 
mouvement orthogonal d’amplitude proportionnelle (mode de torsion) que l’on 
detecte par 1’effet piezoelectrique direct (apparition de charges sous 1’effet d’une 
contrainte). 
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In dotection mod© 


Figure 2.3 — Le principe de function du gyrometre 


Le signal de sortie du gyrometre est donne par : 

u G = u + v + rj G (2.1) 

ou c o est la vitesse reelle de rotation, v est une fonction lentement variable dans le 
temps appele le biais gyrometrique et r/c est un bruit blanc lie a la resolution du 
capteur. 

Lorsqu’il n’y a pas de rotation, le signal de sortie du gyrometre est compose d’un 
bruit blanc et d’une fonction lentement variable dans le temps. La dynamique de 
cette fonction peut etre modelisee par un processus de Gauss-Markov [43] decrit 
par l’equation : 

1 , . 

v — —v + 7) v (2.2) 

r 


Accelerometre 

Un accelerometre est un capteur qui, fixe a un mobile, permet de mesurer l’acce- 
leration. Les principaux parametres d’un accelerometre sont : l’etendue de mesure 
generalement donnee en g (l’acceleration de la gravite standard au niveau de la 
mer pour une latitude de 45°, g = 9, 80665m/s 2 ), la bande passante (en Hz), la 
precision, la sensibilite ( mV/g ) la gamme de temperature d’utilisation (en °C), 
la masse du capteur, la sensibilite transversale, le nombre d’axes (1 a 3 axes), la 
construction mecanique et l’electronique integree. La principale source de bruit 
mecanique presente dans un accelerometre est due au mouvement Brownien des 
molecules de gaz autour de la masse d’epreuve et au mouvement Brownien du 
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ressort. Le bruit resultant “TNEA” ( the Total Noise Equivalent Acceleration) est 
donne en [m/s 2 /y/Hz\. 

Le signal de sortie d’un accelerometre est decrit comme l’addition d’une acce¬ 
leration a, de la gravite g et d’un terme de bruit tja, considere blanc et gaussien : 

s b A = a- g + r) A (2.3) 

Si trois accelerometres sont montes en triaxe orthogonal dans un corps rigide, avec 
leurs axes sensibles qui coincident avec les principanx axes d’inertie dn corps rigide, 
la sortie des accelerometres pent etre ecrite comme le vectenr de mesnres : 

s a = R T s{ + f} A (2.4) 

avec R est la matrice de rotation entre les coordonnees inertielles et les coor¬ 
donnees attachees au corps rigide, s A — (a + g), ou g — [ 0 0 —g ] et a 6 l 3 
represented respectivement le vecteur de gravite et deceleration propre dn corps, 
donnes dans le systeme de coordonnees inertiel E^, g est la constante gravitation- 
nelle et tJa G K 3 est le vectenr de bruits supposes blancs et gaussiens. 

Si 1’acceleration a est faible par rapport a la gravite, ||a|| <C ||p||, les mesnres 
des accelerometres donnent principalement une projection dn vectenr de gravite 
dans le systeme de coordonnees inertiel E 6 . En consequence s A ~ge t : 

s!4 = rT 9 + Va (2.5) 

Dans ce cas, tja represente l’addition dn vectenr de bruit et des perturbations 
induites par des vibrations ou des accelerations parasites. Le vecteur g a ne peut 
pas alors etre considere comme un vecteur de bruits blancs et gaussiens. 


Magnetometre 

Parmi les instruments de navigation qui permettent de s’orienter par rapport a un 
systeme de reference absolu, nous trouvons le compas magnetique. Cet instrument 
utilise le champ magnetique terrestre comme reference, en consequence il permet 
de suivre un cap. Le champ magnetique terrestre interne trouve son origine au 
centre de la Terre ou les atonies du magma portes a haute temperature (5000°C') 
sont "depouilles" d’une partie de leurs electrons et sont done ionises. Ce noyau 
liquide, riche en fer et en nickel, est done conducteur. La rotation de la Terre 
met ce magma en mouvement ce qui constitue un gigantesque anneau de courant 
transformant la Terre en un immense solenoide avec un champ magnetique. En 
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Figure 2.4 — Champ magnetique de la Terre 


X 



(down) 


Figure 2.5 — Vecteur du champ magnetique 


surface, le champ magnetique varie d’un lieu a un autre (figure 2.4). 

Du point de vue theorique, le champ magnetique terrestre est oriente dans l’espace 
et presente deux composantes vectorielles en un point donne de la surface de la 
Terre : une composante d’induction magnetique He z dirigee verticalement et une 
composante horizontale Heh (voir hgure 2.5). La resultante represente l’intensite 
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du champ magnetique (ou le vecteur champ magnetique He). Aux poles magne- 
tiques la composante horizontale s’annule. L’angle forme par le champ magnetique 
He et le plan horizontal est appele l’inclinaison magnetique. Elle est positive vers 
le bas. Sa valeur augmente lorsqu’on se rapproche des poles et atteint 90° au pole 
magnetique nord (—90° au pole magnetique sud). L’angle forme entre la direction 
du pole nord geographique et le pole nord magnetique est appele declinaison (en¬ 
viron 11,5°). 


Le signal de sortie d’un magnetometre est decrit comme l’addition d’une projection 
du champ magnetique terrestre s M dans E 6 et un terme de bruit t)m considere blanc 
et gaussien. Si trois magnetometres sont montes en triaxe orthogonal dans un corps 
rigide, tels que leurs axes coincident avec les principaux axes d’inertie du corps 
rigide, la sortie du magnetometre triaxe dans le systeme de coordonnees mobile 
peut etre ecrit sous la forme d’un vecteur de mesures : 

= rT s{i + Vm (2.6) 


Un modele theorique du champ magnetique, tres proche de la realite consiste a 
considerer le vecteur de champ magnetique unitaire avec un angle d’inclinaison 


1 T 


I = 60° : s{, r = [ sm x 0 sm z ] = ^ 0 — ^ (lorsque la direction x des 

coordonnees inertielles est choisie suivant la direction du nord magnetique). 


Centrale Inertielle 

La Centrale Inertielle 3DM-GX3-25 (Microstrain) est utilisee pour obtenir l’at- 
titude. C’est un systeme miniature, de reference d’attitude et de cap (Attitude 
Heading Reference System (AHRS) en anglais). La 3DM-GX3-25 combine un ac- 
celerometre triaxial, un gyroscope triaxial, un magnetometre triaxial, un capteur 
de temperature et un processeur effectuant un algorithme de fusion de capteurs 
pour fournir l’orientation statique et dynamique. On peut egalement avoir acces 
aux mesures inertielles brutes. La Centrale Inertielle donne des mesures inertielles 
entierement calibrees : l’acceleration, la vitesse angulaire, le champ magnetique 
et des vecteurs des deltaAngle et deltaVelocity (correspondant respectivement a 
l’integration de la vitesse angulaire et de I’acceleration). Les donnees comprennent 
aussi des estimations de l’orientation calculee : l’angle de tangage, roulis, et le 
cap (l’angle de lacet) ou la matrice de rotation. Toutes les mesures sont totale- 
ment compensees en temperature et alignees mathematiquement par rapport a 
un repere orthogonal. L’architecture du systeme a ete soigneusement congue pour 
eliminer les sources d’erreurs communes telles que 1’hysteresis induite par les chan- 
gements de temperature et la sensibilite aux variations de tension d’alimentation. 
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La compensation permet d’utiliser des donnees de sortie avec une frequence faible 
en maintenant un taux d’echantillonnage interne rapide. 


O 

I 

mm I...mim...I 

0 20 40 

Figure 2.6 - La Centrale Inertielle 3DM-GX3™-25{ Microstrain) 



Quelques caracteristiques et avantages de la Centrale Inertielle 3DM-GX3-25 : 

• Centrale Inertielle tres petite et tres legere, seulement 11.5g. 

• Totalement compensee en temperature entre —40°C' to 75°C 

• Etalonnee pour les non alignements, pour l’erreur liee a la pesanteur, et la 
non linearite des gyrometres. 

• Adaptee pour la navigation avec des vibrations, des echantillonnages a 30 
kHz, et les donnees calculees a une frequence de 1 kHz 

• Disponible avec les interfaces de communication serie RS-232, USB 2.0 et 
TTL ( Transistor- Transistor Logic) 

• Frequence des donnees ajustable par l’utilisateur : jusqu’a 1 kHz 

• Les sorties : les angles d’Euler, la matrice de rotation, le vecteur de qua¬ 
ternion, les deltaAngle et les delta Velocity, les accelerations, les rotations 
angulaires et le champ magnetique... 


2.2.2 Capteurs utilises pour l’estimation de l’altitude 

Capteurs Ultrasons 

Le principe du temps de vol TOF (time of flight) est utilise dans les capteurs 
ultrasons. II correspond au temps de transmission des ondes et permet de mesurer 
la distance a un obstacle (voir figure 2.7). Cette methode est principalement ap- 
pliquee dans les appareils utilisant les ondes ultrasonores parce que la vitesse de 
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transmission des ondes ultrasonores dans Pair et dans d’autres materiaux est relati- 
vement lente (environ 343m/s dans Pair) et done nous pouvons mesurer la distance 
avec des erreurs mineures. Cette methode n’est pas appliquee dans des appareils 
utilisant des ondes electromagnetiques car la vitesse des ondes electromagnetiques 
est tres elevee (egale a la vitesse de la lumiere, 300000km/s). 



2 


Figure 2.7 - Le principe du temps de vol TOF 

Cette methode a un inconvenient : lorsque les ondes sont emises, nous ne pouvons 
pas gerer les ondes reflechies. Nous savons que la plupart de Penergie des ondes 
reflecliis se concentre aux rayons reflechies conformement au principe : Pangle d’at- 
taque est egal a l’angle de reflexion (figure 2.8). Pourtant, en dehors des capteurs 
ultrasons speciaux, les capteurs ultrasons generaux ont un grand angle ouvert, ce 
qui entraine des erreurs dans la mesure de la distance : 

• Le phenomene " forecasting " : Conformement au principe TOF, pour 
avoir une mesure correcte, le capteur ultrason doit etre place perpendiculaire- 
ment a la surface des objets. Pourtant la surface des objets n’est pas toujours 
plane, le rayon reflechi ne correspond pas a l’angle d’attaque. Les faisceaux 
des rayons reflecliis ont une energie de reflexion plus faible. Neanmoins, le 
capteur peut detecter et enregistrer ces signaux reflecliis. Par consequent, les 
mesures obtenues ne sont pas tout a fait correctes a cause de l’angle ouvert 
tres grand des capteurs ultrasons. De plus, lorsque l’angle ouvert est tres 
grand, non seulement la position de Pobstacle est detectee incorrectement, 
mais la distance est aussi inexacte (figure 2.9). Pourtant, ce type d’erreur est 
controlable, nous pouvons proposer des filtres pour compenser ces erreurs. 

• Le phenomene " crosstalk " : e’est un phenomene indesirable dans lequel 
un capteur regoit les signaux ultrasonores reflecliis ou directs emis par un 
autre capteur ; ou bien les ondes ultrasonores sont transmises, se reflechissent 
via plusieurs surfaces et retournent a ce capteur de maniere indesirable. Le 
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Figure 2.8 — La reflexion des ondes ultrasonores sur une surface 



Figure 2.9 — Le phenomene forecasting 


phenomene a lieu quand plusieurs robots fonctionnent dans un meme endroit 
et un capteur d’un robot regoit les signaux emis par un capteur d’un autre 
robot directement ou indirectement ou lorsqu’un capteur regoit les signaux 
emis par un autre capteur du meme robot apres la reflexion des ondes. Dans 
le premier cas, les erreurs sont aleatoires et impredictibles et sont done fa- 
cilement filtrees. Pourtant, pour les erreurs de crosstalk de deux capteurs 
attaches sur un meme robot, les resultats obtenus ne permettent pas de dis- 
tinguer quels sont les bruits, et quels sont les signaux emis de chacun des 
capteurs. 

Dans [44], une methode appelee EERUF (Error Eliminating Rapid Ultrasonic 
Firing ) est developpee ponr utiliser plusieurs capteurs ultrasons dans les 
applications de robot mobile. L’idee principale de cette methode est qu’en 
considerant les signaux obtenus comme des bruits, elle utilise un ordre de 
transmettre et de recevoir des ondes pour eliminer les erreurs inconvenables. 
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Figure 2.10 - Le phenomene crosstalk [44] 


La methode EERUF permet aux capteurs ultrasons d’emettre des ondes 
a des taux qui sont cinq a dix fois plus rapide que dans les applications 
classiques. Par ailleurs, la methode ERRUF a la capacite unique de rejeter 
les bruits et elle permet a plusieurs robots mobiles de collaborer dans le 
meme environnement, meme si leurs capteurs ultrasons fonctionnent aux 
memes frequences. 


Capteur Ultrason SRF10 : Le capteur ultrason SRF10 cliez Devantech est 
utilise dans notre systeme embarque pour la mesure de l’altitude. La communica¬ 
tion avec le SRF10 est via le bus I2C. Le faisceau d’onde du SRF10 est conique 


+&/ Power 
SDA 
SCL 

Do Not Connect 


Ov Ground 




Figure 2.11 — Le Capteur Ultrason SRF10 

avec la largeur etant une fonction de la surface des transducteurs et etant fixe. Le 
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faisceau des transducteurs utilises dans le SRF10 est montre dans la figure 2.12 
avec l’unite sur l’axe vertical dans le diagramme de faisceau en dB. 



Figure 2.12 - Le faisceau des transducteurs utilises dans le capteur ultrason SRF10 


2.3 Motorisation 


II existe differents types de moteurs : electriques, thermiques, pneumatiques, hy- 
drauliques... Seuls les moteurs electriques et thermiques sont utilises pour la pro¬ 
pulsion d’engins volants de type drone. Les moteurs thermiques sont trop encom- 
brants et trop lourds pour les drones que nous utilisons. Ils sont aussi plus puissants 
et done plus dangereux. Dans le cadre de cette these, nous ne nous interessons 
qu’aux moteurs electriques. Parmi les differents types de moteurs electriques a 
courant continu, le moteur sans balais, (Brushless Direct Current (BLDC) motor 
en anglais), a ete retenu pour la plateforme experimentale. Ces moteurs sont en 
effet de plus en plus utilises en modelisme et se trouvent facilement dans le com¬ 
merce. Leur electronique de commande (appelee driver) simplifie l’electronique de 
puissance a concevoir sur le drone. II est aussi possible de mesurer la vitesse de 
rotation du moteur sans acheter un capteur supplemental, ce qui est interessant 
pour asservir les moteurs en vitesse. 
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2.3.1 Moteur sans balais 

Le moteur BLDC est utilise pour les consonuuateurs et les applications indus- 
trielles grace a sa taille compacte, sa controlabilite et sa haute efficacite. En effet, 
le moteur sans balais est un moteur electrique de la categorie des moteurs syn- 
chrones, le rotor est const it ue d’un ou de plusieurs aimants permanents et d’un 
capteur de position rotorique (le capteur a effet Hall, le synchro-resolveur, le co- 
deur incremental par exemple). Vu de l’exterieur, il fonctionne en courant continu. 
Son appellation (de l’anglais Brushless ) vient du fait que ce type de moteur ne 
contient aucun collecteur tournant et done pas de balais. Par contre, un systeme 
electronique de commande doit assurer la commutation du courant dans les enrou- 
lements statoriques. Ce dispositif peut etre, soit integre au moteur pour les petites 
puissances, soit exterieur sous la forme d’un convertisseur de puissance (onduleur). 
Le role de Lensemble capteur plus electronique de commande est d’assurer l’auto- 
pilotage du moteur e’est-a-dire l’orthogonalite du flux magnetique rotorique par 
rapport au flux statorique, role autrefois devolu a l’ensemble balais-collecteur sur 
une machine a courant continu. 

Les moteurs reduits sans balais sont tres utilises en modelisme dans des modeles 
d’avions, d’helicopteres. Ils sont moins bruyants que les moteurs avec balais, ce- 
pendant leur cout est encore superieur. Dans les moteurs sans balais utilises pour 
les modeles reduits radioconnnandes, le stator (bobine) est a l’interieur et le rotor 
(comportant les aimants) a l’exterieur (voir hgure 2.13). Le moteur sans balais est 
generalement utilise avec un ou plusieurs capteurs de position du rotor car l’exci- 
tation electrique doit etre synchrone avec la position du rotor. Pour des raisons de 
cout, de habilite, de conditionnement mecanique et surtout si le rotor tourne im- 
merge dans un fluide, il est souhaitable de faire marcher le moteur sans capteur de 
position, ce qui est connu comme le fonctionnement sans capteur ( sensorless ). En 
effet, certains circuits de commande de moteurs brushless utilisent les bobines non 
alimentees a un instant donne pour mesurer la force contre-electromotrice (FCEM) 
{Back-Electromotive Force (BEMF) en anglais) et pour deduire le moment pour 
declencher la commutation. Cette solution perrnet d’eviter l’emploi de capteurs de 
position a effet Hall, et done de reduire le prix du moteur brushless. 


Drivers YGE30i 

Les drivers sont l’electronique de commande qui permettent de controler les mo¬ 
teurs. Certains drivers mesurent la vitesse du moteur et cette valeur de vitesse 
peut etre lue par le microcontroleur principal. C’est le cas des drivers utilisant une 
communication digitale (le protocole UART, I2C, SPI, CAN par exemple). A bin- 
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Figure 2.13 — Schema en coupe d’un moteur sans balais de faible puissance a rotor 
externe utilise en modelisme 


verse, les drivers utilisant une communication analogique (le signal de PWM) ne 
permettent pas de lire la vitesse de rotation, le driver. Dans le commerce, le driver 
YGE30i (figure 2.14) qui utilise la communication I2C nous permet d’acceder a 
cette valeur et aussi aux autres parametres, par exemple le voltage, le courant ou 
la temperature du driver. L’adresse des ces variables sont donnees dans le tableau 
2 . 1 . 



Figure 2.14 - Driver YGE30i 
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Table 2.1 — Les adresses des variables du driver YGE30i 


Adresse 

Variable 

0x00 

Octet 0 de la vitesse 

0x01 

Octet 1 de la vitesse 

0x02 

Octet 2 de la vitesse 

0x03 

Voltage 

0x04 

Courant 

0x05 

Temperature 

0x06 

Derniere erreur 

0x07 

Adresse I2C de l’esclave 


Commande de la vitesse de moteur 

Le temps entre deux passages de la force contre-electromotrice (FCEM) a zero 
est considere comme proportionnellement inverse a la vitesse de rotation. Nous 
l’appelons la valeur reelle de la vitesse de moteur. Nous pouvons constater que 
quand le rotor du moteur brushless tourne plus vite, ce temps est reduit parce que 
la frequence de passage au point de zero de la FCEM est plus elevee. Pour faciliter 
la commande de la vitesse du moteur, on utilise la relation suivante : 

la valeur nominale mesuree = 10000 — la valeur reelle (2-7) 

ou la valeur nominale mesuree est la valeur nominale entrant de la loi de com¬ 
mande PI de la vitesse angulaire, c’est aussi la valeur que les lois de commande 
principales fournissent; la valeur reelle est la valeur lue du microcontroleur du 
driver , c’est en fait le temps entre deux passages au point de zero de la FCEM. En 
utilisant la valeur nominale mesuree, la vitesse angulaire indiquee est maintenant 
directement proportionnelle a l’amplitude de la vitesse de rotation reelle du rotor. 

La figure 2.15 montre l’adaptation ouverte de la vitesse de moteur lorsque la 
consigne PWM varie de 0 a 120. Notez que sur la figure 2.15(a) et 2.15(b), la 
valeur PWM utilisee pour l’acceleration est plus elevee que la deceleration. C’est a 
cause de l’inertie de moteur. Ainsi, a l’aide de Matlab, un polynome approximant 
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au mieux (au sens des moindres carres) la relation entre la valeur PWM et la vi- 
tesse nominale de moteur est constrnite. Un ordre 5 a ete retenu (la figure 2.15(b)) 
et nous obtenons l’equation suivante : 

y = 0.0000038268x 5 —0.0014304x 4 + 0.21348x 3 —16.336x 2 +679.45x —4527.2 (2.8) 
ou x est la valeur PWM et y est la vitesse nominale du moteur 




(a) Sans l’helice (b) Avec l’helice 

Figure 2.15 — Adaptation des vitesses du moteur (boucle ouverte, la valeur PWM varie 
de 0 a 120) 

Lorsque l’operateur change manuellement le gaz, il donne la consigne de PWM, 
puis cette consigne PWM est convertie vers la valeur de la vitesse nominale desiree 
du moteur suite a l’equation (2.8). Le controleur PI de la vitesse du moteur va 
donner ensuite la valeur PWM au driver pour que le moteur tourne avec la vitesse 
desiree. Notons que Loperateur ne donne pas directement la consigne PWM au dri¬ 
ver , cette procedure est pour eviter le probleme d’avoir deux moteurs qui tournent 
avec des vitesses differentes malgre que la meme consigne PWM est fournie. 

Une loi de commande de type PI est utilisee pour commander la vitesse de rotation 
du moteur. Le controleur prend directement la vitesse nominale de rotation lue via 
le bus I2C en entree et donne la valeur PWM en sortie. La Figure 2.16 montre la 
reponse du moteur en boucle fermee, pour une consigne qui varie. Dans la figure 
2.16(c), les valeurs PWM fournies aux deux drivers sont legerement differentes pour 
commander la meme vitesse de rotation, c’est a cause de la difference entre les deux 
systemes de motorisation (driver + moteur). Pourtant, les valeurs nominates de 
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rotation de deux moteurs sont identiques car c’est la valeur calculee par le timer 
de microcontroleur de driver et ces microcontroleurs fonctionnent avec la meme 
frequence. 
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(a) Moteur 1 


(b) Moteur 2 



(c) Les valeurs PWM des deux drivers 

Figure 2.16 — Adaptation des vitesses de moteur et les valeurs PWM generees par la 
commande PI 


2.3.2 Servomoteur 

Un servomoteur est un moteur a courant continu, dont 1a, position est verifiee et 
corrigee en fonction de 1a, mesure. Alors, un servomoteur est un systeme motorise 
capable d’atteindre des positions predetermines, puis de les maintenir. Pour un 
ajustement precis de la position, le moteur et son reglage sont equipes d’un systeme 
de mesure qui determine la, position courante du moteur (par un potentiometre qui 
genere une tension variable proportionnelle a Tangle de l’axe de sortie). Lorsque le 








































39 


Presentation des plateformes 


moteur tourne, l’axe du servomoteur change de position, ce qui modifie la resistance 
du potentiometre. Un dispositif electronique d’asservissement va commander le 
moteur pour que la position de l’axe de sortie soit conforme a la consigne regue. 
Les servomoteurs servent a actionner les parties mobiles d’un avion comme : les 
ailerons, l’elevateur et le gouvernail. Nous envoyons done des signaux PWM par 
rapport aux angles de deflexion avec la frequence de 50Hz (voir la figure 2.17). 



Figure 2.17 — Le servomoteur et le signal PWM par rapport a l’angle de rotor 


2.4 Energie 

Etant donne que les moteurs retenus sont electriques, cette section portera sur 
les batteries electriques. II existe un grand nombre de technologie de batteries 
diflerentes, par exemple : 

• les batteries au plomb ( Pb ) : tres repandues dans les automobiles ou les 
motos 

• les batteries nickel-cadmium ( NiCd ) et nickel-metal-hydrure ( NiMh ) : 
utilisees pour les piles rechargeables 

• les batteries lithium-ion ( Lilon ) : utilisees dans les batteries plus recentes 
telle que celles des telephones portables, des ordinateurs, portables, les lec- 
teurs MP3, etc. 
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• les batteries lithium-polymere ( LiPo) : repandues dans le modelisme 


Une batterie est caracterisee par sa tension et sa charge. La tension est exprimee 
en volts et correspond a la tension nominale de la batterie. La charge ( C) cor¬ 
respond a la capacite de la batterie, elle est exprimee en mAh. C’est la quantite 
moyenne d’electricite qu’elle peut debiter en line heure. Une batterie est en fait 
un pack, compose de plusieurs elements (ou "cells" en anglais). Chaque element a 
une tension nominale dependant de la technologie. Par exemple un element NiCd 
a une tension nominale de 1,2V , un element Pb 2,1V et un element LiPo 3,7V. La 
tension de la batterie est egale a la somme de la tension de chacun de ses elements 
(s’ils sont montes en serie), sachant que les elements peuvent etre charges a une 
tension legerement superieure a la tension nominale. 

Seules les batteries ayant un rapport puissance/poids important nous interessent 
ici, les batteries de type Pb, NiCd et NiMh ne seront done pas traitees. En effet, 
les batteries au lithium proposent des capacites equivalentes aux batteries tradi- 
tionnelles mais pour le tiers du poids. En plus, les batteries LiPo ne presentent pas 
d’effet memoire. Cependant, leurs principaux inconvenients se situent au niveau 
de leur capacite de charge, de leur fragilite et de leur prix. En effet, le taux de 
charge de ces batteries est de 1 C ; done il faut toujours une heure pour la charger. 
Le taux de decharge maximum est limite a 15 fois de la capacite (ISC'). 


rrFu-TroniK 


2200 ,. 


Figure 2.18 - Batterie lithium-polymere (Pro-Tronik, LiPo 2200mAh) 


Les batteries LiPo sont fragiles mecaniquement et electriquement. En effet un choc 
sur la batterie peut provoquer un court circuit interne, et la batterie peut alors 
prendre feu (meme plus de 10 minutes apres le choc). La batterie peut aussi prendre 
feu apres une mauvaise charge. II faut done veiller a toujours la charger au bon 
voltage (4,2V par element) et a la bonne capacite. Lors de la charge d’une batterie, 
il est fortement conseille d’utiliser un equilibreur de charge. Celui-ci permet de 
charger les elements independamment en controlant leur tension. Il faut aussi ne 
pas trop decharger une batterie LiPo. Si un element tombe en dessous de 2,5V, il 
devient alors impossible de le charger et done inutilisable (rendant alors souvant 
le pack entier inutilisable). Les microcontroleurs dans les drivers pour les moteurs 
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sans balais permettent d’ailleurs de couper le moteur si la tension d’alimentation 
devient trop faible, pour eviter de decharger en dessous de ce seuil critique. 


2.5 Communication air/sol 

Le drone, bien qu’autonome, doit toujours garder une liaison avec le sol. Cela 
per met avant tout de pouvoir reprendre la main en cas de problemes : si le com- 
portement du drone devient instable, l’operateur au sol peut couper les gaz ou 
revenir en mode manuel pour eviter tout accident. Cette liaison unidirectionnelle 
est effectuee grace a une radio commande (figure 2.19) issue du modelisme. Cellc-ci 
est associee a un recepteur radio integre au drone. La liaison sans fil s’effectue a 
2,4 Ghz et a une portee d’environ 2km. En effet, les radios a 2,4 Ghz utilisent so it 
la technologie FHSS ( Frequency Hopping Spread Spectrum ) de saut de frequence, 
soit le DSSS ( Direct Sequence Spread Spectrum) d’etalement direct. Ce type de 
radio possede en general de quatre a six canaux et les donnees transmises sont 
analogiques. 

Le signal de sortie du recepteur est represente sur la figure 2.20, il s’agit d’un signal 
PPM (Pulse Position Modulation). La valeur de chacun des canaux est codee par 
la duree entre deux fronts montants successifs; cette duree etant comprise entre 
1 et 2ms. La duree des etats hauts est fixe, de l’ordre de 0,5ms. Un etat bas de 
plus de 5ms indique le debut d’une nouvelle frame. Normalement la duree totale 
d’une frame est de 20ms ou une frequence 50Hz. Afin de recuperer les valeurs des 
canaux, une simple interruption de type " input capture " sur le microcontroleur 
suffit. Cette interruption est en effet associee a un timer, et permet done de cal- 
culer le temps separant deux fronts montants successifs. La precision de la mesure 
depend alors de la precision du timer utilise. Les principaux inconvenients de ce 
mode de communication sont done son aspect unilateral et que la communication 
ne possede aucun outil pour controler et verifier les signaux regus. Cependant elle 
reste tres simple a mettre en place et offre un moyen de commande leger, la radio- 
commande etant facilement transportable. 

Pourtant, cette liaison unilateral est insuffisante car la radiocommande ne permet 
que d’envoyer des parametres de vol vers le sol, ce qui reste tres limite. Ainsi, une 
liaison numerique bilaterale par modem sans fil a ete mise en place. Associee a 
une interface homme machine sur ordinateur, elle permet d’envoyer les etats du 
drone a la station sol et l’operateur peut changer et envoyer les parametres du 
programme plus rapidement. Le modem retenu pour le drone est un XbeePro de 
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Throttle/rudder 
control stick 
(Mode 2) 

Throttle trim 
lever (Mode 2) 


Aileron/elevator 
control stick 
(Mode 2) 

Elevator trim lever 
(Mode 2) 

Aileron trim lever 

Charging jack 

On*off switch 


DATA INPUT lever 
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Neck strap hook 


Aileron, Elevator & Rudder 
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A:Flap switch/CH.6 
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Figure 2.19 — Le radiocommande et le recepteur radio (Futaba, T6EX) 


«- 



CHI 2.00V M 1.00ms 

Figure 2.20 — Une trame d’une radiocommande 5 canaux vue a l’oscilloscope 


chez Maxstream (voir figure 2.21). II possede un mode " transparent " : place dans ce 
mode, le modem gere automatiquement les transmissions et il suffit de le brancher 
au port serie d’un ordinateur ou d’un microcontroleur. La communication est alors 
tout de suite operationnelle. Cependant, ce mode " transparent " reste tres limite. 
En effet, lorsqu’une donnee du microcontroleur est fournie au modem emetteur, 
elle est mise dans un buffer d’entree. La donnee sera alors transmise au modem 
recepteur si Tune des conditions suivantes est satisfaite : 

• si le buffer est plein (sa taille etant parametrable) 

• ou si un delai (lui aussi parametrable) s’est ecoule sans que le modem emet- 
teur ne regoive de donnees du modem recepteur 

Une trame Xbee contenant les donnees du buffer est alors envoyee et le modem 
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Figure 2.21 - Modem XbeePro (Digi International®, XbeePro) 


recepteur envoie un ACK (acknowledgment) pour signaler qu’il a bien regu la 
trame. Si l’emetteur ne regoit pas cet ACK apres un certain temps, il reitere la 
transmission (jusqu’a six fois). Lorsque le recepteur regoit une trame, il la transmet 
via sa broche de sortie. Cependant, si la trame n’a pas pu etre regue elle est perdue. 
Dans ce cas, le modem n’avertit pas le microcontroleur (ou l’ordinateur) car en 
mode transparent aucune donnee autre qu’une donnee envoyee par l’emetteur ne 
sort du modem. Il apparait done necessaire de verifier les donnees regues grace 
a un checksum par exemple, car des donnees peuvent etre perdues. Par ailleurs, 
une trame Xbee perdue peut correspondre a plusieurs trames de la communication 
ordinateur-microcontroleur. De plus, etant donne que le modem renvoie les donnees 
non regues, il est possible que les trames n’arrivent pas dans l’ordre ou elles ont 
ete envoyees, voir figure 2.22. 



Figure 2.22 — Schema d’une communication Xbee. La trame Xbee n°2 a ete envoyee trois 
fois sans succes et est perdue. La trame Xbee n °4 n’a pas ete regue la premiere fois ; elle 
a ete regue a la secondee tentative mais entre temps le trame Xbee n °5 a ete regue 


Enfin, ce mode n’assure done pas qu’une donnee envoyee arrive a destination. Cela 
peut done etre problematique, car certaines trames doivent etre regues, tel que 
des ordres d’arret des moteurs en cas de problemes. En general, les trames sont 
perdues a cause de collisions. En effet, la liaison radio frequence n’est pas biclirec- 






















2.6 Microcontroleur 


44 


tionnelle. Un seul des deux modems peut transmettre a la fois, sinon il y a collision 
et la trame est perdue. Pour eviter cela, les modems, possedent des algorithmes 
anti-collisions plus ou moins evolues. Cependant, plus les frequences d’envoi sont 
elevees, plus de risque de collision et eleve. Ainsi, dans l’application de station sol 
environ 10% des paqnets sont perdus. 

Le modem possede cependant une autre mode, dit API (Application Program¬ 
ming Interface), beaucoup plus complet mais aussi plus complexe. Dans ce mode, 
l’utilisateur envoie les donnees au modem sous une syntaxe bien precise. Cela per- 
met entre antres de former une frame qui sera envoyee en totalite dans une senle 
frame X'Bee, et de recevoir YACK ou le non-ACK (NACK). En cas de NACK, le 
logiciel pent decider de renvoyer la frame si elle est importante. En plus, le proto¬ 
cole de communication inclut un checksum et donne pour chaque trame regue la 
puissance du signal. 


2.6 Microcontroleur 

Le microcontroleur est le "cerveau" du drone. C’est lui qui va recevoir les infor¬ 
mations des differents captenrs, il les filtre, il effectue les calculs lies aux lois de 
commandes et il pilote les actionnenrs. Pour choisir le microcontroleur de la pla- 
teforme, un cahier des charges englobant differentes architectures possibles a done 
ete dehni : 

• controle par radiocommande et/ou modem 

• moteurs sans balais avec le controle de vitesse 

• centrale inerticlle du commerce par port serie, ou capteurs inertiels (analo- 
gique ou numerique) 

• telemetres analogiques ou numeriques (pour mesurer la distance au sol ou 
aux obstacles) 

• possibility de communiquer avec une antre carte effectuant les calculs de 
vision embarqnee (ponr les developpements futurs) 

Ponr repondre a cela, les peripheriques snivants sont done necessaires : 

• 4 PWM pour les servo moteurs 

• 1 input capture pour la radiocommande 

• 1 port serie ponr le modem 

• 1 port serie ponr la centrale inertielle 
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• 3 entrees analogiques pour les mesures du voltage de la batterie et des cap- 
teurs analogiques (barometre, telemetre) 

• 1 bus SPI pour les capteurs numeriques (inertiels, barometre, telemetres) 

• 1 bus I2C pour les capteurs numeriques (inertiels, barometre, telemetres) 

• 1 bus I2C pour communiquer avec les drivers YGESOi des moteurs 

• 1 port de communication pour recevoir les donnees du systeme de vision 
embarque (a utiliser dans le futur, soit le port serie , le bus I2C, SPI ou 
CAN ) 

II faut noter que dans le cas des capteurs numeriques le bus de communication peut 
etre UART , SPI ou I2C, suivant le constructeur. Par ailleurs, l’existence d’outils 
de developpements gratuits et d’outils de programmation libres ont fait partie des 
criteres de selection. Le microcontroleur dsPIC33FJ256GP710 chez Microchip a 
ete choisi pour ses caracteristiques suivantes : 

• 16 bits, 80 MHz, 40 MIPS, 256 Ko de memoire flash , 30 Ko de memoire ram 

• 9 timers de resolution 16 bits, 4 timers de resolution 32 bits (pouvant servir 
aux PWM et aux input capture ) 

• 2 ports serie 

• 2 bus I2C 

• 2 bus SPI 

• 32 entrees analogiques {ADC, Analog to Digital Converter ) 

• 1 module de DSP (Digital Signal Processor ) 

• 8 canaux hardware DMA (Direct Memory Access) 

• 2 bus ECAN {Enhanced Controller Area Network ) 

En fait, dans le systeme electronique embarque le microcontroleur va effectuer pa- 
rallelement les travaux suivants : Lacquisition de donnees d’attitude de la Centrale 
Inertielle, la lecture des entrees de commande du recepteur radio, la realisation de 
l’algorithme de commande, la generation des signanx PWM correspondants pour 
commander les angles de deflexion des ailerons, de l’elevateur et du gouvernail, et 
la transmission sans fil des parametres de vol a la station sol. 


2.7 Carte electronique embarquee 


Une carte electronique (voir figure 2.23) a ete congue et realisee afin de reunir les 
differents elements vus precedennnent autour du microcontroleur. 
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Figure 2.23 - Carte electronique 


Les travaux de conception de la carte electronique ont ete realises soigneusement 
pour garantir une fonctionnement correcte et une taille non encombrante. Dans 
la conception, pour eviter le bruit electronique induit par les moteurs, la partie 
electronique du signal et de puissance ont ete totalement decouplees. Cela permet 
d’avoir notamment des mesures moins bruitees sur les capteurs inertiels (accelero- 
metres, gyrometres et magnetometres). 

Toutes les lois de commandes (y compris celles des moteurs) sont programmees 
dans ce microcontroleur dsPIC33FJ256GP710. II est relie au modem XbeePro et 
a la centrale inertielle 3DM-GX3-25 par les ports series, au recepteur radio, au 
capteur ultrasons et aux drivers YGESOi par le bus I2C et aux moteurs servos 
par les broches PWM. Les modules SPI et EC AN sont conserves pour des deve- 
loppements futurs. L’architecture complete est representee sur la figure 2.24. 
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Radio 2.4GHz 


CAN, SPI 


Pilote 


Modules 

Supplemenlaires 


ZigBee 2.4GHz, IEEE 802.15.4 


PPM: Pulse Position Modulation 

PWM: Pulse Width Modulation 

I2C: Inter Integrated Circuit 

UART: Universal Asynchronous Receiver Transmitter 

DMA: Direct Memory Access 

CAN : Controller Area Network 

SPI: Serial Peripheral Interface 



Figure 2.24 — Architecture electronique embarquee 


2.8 Station sol 


Dans la section 2.5, une liaison sans fil entre le drone et le sol est presentee, ce qui 
permet de transferer des donnees de vol au sol. Dans cette section, le logiciel de 
la station sol est presente brievement. La station sol a ete ecrite en Visual C++ 
(figure 2.25), elle permet d’afficher a l’aide de graphiques tous les etats du drone : 
les consignes de l’attitude et de l’altitude, les angles d’Euler, les vitesses angulaires, 
les vitesses de rotation des moteurs (ainsi que leurs consignes), l’altitude actuelle, 
les angles des ailerons, de l’elevateur et du gouvernail. Toutes ces donnees sont sau- 
vegardees et horodatees dans un fichier execute par le logiciel Matlab pour pouvoir 
afficher les graphes de chaque essai. La figure 2.25 affiche les graphes des donnees 
regues au sol. Les courbes en blanc represented les consignes tandis que celles en 
bleu sont les valeurs mesurees. On peut voir qu’il y a cinq ecrans d’affichage per- 
mettant de visualiser l’etat complet, les vitesses et les entrees de commande. Sur 
la figure, on peut voir les valeurs mesurees et les consignes de 1’altitude, du gaz, 
de Tangle de roulis, de tangage, de lacet et les valeurs PWM de chaque moteur. Le 
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bandeau inferieur contient les angles qui doivent etre observes pendant toutes les 
phases de vol. Un bouton "STOP" permet d’envoyer une commande au drone afin 
d’arreter les moteurs en cas de probleme. Les autres boutons servent a demarrer 
et arreter 1’ecriture des donnees, pour eviter de surcharger la memoire quand le 
stockage de donnees n’est pas necessaire. 



Figure 2.25 - Station sol 

La station sol permet aussi de changer rapidement les gains des lois de commandes 
embarquees, et de les sauvegarder dans la memoire non volatile du microcontroleur 
(la memoire flash). En effet, lors de Limplementation d’une loi de commande, il 
est toujours necessaire d’ajuster les gains manuellement en faisant des essais et en 
observant le comportement du drone entre chaque changement et chaque essai. Le 
fait de pouvoir les modifier instantanement apporte ainsi un grand avantage, cela 
permet d’ajuster et d’essayer differentes valeurs rapidement. Sans la station sol, il 
aurait fallu reprogrammer le microcontroleur pour changer les gains. 

La station sol a ete developpee dans le contexte de l’utilisation du modem XBee 
qui utilise le mode API (comme le programme embarque dans le drone). Elle nous 
permet de controler chaque frame envoyee et regue afin de savoir quelle frame 
est perdue pendant la transmission et le taux de la transmission reussi. Ce taux 
nous indique la qualite de la communication pour changer les reglages (frequence 
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d’emission et le nombre maximum d’octets transmises par unite de temps) afin 
d’arriver a une meilleure performance. De plus, en utilisant ce mode, on peut sa- 
voir quelle frame a ete perdue. S’il s’agit d’un message prioritaire ou importante, 
la station le renvoie. 
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Chapitre 3 


Modele mathematique 


Ce chapitre aborde le modele mathematique du mouvement d’un corps rigide dans 
l’espace a trois dimensions. On pent considerer le drone coninie un corps rigide 
snr leqnel agissent des forces et des couples exterieurs. Ces forces et couples d’une 
part contrent la force de gravite, d’autre part stabilisent l’attitude, ce qui permet 
au drone de "voler" dans l’espace a trois dimensions. En changeant les valeurs des 
forces, on pent changer certains couples. Par consequent, l’orientation est modihee 
et done le corps rigide (en general) ou le drone (pour un cas precis) peut se depla¬ 
cer dans n’importe quelle direction de l’espace. Le modele mathematique presente 
est decrit par les equations dynamiques de Newton-Euler. Quelques categories de 
forces et de couples externes qui actionnent le corps rigide sont examinees : les 
forces et les couples a l’origine des actionneurs (les moteurs-helices dans nos cas), 
a l’origine de l’effet aerodynamique (la portance et la trainee de l’aile quand le 
drone vole a vitesse elevee), a l’origine de la gravite et a l’origine de la force 
d’Archimede (quand le corps rigide est plus leger que l’environnement). Ensuite, 
les forces et les couples qui agissent sur quelques modeles sont analyses. 11 s’agit 
d’un dirigeable, d’un helicoptere a quatre rotors, d’un avion convertible et d’un 
avion basculant tri rotors. Les resultats de ces analyses sont ensuite utilises pour 
construire des lois de commande dans le prochain chapitre. 

Ce chapitre presente aussi les outils mathematiques utilises pour representer l’at¬ 
titude, que Lon doit controler et stabiliser. En general, en temps reel, la boucle 
d’asservissement de 1’attitude est plus rapide que le controle de la trajectoire. 
L’attitude est ainsi stabilisee a une consigne donnee par la direction desiree du 
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mouvement. On peut citer ici quelques methodes pour representer l’attitude : la 
matrice de rotation, les angles d’Euler, le vecteur des quaternions ... La matrice 
de rotation, qui a 9 elements, est plus compliquee a utiliser que les angles d’Euler 
(angles de roulis, de tangage, de lacet) qui permettent un affichage de 1’attitude 
plus simple. Cependant, ce formalisme presente le probleme de singularite quand 
1’angle de tangage s’approche de ±90°. Le vecteur de quaternion ayant 4 elements 
est capable d’eviter le probleme de singularity et de representer l’attitude dans 
n’importe quelle orientation. L’inconvenient du vecteur de quaternion est que sa 
representation n’est pas facile a lire. En general, on doit le transformer vers les 
angles d’Euler pour afficher l’attitude. 


3.1 Representation de l’attitude 



Figure 3.1 — Deux systemes de coordonnees 


Pour analyser le mouvement d’un corps rigide dans l’espace, nous definissons deux 
systemes de coordonnees comme ceux montres sur la figure 3.1 : 

• E / = [ x-F yf z-f' ] est un ensemble de trois vecteurs orthogonaux repre- 
sentant le systeme de coordonnees inertiel, situe a un certain point dans 
l’espace. 

• E 6 = [ x b y 6 5* ] est un ensemble de trois vecteurs orthogonaux represen- 
tant le systeme de coordonnees fixe au corps, situe au centre de la masse du 
corps rigide. 
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3.1.1 Matrice de rotation 

La matrice de rotation R appartient au sous-espace des matrices orthogonales de 
dimension trois, appele gronpe orthogonal special, denote par SO(3) et defini par : 

SO (3) = {R\Re M 3x3 , R t R = / 3 ,det (R) = 1} (3.1) 

La matrice R identifie completement l’orientation du repere mobile E 6 par rapport 
au repere inertiel et clle perrnet de faire la transformation de coordonnees d’un 
vecteur d’un systeme de coordonnees a un autre. 


La vitesse de rotation 

Dans le cas ou le systeme de coordonnees E 6 tourne par rapport au repere E^" avec 
une vitesse angulaire u = [ co x co y co z \ , la matrice R devient variable dans le 
temps et elle est denotee R(t). La relation entre Lattitude du corps et la vitesse 
angulaire Co s’ecrit : 

R = — lo x R = —S(lo)R (3.2) 

oil x denote le produit vectoriel, S(.) represente un tenseur antisymetrique associe 
au vecteur C, denote par : 


/ Cl \ / 0 -Cs c 2 \ 

S{ 0 = sU = c 3 0 -Cl (3.3) 

V c 3 ) V -c 2 Ci o J 

L’equation 3.2 est connue comme l’equation cinematique, aussi nominee equation 
“Strapdown”. Elle est utilisee dans les systemes actuels de navigation. 

Remarque 3.1 : Notons qu’en utilisant I’operateur de matrice antisymetrique 
S( C) = —S(() T , le produit croise du vecteur en 3D peut etre exprime sous une 
forme matricielle comme S(a)b = a x b pour deux vecteurs a, b e M 3 . 


3.1.2 Les angles d’Euler 

Les angles d’Euler decrivent la sequence de rotation la plus utilisee dans l’aeronau- 
tique et la robotique. En utilisant les angles d’Euler, le passage d’un systeme de 
coordonnees E^' a un systeme de coordonnees E 6 s’effectue toujours en 3 etapes : 

1. Rotation de if autour de Zf (angle de lacet avec — n ^ if ^ tt) 

2. Rotation de 6 autour de yf (angle de tangage avec — | ^ 0 ^ |) 
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3. Rotation de (p autour de Xf (angle de roulis avec — tt ^ cp ^ tt) 

De la meme fagon, si nous allons du systeme de coordonnees E 6 vers le systeme de 
coordonnees E , la sequence (cp,9,ip) doit etre respectee. Les matrices de rotation 
qui represented chaque rotation sont donnees par : 


Ra ;,</> 


/1 

0 

0 ^ 


( cos 9 

0 

— sin 6 \ 

0 

cos <p 

sin0 

5 Ry,0 | 

0 

1 

° 

\0 

— sin (p 

sin0 j 


( sin 9 

0 

cos 6, ) 


Rz,ip 



sin ip 0 
cos ip 0 
0 1 


(3.4) 


En consequence, la matrice de rotation R qui decrit le passage du systeme de 
coordonnees E 6 au systeme de coordonnees E^" devient : 


R = 


R Rz,i^Ry } oRx,(f) (3.5) 

cos 9 cos ip — cos (p sin ip + sin <p sin 9 cos pj sin <p sin ip + cos <p sin 9 cos ip 
cos 9 sin ip cos <p cos ip + sin <p sin 9 sin ip — sin <p cos ip + cos (p sin 9 sin ip 
— sin 6 sin <p cos 9 cos <p cos 9 

(3,6) 


La representation de Pattitude a travers les angles d’Euler est la plus utilisee dans 
les communautes robotique et aeronautique grace a la simplicite pour decrire une 
rotation dans l’espace. Neanmoins, ellc a le desavantage de presenter des singulari¬ 
ties geometriques lorsque 9 — ±|. De plus, cette singularite geometrique introduit 
une singularite dans l’equation cinematique associee a cette representation [27] : 

( <p \ ( 1 tan 6 sin (p 

9 ] = I 0 cos (p 

i>J \o SI 


tan 9 cos 

— sin 

COS <j) 

cos 0 



(3.7) 


La singularite peut etre eliminee si l’orientation du corps dans Pespace peut etre 
limitee. Par consequent, cette representation n’est pas adequate lorsque le corps 
evolue dans tout Pespace. 


3.1.3 Quaternion 

Une autre representation de Pattitude est le quaternion unitaire aussi appele pa- 
rametres d’Euler. Le quaternion est une solution alternative au theoreme d’Euler, 
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lequel enonce qu’une rotation generate dans 1’espace peut etre realisee par une 
simple rotation 0 autonr d’nn axe de rotation e. Le quaternion unitaire est com¬ 
pose d’un vecteur unitaire e, nomine axe d’Euler et d’un angle de rotation 0 autour 
de cet axe. II est defini par : 


Q = 


cos f 
e sin | 



(3.8) 


oil 

^-= [q\ql + (fq = l,q= [q 0 f ] T , go G K, q G M 3 | (3.9) 

q — \ qi q 2 qz ] T et go sont la partie vectorielle et la partie scalaire du quater¬ 
nion. 


Le quaternion identite et le quaternion conjugue sont definis par : 

Qid=[ 1 0 T ] T ,q=[q 0 -f] T (3.10) 

Par consequent, connne le quaternion est unitaire, on a : g -1 = q. 

La multiplication de deux quaternions quelconques qi = [ qi 0 q[ ] T et q 2 = 
[ q- 2 0 qJi ] T es t definie par : 


q\ ® q 2 


Qio ~Qi 

qi hqi 0 + S (gl) 



(3.11) 


La matrice de rotation R peut etre exprimee, en tenues de quaternion par : 


done : 


R = R ( q ) = [ql - g T g) J 3 + 2 (gg T - q 0 S (q)) 


R = 


1 - 2 {ql + ql) 

2 {Q1Q2 + %%) 
2 (gift - % ( h) 


2 {Q1Q2 - %%) 

l- 2 (ql + 4) 

2 {QoQi + q.2%) 


2 {%q.2 + <h%) 

2 (Q2Q3 - Q0Q1) 

1 - 2 (ql + ql) 


(3.12) 


(3.13) 


La conversion entre le vecteur de quaternion et les angles d’Euler est donne par : 

cos {(j)/ 2) cos {6/2) cos (-0/2) + sin {</>/2) sin {6/2) sin (0/2) 

_ sin (0/2) cos {9/2) cos (0/2) — cos (0/2) sin {9/2) sin (0/2) 

’ / cos (0/2) sin {9/2) cos (0/2) + sin (0/2) cos {9/2) sin (0/2) 
cos (0/2) cos {9/2) sin (0/2) — sin (0/2) sin {9/2) cos (0/2) 


(3.14) 
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et 


' 0 " 


t 2(<?09l+929 3 ) 

1-2 (ql+ql) 

6 

= 

arcsin (2 (q 0 q 2 - q 3 q J) 

. 0 . 


arctan 



1 _ 2 (< ? 2 + 9 2 ) 


(3.15) 


Remarque 3.2 : De Vequation (3.12), il est clair que nous obtenons la meme 
matrice R(q) pour le quaternion q et pour le quaternion —q, i.e. R(q ) = R(—q). 
En effet, les quaternions q et —q representent la meme attitude physique. Cepen- 
dant pour q, la rotation est effectuee avec un angle (3, tandis que pour — q elle est 
effectuee avec un angle 2n — f3, autour du vecteur e. 

La representation de Pattitude en utilisant le quaternion est parfaitement adaptee 
car avec quatre parametres et une contrainte, il est possible de representer l’at- 
titude d’un corps rigide de fagon globale et sans singularity. En fait, Pequation 
cinematique (3.7) maintenant en termes de quaternion est donnee par : 

(io\ = i ( -<f 

V Q ) 2 \ hqo + S (q) 


u 


(3.16) 

(3.17) 


3.1.4 Erreur d’attitude 

Considerons deux attitudes d’un corps rigide, parametrees par les matrices de 
rotation Ri et R 2 respectivement. Alors, Pattitude relative entre ces deux attitudes 
est calculee par : 

R r = R( 1 R 2 (3.18) 

En fait, R r represente un operateur d’attitude qui fait tourner R - 2 en R.\. A partir 
de cette observation, en general, l’attitude relative est utilisee dans le cadre de 
l’estimation et de la commande d’attitude comme une mesure d’erreur d’attitude. 
Dans ce contexte considerons Rd = R\ l’attitude desiree du corps rigide et R = R - 2 
l’attitude reelle du corps. En consequence, Perreur d’attitude est calculee par : 

(3.19) 


R e = R d l R 

Si Perreur d’attitude vaut zero, alors, R e = I 3 . 
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Lorsque le quaternion unitaire est utilise pour representer 1’attitude du corps, 
l’attitude relative entre q\ et q 2 s’exprime par : 


Qr = Qi 1 


Q2 = 


Qio 

~Q\ hqi 0 - S (gl) 


<? 2 0 

Q2 


= qi ® Q2 


(3.20) 


Dans le cadre de la commande d’attitude, si nous notons qd = q\ 1’attitude desiree 
du corps et q = q 2 l’attitude reelle du corps, l’erreur d’attitude est representee 
par : 

q e = Qd 1 ® q = q d ® q (3.21) 

Dans le cas on l’erreur d’attitude vaut zero, le quaternion d’erreur a deux valeurs 
possibles : 

q e =[±l 0 0 Of (3.22) 

Ceci est du a la non bijection du quaternion avec le groupe SO(3). 


3.2 Forces et Couples Aerodynamiques 

3.2.1 Forces et Couples Aerodynamiques d’une aile 


Bord Corde 

d'attaque moyenne 



Corde 


A 


Direction de Flux 




■►Trainee 


Figure 3.2 - Notation du profil d’aile 


Dans cette partie, la force et le couple generes par une reaction dynamique du 
mouvement de Pair seront etudiees. Principalement, cette reaction est le facteur 
contribuant des forces aerodynamiques de l’avion et des helices de l’actionneur. 
La section transversale d’une aile montree dans la figure 3.2 est appelee profil de 
l’aile. Le profil est utilise pour analyser la reponse aerodynamique de l’aile. 
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L’air aborde le profil par le bord d’attaque et le quitte par le bord de fuite. La par- 
tie superieure du profil est appelee extrados et l’interieur intrados. Le segment qui 
joint le bord d’attaque et le bord de fuite est appele le corde du profil. La distance 
entre l’intrados et l’extrados est l’epaisseur. L’epaisseur relative est la valeur de 
bepaisseur maximale divisee par la longueur de la corde, e re i = emax /^ . II existe 
aussi la corde moyenne qui est la ligne equidistante entre l’intrados et l’extrados. 
Elle donne la cambrure desiree, c’est pourquoi elle est connue comme la ligne de 
cambrure. La distance maximale entre la corde et la corde moyenne est appelee 
fleche du profil. Le rapport entre la fleche et la longueur de la corde est appelee la 
courbure (ou cambrure relative). 

L’ecoulement de l’air autour de l’aile est represente par lignes de flux. L’angle 
d’attaque a est bangle entre la corde et le vecteur de vitesse de l’air ou direction 
d’ecoulement. Cet angle est aussi connu comme angle d’incidence. Comme nous 
pouvons remarquer sur la figure 3.2, il y a deux forces et un couple qui sont le 
resultat du mouvement longitudinal. II s’agit de la portance, P, la trainee, T, et le 
couple de tangage m. Ils seront examines dans la suite. 


Force de Portance : 

La portance est definie comme la composante de force dans le plan de symetrie 
avec une direction perpendiculaire a la ligne de flux, elle est donnee par : 

P = I P V 2 SC P (3.23) 

ou p est la masse volumique de l’air, V est la vitesse de l’air, S est la surface de 
reference de l’aile et C p est le coefficient de portance (sans dimension) calcule par : 

C p = C Pa + C Ps (3.24) 

• C Pa depend de la configuration du profil de l’aile et sa variation par rapport 
a bangle d’attaque, totalement non lineaire et est identifiee par l’experimen- 
tation (voir la figure 3.3). 

• C Ps = C p 6 est le coefficient de portance dependant de la gouverne de l’aile, 
on peut considerer qu’il est proportionnellement lineaire a la deflection de la 
gouverne quand bangle d’attaque est inferieur a 20° (voir la figure 3.4). 

Dans la figure 3.3, on peut constater que la variation de C Pa par rapport a bangle 
d’attaque est presque lineaire jusqu’au point de decrochage (a < 10°), puis il 
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cjjdeg] 

Figure 3.3 — Coefficient de portance 

diminue quand a augmente. La pente de la courbe du coefficient de portance ne 
depend pas du nombre de Reynolds R e quand il est situe entre R e ~ 10 6 et 10 7 , 
mais elle est reduite quand le nombre de Reynolds est voisin de R e & 10 5 , comme 
c’est le cas d’un avion miniature. 



Figure 3.4 — Influence de la deflection de la gouverne au coefficient de portance 


Force de Trainee : 

En plus de la force de portance, il existe aussi une force appelee la trainee qui s’op- 
pose directement au mouvement de l’aile dans Lair. La trainee est une combinaison 
du frottement de Fair sur la surface, de Fintegrale de la pression sur la surface et 
de la trainee induite par des effets de tourbillons. La trainee est une fonction non 
lineaire qui depend fortement de Tangle d’attaque. En aviation, la trainee totale 
est la sonnne de la trainee parasite (non liee a la portance) et de la trainee in¬ 
duite par la portance. On peut quantifier la trainee totale par un coefficient global 
rapporte a la surface projetee de l’aile : 

T = ^ pV 2 SC T 


( 3 . 25 ) 
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ou C T est le coefficient de trainee (sans dimension) calcnle par 

C T = C Ta + C T ° (3.26) 

• C T ° est le coefficient constant de la trainee parasite (ou trainee de frottement) 

• C Ta est le coefficient de la trainee induite par la portance, varie non lineaire- 
ment par rapport a Tangle d’attaque et identifie par l’experimentation (voir 
la figure 3.5). 



«[deg] 


Figure 3.5 - Coefficient de trainee 


Couple de Tangage : 

II existe un couple agissant sur l’aile, autour d’un pivot dans le plan de symetrie, 
appele couple de tangage. Sa direction positive est montree dans la figure 3.2. Ce 
couple a tendance a retablir l’angle d’attaque. II est defini par l’equation suivante : 

m = l -pV 2 ScC m (3.27) 

ou c est la corde aerodynamique moyenne et C m est le coefficient du couple de 
tangage. Ce coefficient est une fonction de Tangle d’attaque et de la position de la 
gouverne de l’aile, il s’ecrit comme suit : 

C m = CZ + C r f l a + (3.28) 

• C™ est le coefficient du couple de tangage au centre aerodynamique ou foyer 
de l’aile 

• C™ est le coefficient du couple de tangage de l’aile par rapport a Tangle 
d’attaque a 
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• C™ represente les variations du couple de tangage par rapport a Tangle de 
la gouverne de l’aile 5 

Pourtant, ce couple n’existe que dans le cas ou le profil asymetrique de l’aile est 
utilise. Dans notre cas tous les profils utilises sont symetriques, done on lie prend 
pas en compte ce couple dans les equations dynamiques de l’avion. 


Modele longitudinal de la force aerodynamique : 



On examine ici le modele longitudinal de la force aerodynamique qui agit sur 
l’aile. Le modele dynamique longitudinal est montre dans la figure 3.6. Le flux 
d’air s’approche de l’aile avec une vitesse V air . Par l’effet aerodynamique il cree 
la portance P perpendiculaire a la direction du flux et la trainee T opposee a la 
direction du flux. A partir de cette portance P et de cette trainee T, on doit definir 
les forces de la portance P b et de la trainee T b dans le repere E^. Ils sont calculees 
par : 


P b = P cos a + T sin a 
T b = —P sin a + T cos a 


(3.29) 


ou : 

• a = 6 — A est Tangle d’attaque 

• 6 est Tangle de tangage 

• A = arcsin ^^j = arctan es t Tangle du vecteur de la vitesse du flux 
d’air par rapport a Thorizontale 

• P, T sont calcules a partir des formules aerodynamiques dans la table 3.1 


ou : 
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Table 3.1 — Forces aerodynamiques 


Portance 

Trainee 

P = - P V 2 SC P 

T = \pV 2 SC T 

2 

2 

C T = C Ta + C To 

C P = f P {a) 

C T " = f T (ot) 


• S est la superficie de l’aile 

• les functions f p (ct) et / T (a) sont les fonctions non lineaires des coefficients 
de portance et de trainee dependant de Tangle d’attaque et du profil des ailes 

• C T ° est le coefficient constant de la trainee parasite (ou trainee de frottement) 


3.2.2 Forces et Couples Aerodynamiques d’un actionneur 

L’objectif principal de cette partie est la determination des expressions des diffe- 
rentes forces et couples aerodynamiques produits par chaque actionneur. Ils sont 
exprimes, principalement, en fonction de la vitesse de rotation du rotor et de la 
geometrie de la pale. Nous considerons que les forces et couples de chaque ac¬ 
tionneur sont independants. Comme tous les actionneurs ont des caracteristiques 
similaires, on s’est interesse a un seul. 


La force de poussee totale generee par chaque actionneur peut etre calculee comme 
suit [23] : 


avec : 


NpacMS 2 R 3 M + -p ^^ <fo — (l + p 2 ) ^--A 


T = 

= b 2 s 2 + bi s + b { 


b-2 = Npac M R 3 M ( gV^o - g <Pt 


3 1 Fnci 


b\ = -Npac M R M - A D 

4 K m 

(\ 1 \ v? 

b 0 = Npac M R M ( |</?o - g <Ptw J ^2r 


M 


(3.30) 


(3.31) 
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Figure 3.7 — Composantes du couple resistant 


Le couple resistant ou couple de trainee genere par chaque actionneur est calcule 
par [23] : 


Q = Npac M s 2 R a m - 1 - (l + p 2 ) c d + A ( -<p 0 - ^-<p tw - 7 A 
8a \d 8 4 


= k 2 s 2 + hs + k 0 


avec : 


k 2 = Npac M R Ma¬ 
ttel 


k\ = -Npac M R M - 


lV f 


ind 


4 R 


k 0 = -Npac M R 3 M 


M 

4 R 2 m 


Po Q ^Ptv 


ou : 


(3.32) 


(3.33) 


• N est nombre de pales de l’helice 

• p est la masse volumique de l’air 

• a est la pente de la courbe de la force de sustentation 

• cm est la longueur de la corde de la pale 

• s est la vitesse de rotation 


















3.2 Forces et Couples Aerodynamiques 


64 


• Rm est le rayon de Thelice 

• p est la valeur adimensionnelle qui relie la vitesse de translation du drone 14 
a la vitesse de rotation de l’axe de l’actionneur est donnee par : 


0 


V h 

sRm 


(3.34) 


• (p est l’angle de pas ou de tangage qui varie lineairement par rapport a la 
position radiale r e selon la relation : 

T e 

P = Po ~ s—Vtw (3.35) 

Km 


ou p est Tangle de pas au voisinage de 1’axe de Tactionneur et p 0 est le 
coefficient de torsion. 

• A est le facteur qui lie la vitesse du vent induit 1 considere parallele a 
1’axe de Tactionneur a la vitesse de rotation de l’axe de Tactionneur : 


A 


Vi 


ind 


sRm 


(3.36) 


• Cd est le coefficient de trainee qui peut etre approxime en fonction de Tangle 
moyen d’attaque et du nombre de Mach a 75% de la longueur de la pale. Avec 
Tactionneur utilise, le nombre de Mach : M = — = (ou v s = 340.6 m/s 
est la vitesse du son au niveau de la mer) peut varier entre 0.0210 et 0.0273 
pour la valeur 75% mentionnee et pour les amplitudes de vitesse de rotation 
entre 200rad/s et 260 rad/s, ce qui donne un coefficient environ q = 0,01 
pour des angles d’attaque inferieure a 8°. 

Pour les actionneurs conime ceux utilises dans un systeme de drone, la vitesse 
moyenne de rotation peut atteindre une vitesse tres elevee (plus de 200 rad/sec). 
Dans les equations (3.30) et (3.32) le terme quadratique de la vitesse de rotation 
est predominant. Par consequent, celles-ci peuvent etre approximees par : 


T = b 2 s 2 
Q = k 2 s 2 


(3.37) 


Dans la suite, nous allons utiliser la notation b — b 2 et k — k 2 . Les valeurs de ces 
parametres varient selon la densite d’air, de rayon, de forme, d’angle de lancement 
de la lame et d’autres facteurs. 
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Parametre 

Description 

Valeur (USI) 

b 

parametre de la force de poussee 

29.1 x 10 -5 

k 

parametre de trainee 

1.14 x 10 -6 


Remarque 3.3 : 


II existe un couple de rotor qui apparait quand la vitesse de rotation est changee. 
Ce couple est appele le couple de reaction et est defini par : 


Q react Jrotor & 


(3.38) 


ou Jrotor est I’inertie des compos ants de rotation comme I’helice, l’axe du moteur 
et I’adaptateur d’helice. La direction de ce couple est la meme que celle du couple 
resistant. Pourtant, en general la vitesse du moteur ne change pas brusquement 
pendant le vol et en plus I’inertie des composantes de rotation J ro tor est tres petite. 
Done, l’amplitude de ce couple n’est pas importante pour etre prise en compte dans 
le modele mathematique du drone. 


Remarque 3.4 : 

Dans le cas ou Vavion convertible vole verticalement, il utilise les deflecteurs de la 
surface de controle des ailes qui sont immerges dans le flux d ’air cree par l ’helice 
("slipstream, flow" en anglais) pour stabiliser Vavion. Alors, nous devons estimer 
la vitesse de Pair creee par I’helice (inflow ou induced velocity en anglais) pour 
deduire la force necessaire pour la commande. 


L’helice est modelisee comme un disque d’actionneur sur lequel il y a une continuity 
dans la vitesse, mais il n’y a pas de continuity de pression. La vitesse du flux d’air 
lorsqu’elle traverse le disque est consideree uniforme a tous les points sur le disque. 
En utilisant la conservation de la masse et la conservation de l’energie a travers le 
disque (Theoreme de Bernoulli), on a : 


Po + Py = Pi + Py 


(3.39) 


La poussee de l’actionneur i est identifiee par : T x = (p 0 — pf) Am avec Am = vr R 2 M , 
Rm le rayon de l’helice, on obtient alors : 


—>■ 


T 

A M 



P 


2 


(3.40) 
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Figure 3.8 — Forces appliquees a un rotor 


2 i i 2 

vt = -b Vq 


in 


pA-M 


Done, la vitesse de Pair creee par la rotation d’helice est identifiee par 


^ In. 


2 Ti 9 / 26s? 

+ v 2 0 = 


pA-M 


pAm 


+ v o 


Nous avons deux cas de v 0 , la vitesse d’air devant de l’helice : 

• Le cas ou il n’y a pas de flux d’air devant l’helice, on a vq = 0 


-> Vin — 


2 T r 


26s? 


pAm V pA M 

Le cas ou il y a un flux d’air devant l’helice, on a v 0 — V 




2Ti + F 2 = \I^ + V 2 


pAm 


pAm 


La force Pf de la surface de controle immergee dans le flux d’air est 
partir de la formule aerodynamique : 


P- = - 2 ^'ls-cp 


(3.41) 


(3.42) 


(3.43) 


(3.44) 

calculee a 


( 3 . 45 ) 
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ou, Cp = C[ g 5i est le coefficient de portance par rapport a Fangle de deflexion et 
Cf est un parametre constant positif. Done, on a : 


->• 


Pi = ^stersi si v 0 = o 


p,‘ = {^ + p P)stcrfisiv„ = v 


->■ 


p t = « ”0 = 0 


Pt = (I S‘CfA si v„ = V 


l is 


(3.46) 

(3.47) 

(3.48) 


3.3 Modele dynamique d’un corps rigide 

II existe differentes approches pour decrire les equations dynamiques d’un corps 
rigide, a savoir, l’approche de Newton-Euler, l’approche de Lagrange ou autres... 
L’approche de Newton-Euler exploite les concepts de force et couples et e’est celle- 
ci que nous allons utiliser. Le mouvement d’un corps rigide dans l’espace 3D avec 
le systeme de coordonnees attache E & par rapport a la repere fixe Ef est decrit 
par les equations de Newton-Euler : 

rnv + mu) x fy; + mu) x v + mu) x (u) x re) = Pb (3.49) 

Iu) + mfc x v + mfc x (u) x v) + u) x (Iu)) = Q B (3.50) 

oil : 

T 

• fc — [ Xq Ug z g ] est le centre de gravite du corps rigide dans le repere 
E b 

• v — [ v x v y v z ] et u) — [ uj x uj y uj z ] sont respectivement les vitesses 
lineaires et angulaires 

• / G M 3x3 est la matrice inertielle par rapport a l’origine du repere fixe E b et 
m est la masse du corps rigide 

• Tb et Qb sont respectivement les vecteurs des forces et des couples exterieurs 
agissants sur le corps rigide 

Ces equations peuvent etre exprimees sous une forme plus compacte : 

M b V + C B {V)V = Tb (3.51) 


ou : 
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• V = [ vF u T ] 1 = [ v x v y v z u x u y uj z ] T est le vecteur de vitesse 

• Tb = [J^B Q T b ] T = [ F x Fy F z Ty r 2 ] T est un vecteur des 
forces et des couples externes 

• Mb et Cb(V) sont respectivement la matrice inertielle du corps rigide, et la 
matrice de Coriolis et centrifuge. 


3.3.1 La matrice inertielle de corps rigide 


La matrice inertielle de corps rigide Mb peut etre determinee des equations (3.49) 
et (3.50) : 


M r V = 


mv + mu x fc 


mv — mrc x u 



Iu + mrc x v 


Iu + mrc x v 



mv — mS(rc)u 


ml 3x3 -mS(f G ) 

V 

Iu + mS(fG)v 


[ mS(r G ) I 

J 

u 


(3.52) 


Done on a : 


ou : 


M b = 


ml 3x3 -mS(r G ) 
mS{fc ) I 


(3.53) 


m 

0 

0 

0 

mz G 

—my G 

0 

m 

0 

— mz G 

0 

mx G 

0 

0 

m 

my G 

—mx G 

0 

0 

—mz G 

my G 

IXX 

— Ixy 

Ixz 

rnze 

0 

—mx G 

~Iyx 

lyy 

— lyz 

-my G 

mx G 

0 

1zx 

— Izy 

Izz 


(3.54) 


3.3.2 Les matrices des forces de Coriolis et centrifuges 

La partie Cb(V)V sont les forces et les couples a l’origine de la force de Coriolis et 
de la force centrifuge : 


C B (V)V 


mu x v + mu x (u x fa) 
uj x (Jen) + mrc x (tJ x v) 


(3.55) 


ou : 
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• La partie de la force de Coriolis : 


mu x v 
rnfc x (u x v) 

La partie de la force centrifuge : 

mu x (u x f G ) 

Ld X (. Iu) 

mu x v = —mv xu = —mS(v)u 


(3.56) 


(3.57) 


0 

-v z 

1 


0 

mv z 

—mv y 

Vz 

0 

~v x 

u = 

—mv z 

0 

mv x 

. ~ v v 

V x 

0 


mv.y 

—rnv x 

0 


—> mu x v = — m 

• mu x(ux r G ) = mS(u)(u x fa) = —mS(u)S(rG)u 
—» mu x (u x fa) = — m 


u 


(3.58) 


1 

O 

1 

£ 

e 

_i 


0 ~z G y G 


uj z 0 — (jj x 


z G 0 -x G 

u 

LOy (jJ X 0 


_ -y G x G 0 


[y G Uy + z G u z ) 

-mx G u y —mx G u z 


—my G u x m(z G u z + x G u x ) —my G u z 
—mz G u x —mz G u y m(x G u x + y G u y ) 

(3.59)' 


Done : 


—> mu x v + mu x (u x f G ) = 

m(y G u y + z G u z ) -m(x G u y - v z ) -m(x G u z + v y ) 
-m(y G u x + v z ) m(z G u z + x G u x ) -m(y G u z - v x ) 


u 


(3.60) 


—m(z G u x - v y ) —m(z G u y + v x ) m{x G u x + yc^y) 
mf G x (u x v) = mS(r G )S(u)v — mS(v)v car mS(v)v = mv x v = 0 



0 

— ZG 

Vg 


0 

—u z 

Uy 

—> mr G x (u x v) — m 

Z G 

0 

-xg 


u z 

0 



Vg 

x G 

0 



^X 

0 


\ 0 

-v z 

Vy 

1 





-m 


'Vy V 


0 -v x 
0 


(3.61) 
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—m(y G ujy + z G oo z ) my G u x mz G u x 

mx G uj y ~m(z G uo z + x G co x ) mz G uj y 

mx G u z my G uj z -m(x G u x + y G u y ) 

; 0 mv z —mVy 1 


+ 


—mv z 0 


mv x 


v 


mv y —rnv x 0 


v 


(3.62) 


—* mf G x(Coxv) 


- m(y G u y + z g uj z ) 
m(x G u y - v z ) 
m(x G uj z + Uy) 


• Co x ( ICo ) = —S(ICo)Co 
car : 


m(y G uo x + n 2 ) m(z G io x - u y ) 

—m(z G io z + a: G o; x ) m(z G u y + t> x ) 
m(y G u z - v x ) —m(x G uj x + 2/g^) 

(3.63) 



Ixx Ixy Ixz 


^ X 


Ixz^z Ixy^y Ixx^x 

ICo = 

i 

a 

1 


UJy 

= 

lyz^z l-xy^x H - ^-yy^y 


Izx Izy Izz \ 


. UZ . 


lyz^y Ixz^x Izz^z 


i_ «" _j i_ _i i_ _i 

(3.64) 

Done : 


—^ ca x (Iuo) — 


lyz^y 4"~ 


0 

XZ^X J-ZZ^Z 


-I yz LO z 


Ixy^x 


+ I-\ 


lyz^Jy Ixz^x 


^x + -Tzz^z 


yyiOy 


Ixz^z C I. 


xyOOy 


r.OJx 


lyz^z 4“ ^xyX'x ^yy^y 


I x z LO 2 


Ixy^y 4~ ^x 


rZOr 


(3.65) 


ca 


Enfin, on a : 


mw x n + m3 x (ca x r G ) 
mf G x (Co x v) + Co x (ICo) 


C B (V)V = 

C B (V)V = 

Done, Cb(V) est donnee par : 

C B (V) = 


—mS(v)Co — mS(Co)S(r G )Co 
mS(r G )S(Co)v — mS(v)v — S(ICo)Co 

(3.66) 


0 

—mS(v) —mS(Co)S(r G ) 

V 

— mS(v) + mS(r G )S(Co) 

-S(ICo) 

Co 

(3 


0 — mS(v) — mS(Co)S(f G ) 

mS(v) + mS (r G ) S (Co) -S(ICo) 


(3.68) 
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-+C B {V) = 


0 

0 

0 

-m(y G ujy + zqu z ) 
m(x G Uy - v z ) 
m(x G u z + v y ) 


m(y G u y + z G u z ) 
—m(y G u> x + u z ) 
—m(z G u x - Vj,) 

0 


lyz^v “1“ h xz UJj 


L yz^y 

- Iy Z ! xJ Z 


z OJ -/ 


0 
0 
0 

m(y G u x +n 2 ) 
-m(z G UJ z + XG-CUa;) 
m(y G uj z - v x ) 

- m(x G u y - v z ) 
m(z G U z + ^G^z) 
—m(z G u y + i> x ) 

h-XZ^X P -7 22 CU 2 
0 


0 

0 

0 

m(z G uj x - v y ) 

m(z G ujy + i> x ) 

-m(x G u x + t/g^j/) 
—m(x G u> z + Vj,) 
-m(y G uj z - v x ) 
m[x G u x + t/g^d) 


■ Iy Z UJy 


IiizLVz P I.rii'-V 


yz* 

I nr ? 


1-xy^x ""I - lyv^v 1-xz^z “I - 


Lyy^y 


xy^y 


r.^cr 


xy^x 
I X yUJy + I' 

0 


IyyiVy 




(3.69) 


Remarque 3.5 : 


Une autre representation antisymetrique de Cb{V) est obtenue en utilisant Vequa¬ 
tion S{v)tJj = —S(oj)v. Par consequent, le resultat ci-dessus peut etre reformule de 
la maniere suivante : 


C B {V) 


mS{uj) —mS(oj)S(r G ) 
mS(r G )S(uj) —S(hjj) 


Notons que maintenant la matrice C'b(V) ne depend que de u. 


(3.70) 


Remarque 3.6 : 

On peut constater que si le corps rigide n’est pas symetrique et si le centre de gra¬ 
nite ne coincide pas avec le centre de coordonne attache Et,, le modele devient tres 
complique. Par consequent, les lois de commande ne garantissent pas de compenser 
exactement les forces et les couples de la force de Coriolis et centrifuge. En effet, 
le drone est generalement congu symetrique avec le centre de coordonne attache au 
corps rigide se situant au centre de granite. 


Dans ce cas, on a : 


r G = [ x G y G z G ] 7 = [ 0 0 0 ] T 


( 3 . 71 ) 
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et 


I 


0 Iyy 0 

0 0 I zz 


done, les matrices Mb et Cb deviennent : 


C B {V) 


C B {V) = 



m 

0 

0 

0 

0 

0 


0 

m 

0 

0 

0 

0 

M b = 

0 

0 

m 

0 

0 

0 

0 

0 

0 

Ivx 

0 

0 


0 

0 

0 

0 

Iyy 

0 


0 

0 

0 

0 

0 

Izz 




0 

g 

T 


g 

1 _ 

0 


— 

rriS (v) 

-S(L u) \ 


L 0 - 

■S(Icj) 

0 

—mco z 

mujy 

0 

0 

0 

mu z 

0 

—muj x 

0 

0 

0 

—nrnjj y 

moj x 

0 

0 

0 

0 

0 

0 

0 

0 

4z w z 

— IyyUJy 

0 

0 

0 

Izz 

u z 0 

1 xxUJj- 

0 

0 

0 

IyyUJy I XX^ X 

0 


(3.72) 


(3.73) 


(3.74) 


(3.75) 


3.4 Forces et couples agissant sur un corps rigide 

Les forces Tb et les couples Qb exterieurs agissant sur le corps rigide sont divises 
en quatre categories : 

• Les forces et les couples de l’environnement : il s’agit de la force de gravite 
et de la force d’Archimede. La force de gravite est a l’origine de la force 
attractive de la Terre et la force d’Archimede apparait si le corps rigide est 
plus leger que le fluide dans l’environnement il fonctionne (e’est l’air dans 
notre cas). 

• Les forces et les couples generes par les actionneurs : principalement par les 
helices tournantes. On considere que les moteurs sont les actionneurs qui 
donnent une poussee parallele a son axe de rotation et un couple resistant 
oppose a la direction de rotation de l’helice. 
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• Les forces et les couples generes par l’effet aerodynamique : ce sont la por- 
tance et la trainee de l’aile lorsqu’il existe un flux d’air qui approche 1’aile. 

• Les forces et les couples de perturbations de l’environnement et du vent... 
En fait, ces forces et couples sont difflciles a modeliser, done dans la these ils 
sont consideres comme des perturbations au systeme dynamique dn drone et 
ils ne sont pas pris en compte dans le modele dynamique. 

Alors, les forces Tb et les couples Qb sont donnes par : 

T b = W b + A* + F\ + F b (3.76) 

Q B = r G xW b + r A xA b + ? A + ? T (3.77) 

oil : 

• W b , fa x W b sont les vecteurs de la force et du couple issus de la force de 
gravite 

• A b , r A x A b sont les vecteurs de la force et du couple issus de la force d’Ar- 
chimede 

• F\,T A sont les vecteurs de la force et du couple issus de la force aerodyna¬ 
mique 

• F b , IV sont les vecteurs de la force et du couple issus de la force generee par 
les actionneurs 

On va examiner Faction de ces forces et couples exterieurs au corps rigide pour 
quelques exemples : le dirigeable, l’hclicoptere de quatre rotors, l’avion convertible, 
Lavion basculant a trois rotors. 


3.4.1 La force de gravite 

Dans cette section, la force de gravite est examinee. Cette force s’exerce sur tous 
les corps rigides. C’est la force qui tire le corps rigide vers la Terre et l’empeche de 
voler. Pour que le drone puisse voler, il faut avoir une force opposante a la force 
de gravite. On appelle le poids du corps rigide : W = mg. La force de gravite dans 
le repere E b est calculee par : 


La force produite par la gravite : 


r o i 


0 

i 

1 o 

= mgR T 

0 

-1 


W b = mR T g f = mR T 


(3.78) 
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• Si on utilise les angles d’Euler, on a R T : 


R t = 


cos 9 cos 0 cos 9 sin 0 — sin 9 

— cos 0 sin 0 + sin 0 sin 9 cos 0 cos 0 cos 0 + sin 0 sin 9 sin 0 sin 0 cos 9 

sin 0 sin 0 + cos 0 sin 9 cos 0 — sin 0 cos 0 + cos 0 sin $ sin 0 cos 0 cos $ 

(3.79) 


(3.78) ->■ W b = mg 
Si on utilise les quaternions, on a R T : 


sin 9 

— cos 9 sin 0 

— cos 9 cos 0 


(3.80) 


R 1 = 


1-2 (q% + ql) 2 (q x q 2 + q 0 q 3 ) 2 ( qi q 3 - g 0 g 2 ) 

2 (?i? 2 - %%) 1 - 2 (ql + q 2 3 ) 2 (g 0 g 1 + q 2 q 3 ) 

2 (f/o? 2 + 9 i? 3 ) 2 ((/ 2 (/ 3 - %q x ) 1-2 (ql + ql) 


(3.78) ->■ W b = mg 

Les couples produits par la gravite : 


2 (q Q q- 2 - m*) 
-2 (q 0 q 1 + q 2 q 3 ) 
2(ql + ql)-l 


(3.81) 


(3.82) 


fc x W b = tq x mR T gl = mS(r G )R T g f 
Si on utilise les angles d’Euler : 

(3.83) -> f G x W b = S{r G )W b = 


(3.83) 


o -Z G y G 


mg sin 9 

z G 0 -x G 


—mg cos 9 sin 0 

_ -y G x G 0 


—mg cos 9 cos 0 


—> r G x W h = mg 
Si on utilise les quaternions : 

(3.83) ->• f G x W b = S(f G )W b = mg 


z G cos 9 sin 0 — y G cos 9 cos 0 
z G sin 9 + x G cos 9 cos 0 
—yc sin 9 — x G cos 9 sin 0 


(3.84) 

(3.85) 



0 -z G 

Vg 


2 {%q 2 - f/dls) 

mg 

z G 0 

-X G 


-2 (q 0 q i + q 2 q 3 ) 


-yc x G 

0 


. 2 (qj + ql) — 1 


(3.86) 
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-> r G x W b 


mg 


2z G (q 0 qi + <Ms) + Vg (2 {q\ + g|) - 1) 
2z g (q 0 q 2 - q^s) - x G (2 (q\ + gf) - 1) 
-2J/g (9o9 2 ~ 9i 9 3 ) “ 2a; G (g^ + q. 2 q 3 ) 


(3.87) 


3.4.2 La poussee d’Archimede 

La poussee d’Archimede est la force particuliere qu’un corps rigide plonge subit 
en tout ou en partie dans un fluide (liquide ou gaz) sounds a un champ de gravite. 
Cette force provient de l’augmentation de la pression du fluide avec la profondeur 
(l’effet de la gravite sur le fluide) : la pression est pins importante snr la partie in- 
ferieure d’un objet innnerge que sur sa partie superienre, il en resulte une poussee 
globalement verticale orientee vers le haut. C’est a partir de cette poussee qu’on 
dehnit la flottabilite d’un corps rigide. 

Un corps rigide innnerge dans un fluide aura une force de flottabilite vers le haut 
agissant sur celui-ci equivalant au poids du fluide deplace, ce qui lui per met de 
hotter ou au moins avoir 1’air de devenir plus leger. Si la flottabilite depasse le 
poids, 1’objet hotte, si le poids est superieur a la flottabilite, l’objet plonge. Si la 
poussee d’Archimede est egale au poids, le corps rigide a une flottabilite neutre 
et peut rester a son niveau. La decouverte du principe de la flottabilite, qui est le 
resultat de la pression hydrostatique dans le liquide, est attribute a Archimede. 


Auto-stabilite : 

Un objet flottant est stable s’il tend a se rendre a une position d’equilibre apres un 
petit deplacement. En fait, les objets flottants ont generalement la stabilite verti¬ 
cale parce que si l’objet est pousse legerement vers bas, il creera une plus grande 
force de flottabilite. Cette force, inegale avec la force de gravite va pousser l’objet 
pour retourner a la hauteur d’equilibre. 

La stabilite de rotation est aussi un aspect important des corps rigides flottants. 
Etant donne un petit deplacement angulaire, l’objet peut revenir a sa position 
initiate (s’il est stable), s’eloigner de sa position d’origine (s’il est unstable), ou 
rester ou il est (neutre) (voir figure 3.9). La stabilite de rotation depend de la ligne 
relative entre les forces agissant sur 1’objet. La force d’Archimede agit sur le corps 
rigide au centre de flottabilite, qui est au centre du volume du fluide deplace. La 
force de gravite agit sur l’objet a son centre de gravite. Lin objet sera auto-stable 
si une rotation deplagant la ligne d’action de ces forces va generer un couple de 
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correction Yrm (righting moment). 



Figure 3.9 — L’effet de la stability : a) Configuration instable b) Configuration stable 


Si le centre de gravite et le centre de flottabilite du corps rigide font une ligne 
verticale, le corps rigide sera auto-stable. En fait, il existe un couple de correction 
quand le corps rigide est sounds a une petite rotation de Tangle de roulis ou 
tangage a cause des perturbations de l’environnement (voir figure 3.10). La capacite 
d’auto-stabilite de Tobjet flottant (par example : le dirigeable ou le vehicule sous 
marin) est tres avantageuse. Nous n’avons pas besoin d’utiliser des actionneurs 
pour stabiliser le corps rigide. En effet, le couple de correction Yrm est ecrit par : 


r RM = \d{A + W) (3.88) 

ou : d est la distance entre le centre de la force de gravite et la force d’Archimede. 




Figure 3.10 — Le couple de correction Yrm proportionnel a l’inclination de Tangle de 
roulis ou de tangage 

Yrm varie par rapport a Tangle d’inclination (Tangle de tangage ou roulis) : 

Y RM — 7 ^ {A + W) (sin cj) + sin 9) 


(3.89) 
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ou : l est la distance entre le centre de gravite et le centre de flottabilite. 


La force d’Archimede : 

La force d’Archimede, A b , est calculee de la formule : A = pgV, ou p est la 
masse volumique de l’air, V est le volume de fluide (l’air dans notre cas) deplace. 
On definit le centre de flottabilite, Fa = [ xa Da ] T , est la distance entre 
l’origine du repere E h et le centre de flottabilite. On analyse le calcul similaire avec 
la force de gravite : 

' 0 1 [ 0 

A 6 = R T A f = R t 0 =R t 0 (3.90) 

_ A J [ pgV _ 

• Si on utilise les angles d’Euler : 

— sin 6 

(3.90) —> ^ = pg\7 cos 9 sin <f> (3.91) 

cos 6 cos 4> 

• Si on utilise les quaternions : 

I" 1-2 (q* + qf) 2 (g x g 2 + q 0 q 3 ) 2 (q^ - q 0 q 2 ) 1 T 0 

(3.90) A 13 = 2 {q x q 2 - q 0 q 3 ) 1-2 (q\ + qj) 2 (q^ + q 2 q 3 ) 0 

2 (Q 0 Q 2 + Q 1 Q 3 ) 2 i ( h ( h ~ <M \) l-2 (g? + g|)J [pgW 

(3.92) 

r 2 (?i? 3 - %<h)' 

->•& = pgV 2 (q 0 q 1 + q 2 q 3 ) (3.93) 

.1-2(9? + ^) . 

Les couples produits par la force d’Archimede : 


r A x ^ = S(f A )Ah (3.94) 

• Si on utilise les angles d’Euler : 

[" 0 -z A V a 1 T — sin 6 

(3.94) ->■ r A xAh = pgV za 0 — xa cos 9 sin 0 (3.95) 

—yA xa 0 cos 6 cos 0 
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—y Fa x A 1 = pgW 


—za cos 9 sin <f> + ua cos 9 cos <f> 
—za sin 9 — x a cos 9 cos 0 
Ua sin 9 + xa cos 9 sin <f> 


• Si on utilise les quaternions : 


(3.96) 


(3.94) ->• r A x$ 



0 za 

Va 


" 2 (^3 - q 0 q 2 ) ' 

pgN 

za 0 

- x A 


2 (q 0 Qi + 9 2 9 3 ) 


_ -IJa x A 

0 


_ 1-2 (ql + ql) _ 


(3.97) 


—y r a x A 1 = —pgV 


2 z A (■ q 0 Qi + q 2 q s) + Vg ( 2 (q\ + ql) ~ 1) 
2 z A ( q 0 q 2 - <h%) ~ x a (2 (ql + ql) ~ 1 ) 

- 2 VA ( q 0 q 2 - 9 i 9 3 ) - 2 X A (q 0 q 1 + q 2 q 3 ) 


(3.98) 


Remarque 3.7 : 

Grace a la force d ’Archimede, le dirigeable est capable de s ’auto-stabiliser. En effet, 
on ajoute en plus un actionneur parallele a son axe principal pour le mouvement 
de translation et un actionneur pour commander Vangle de lacet. Les actionneurs 
utilises pour stabiliser l’attitude ne sont pas necessaires et sont done enleves. L’en- 
veloppe du dirigeable est consideree comme un actionneur, en changeant le volume 
de I’enveloppe, on pent commander Valtitude. Un autre solution pour stabiliser Val¬ 
titude sans changer le volume de I’enveloppe est d’ajouter un actionneur sur I’axe 
vertical. La capacite d’auto-stabiliser permet d’utiliser le dirigeable sans avoir le 
probleme de stabiliser l’attitude. Pourtant, le dirigeable est encombrant a cause 
de I’enveloppe et done il n’est pas pratique dans les essais. Notons que ce modele 
mathematique du dirigeable est aussi utilise pour les vehicules sous-marins car Us 
ont le meme principe. 


3.4.3 Helicoptere a quatre rotors 

Avec une structure mecanique simple, l’helicoptere a quatre rotors est de plus en 
plus utilise a la place de l’helicoptere classique, qui est complique a construire ou 
reparer a cause de la structure mecanique complexe du plateau cyclique. Comme 
montre dans la figure 3.11, les moteurs arriere et avant (Mi, M 3 ) tournent dans le 
sens contraire des aiguilles d’une montre tandis que les moteurs droit et gauche 
(M 2 , M 4 ) tournent dans le sens des aiguilles d’une montre. Chaque actionneur pro- 
duit une force /* parallele a son axe de rotation, ainsi qu’un couple resistant Qi 
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oppose au sens de rotation. 

La force totale on poussee totale exercee snr l’helicoptere (parallele a l’axe Zb) est 
la somme des qnatre forces generees par chaqne actionnenr (T — f\ + / 2 + /a + /r ) • 
La combinaison des forces f % et des couples resistants Qi donne origine aux mon- 
vements angnlaires antonr des axes principanx de l’helicoptere. Par consequent, 
les forces Fj, et les couples des moteurs sont donnes par : 

= [ 0 0 T] t 

r T = [ n r 2 r, ] T 

La vitesse des rotors pent atteindre des valeurs elevees (plus de 200 rad/sec ), 
done, le couple reactif produit par le rotor i peut etre approxime par Qj = ksj 
et la poussee par /,; = bsf , ou ,s t est la vitesse de rotation du rotor i, k > 0 et 
b > 0 sont deux parametres selon la densite d’air, de rayon, de forme, d’angle de 
lancement de l’helice et d’autres facteurs. 




Figure 3.11 — Helicoptere a quatre rotors 


• Fj. = [0 0 T ] T est la poussee totale exprimee dans le repere E 6 . La 
poussee totale produite par les quatre rotors est approximee par : 

4 

T = (3-100) 

i =1 

• Tt — [ Ti r 2 r 3 ] T est le vecteur des couples de commande. Les compo- 
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santes du couple de commande IV produit par les rotors sont donnes par : 

Ti = d(f 2 -U) = db(sl-sl) (3.101) 

r 2 = d(f 1 -f 3 ) = db(sl-sl) (3.102) 

T 3 = Qi — Q 2 + Q 3 — Qi = k (f — ^2 + S 3 — s^) (3.103) 

ou d represente la distance d’un rotor au centre de la masse (ou centre du 

repere EQ. 

A partir de la combinaison de (3.100) avec (3.101) - (3.103), les forces et les couples 
appliques a l’hclicoptere a quatre rotors sont ecrits sous forme de vecteurs : 



( 

0 

db 

0 

-db 

r ^ _ 


dh 

0 

-db 

0 

T ) 


k 

-k 

k 

—k 


V 

b 

b 

b 

b 



/ V \ 


‘l 


S 2 


„ 2 


S 3 


\ si / 


= N S r 


(3.104) 


3.4.4 Avion convertible 

L’avion convertible est un avion du type " Tail-sitter " . Dans la configuration, deux 
rotors contrarotatifs sont places cote a cote dans la partie superieure du vehicule 
et entre l’aile pour diminuer les risques d’accidents . Les couples de tangage, roulis 
et lacet sont obtenus par des ailerons et le gouvernail. L’avion a ete bien calibre 
pour que le centre de gravite coincide avec le centre aerodynamique de l’aile et ils 
se situent au quart du bord d’attaque de l’aile. 


Les forces Fi? et les couples generes par les actionneurs sont identiques dans 
les deux cas du vol horizontal et vertical. Comme montre dans les figures 3.12 et 
3.13, le moteur gauche tourne dans le sens des aiguilles d’une montre tandis que 
le moteur droit tourne dans le sens contraire des aiguilles d’une montre. Chaque 
moteur produit une force T t parallele a son axe de rotation, ainsi qu’un couple 
resistant Qi oppose au sens de rotation. La combinaison des forces T t et des couples 
resistants Qi est donnee par : 


F h T = [ F+T 2 0 Of 

f T = [Q2-Q1 0 {T 2 -T 1 )l m ] T 


oil l rn est la distance entre les moteurs et 1’axe principal. Par contre, les forces 
F\ et les couples aerodynamiques dependent du mode fonctionnaire de l’avion. 
Lorsque le drone fonctionne en vol vertical, les hypotheses suivantes sont suppo- 


sees : 
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• la magnitude de la force de trainee est negligeable par rapport a la portance 
et la poussee. 

• la force de portance de 1’aile est negligeable du a la combinaison du profil 
symetrique de 1’aile et a la direction de Pair de Phelice qui coincide avec la 
ligne de portance zero de l’aile. 

• on considere qu’il n’y a aucun vent lateral, done la force aerodynamique n’est 
fournie que par la deflection des gouvernes de controle. 

Ainsi, en mode helicoptere (voir figure 3.12), les forces et les couples aerodyna- 
miques sont calcules simplement : 

Fa = [ 0 - P ? PL + Plr ] T « [ 0 0 0 ] T 

?a = [ (PL - PL) I m (PL + PL) h PL. ] T 

ou : 

• Pf = lpvf n S-Cp avec i = al,ar,r ( left aileron , right aileron , et rudder ) : 
les forces aerodynamiques aux surfaces de controle generees par le flux d’air 
cree par l’helice 

• p est la densite de Pair 

• Cp = C[ s 5i est le coefficient de portance par rapport a l’angle de deflexion 
ou Cf s est un coefficient constant positif de portance 

• Si et Si sont des angles de deflexion et le superficie des surfaces de controle 

• v in = \Jptd es ^ v ^ esse d e l’ a i r P ar 1’helice 

• Am = 7T R 2 M) Rm est le rayon de Phelice 
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Figure 3.13 - Vol horizontal (mode d’avion) 


• l s est la distance entre la portance P/ et l’origine du repere Eb 

En vol horizontal, on suppose qu’il n’y a pas de vent. Alors, la vitesse de Pair 
est egale a la vitesse relative du drone, mais dans le sens contraire, c’est-a-dire 
V a i r = —V = — (yl + Vy + v(). Le vecteur de la vitesse de Pair cree avec Pavion 
un angle d’attaque a indique dans la figure 3.13. Dans ce mode, les forces et les 
couples aerodynamiques sont donnes par : 


A=[-r* -p; p' + ptr + Ph [ -T b o p b f 

pa = [ <ra - o (fa + o i. px ] T 


(3.107) 


Le portance P b et la trainee T h sont calcules a partir du modele longitudinal de la 
force aerodynamique (voir section 3.2.1). 


En ecrivant concretement, les forces et les couples agissant sur l’avion convertible 
sont : 

• Les forces P* ; , P 0 s r , P r s sont calculees a partir de la formule (3.44) : 

r2 


ps 
al 


Ti +P V 


PL = 


Pf = 


Am 

t 2 


A 


+ 


2 

pV 2 


M 

bs i I 
, Am 


S s al C^ l5 6 al = 


qs nr x _ 
ars° ar 


Li- + El ) S’fi&a (3.108) 


M 


bsi . pV 2 


hs l _|_ pv 2 \ 


S s r CL5 r si5 r > 0 


+ 


M 


StCLL si 5 r < 0 


5‘ r C ( U3.109) 

(3.110) 


2 J r s 
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Les forces F\ t — F\ + : 


Kt = 


Ti + T 2 - T b 
0 

P b 


(3.111) 


n T = 


Ti + T -2 + P sin a — T cos a 
0 

P cos a + T sin a 


(3.112) 


FXt = 


b (sf + s|) + \pV 2 S (C p sin a — C T cos a) 
0 

\pV 2 S (C p cos a + C T sin a) 


(3.113) 


Les couples T a,t = + IV : 


r a ,t = 


( P 0 - P ar) lm 


1 - 

6 

1 

<N 

s 

( Pal + ^ar) ^ 

+ 

0 

P r % 


(T 2 — Ti) l m 


(3.114) 


Ta,t — 


(P s al - P s ar ) i m + (q 2 - go 
(Pal + Par ) h 
Prh + (T 2 — Ti) l m 


“A f a,T — 


fcgi 1 pF 2 
A m 2 


(3.115) 


S^Sa, ~(^ + ^) Si r Ct,S.r) U + k(4 


K A 


■) + (S + "¥) 30O-) 

^ ^ 05^(0 + bl m {,s 2 - sj) si 5 , 

^ 0P r s CV^ r + bl m (. s I - s?) si <5, 


> 0 


On peut considerer que : C ,p _ = C a p r , = C p = Cf ; ^ ^ : 
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( s 2 S al - s 2 S ar ) + (S al - <5 or )) S s a Cfl m + k(s 2 - s 2 ) ' 

( s\5 al + s 2 2 S ar ) + + <J ar )) S-Cf Z, 

" (S + ¥) <J r + 6Z ro (s| - s 2 ) si S r > 0 

(H + *£) Z s S r -Cf <$ r + bl m (. s 2 - s 2 ) si5 r < 0 

(3.117) 


Remarque 3.8 : 


Notons qu’en vol vertical la vitesse de Vavion est faible (V ~ 0 ), les forces et les 
couples qui agissent sur l ’avion convertible dans le repere Eb sont donnes par : 



b (s 2 + ,1) 


0 

0 



Am 


(sl^al - S^bar) + k (s 


^ (sRi + s*6 a 


jk-S°Cfs 2 6 r + bl m (s 2 -s 2 ) 

bio 


L A 


£-SfCfs 2 S r + bl m (s 2 2 ~s 2 



si 5 r > 0 
si 5 r < 0 


(3.118) 


(3.119) 


En conclusion, les forces Eb et les couples Qb exterieurs agissant sur l’avion conver¬ 
tible sont : 


jr B = W b + F b A + (3.120) 

Qb — t'g x w b + + r T ( 3 . 121 ) 


Ponr eviter le probleme de la singularity, l’attitude est representee par les quater¬ 
nions dans le modele de l’avion convertible. Les forces Eb et les couples Qb sont 
alors donnes par : 

f 2mg (q 0 q 2 — q^) + b (s 2 + s|) + \pV 2 S (C p sin a — C T cos a) 1 


-A- Eb — 


-2 mg (g 0 <?i + q 2 q 3 ) 


mg (2 (qf + q 2 ) — 1) + \pV 2 S (C p cos a + C T sin «) | 

(3.122) 
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— > Qb — 


mg (2 z G (q^ + q 2 q 3 ) + y G (2 (qf + q\) - 1)) + 

{s\5 al - s 2 2 5 ar ) + ^ (S al - <5 or )] S s a C[l m + k(4- sj) 

mg (2 z G ( q 0 q 2 ~ Q 1 Q 3 ) ~ X G (2 (qf + qf) - 1)) 

+ ( S 1 &al + S^ar) + (5 a l + Sar ) j S*C[l s 

mg (-2 y G (q 0 q 2 - q 3 q 3 ) - 2x c (q^ + q 2 q 3 )) 


( 1 pv 2 \ 
\A m " l_ 2 


b4 pv^_\ 

V a m 2 


J l s S s r C[S r 


J l s S s r C[S r 


+ bl r 
+ bl r 


( s 2 ~~ s i) si 5 r >0 


) si 5 r < 0 


{4 ~ s\ 


(3.' 


23) 


Les parametres du modele de l’avion convertible qui seront utilises dans les simula¬ 
tions en boucle fermee dans le chapitre 4 sont detailles dans le tableau 3.2. Notons 
que les coefficients de portance et de trainee sont donnes dans les figures 3.3, 3.4 
et 3.5 en fonction de Tangle d’attaque et que les coefficients concernant les couples 
des moteurs des equations 3.37 ont ete obtenus de la reference [20]. 


3.4.5 Avion basculant a trois rotors 



Figure 3.14 - Vol vertical (mode d’helicoptere) 


Le modele dynamique des forces et des couples qui agissent sur Tavion basculant 
tri-rotor est montre sur la figure 3.14 (cas du vol vertical) et la figure 3.15 (cas du 
vol horizontal). Premierement, on analyse les forces F b T et les couples generes 
par les actionneurs. La combinaison des forces T % et des couples resistants Qi est 
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Table 3.2 — Parametres du modele de l’avion convertible 


Parametres 

Description 

Valeur (USI) 

m 

masse de V avion 

0.85kg 

F b T 

Forces generees par les actionneurs 

equation (3.105) 

r T 

Couples generes par les actionneurs 

equation (3.105) 

T % 

Poussee du moteur i 

equation (3.37) 

Qi 

Couple du moteur i 

equation (3.37) 

b 

Parametre de la force de poussee 

29.1 x HT 5 

k 

Parametre de trainee 

1.14 x 10" 6 


Vitesse de rotation du moteur i 

rad/s 

F\ 

Forces aerodynamiques 

equations (3.106) et (3.107) 

f .4 

Couples aerodynamiques 

equations (3.106) et (3.107) 

Pi = >ls;c? 

Forces aerodynamiques dues 

aux surfaces de controle 

equation (3.45) 

Vln ~ \J pA M 

Vitesse de Fair creee par Fhelice 

equation (3.43) 

Am = ttR m 

Surface de Fhelice 

0.00053m 2 

Rm 

Rayon de Fhelice 

0.013m 

Cf‘ = c?6, 

Coefficient de portance par rapport 
a l’angle de deflexion 


Cfs 

Coefficient constant positif de portance 

0.005 

Si 

Angles de deflexion 


s; 

Superhcie des surfaces de controle 

S s al = 0.045m 2 

S s ar = 0.045m 2 

S s r = 0.032m 2 

Is 

Distance entre la portance Pf 
et Forigine du repere E h 

0.03m 


Distance entre les moteurs et 

Forigine du repere E h 

0.015m 
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Parametres 

Description 

Valeur (USI) 

P 

Densite de l’air 

1.2 kg/m 3 

V air = —V 

Vitesse relative du drone 

m/s 

p = \pV 2 SC p 

Force de portance 

equation (3.23) 

s 

Superficie de l’aile 

0.4 m 2 

c p 

Coefficient de portance 

Figure 3.3 

T = \pV 2 SC T 

Force de trainee 

equation (3.25) 

C T 

Coefficient de trainee 

Figure 3.5 

a = 9 — X 

Angle d’attaque 


9 

Angle de tangage 


A = arcsin j 

Angle du vecteur de la vitesse 

du flux d’air par rapport a 1’horizontale 



donnee par : 


( T?b 


= 


T r = 


[ (Ti + T 2 ) cos 5 0 (Tj + T 2 ) sin 5 + T 3 ] 
(T 2 - Ti) l m sin 6 + (Q 2 - Qi) cos <5 
Tsh-in + tyhsmS 
(Ti - T 2 ) lm cos 5 + (Q 2 - <2i) sin <5 + Q 3 


(3.124) 


ou ImMM sont les distances des moteurs a l’axe x et y dn repere E b indiquees 
dans la figure 3.14, <5 est Tangle de bascnlement de denx moteurs avant, ce qui 
determine le mode de fonction de l’avion : 

• En mode vertical ( <5 = 90°, figure 3.14), l’equation (3.124) devient : 


r F* = [ 0 0 Ti + T 2 + T 3 / 


r r = 


(T 2 - Ti) lm 
T 3 I 2 — (Ti + T 2 ) l\ 
( Q2 — Qi) + Q3 


(3.125) 
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• En mode horizontal (6 = 0°,T 3 = Q 3 = 0, figure 3.15), Fequation (3.124) 
devient : 

f 4= [ (Ti + T 2 ) 0 0 ] T 
l (3.126) 

(f T =[ (Q2-Q1) 0 (Ti —T 2 )l m f 

Les forces F\ et les couples Y a aerodynamiques dependent du mode de fonction- 
nement de l’avion. Dans toutes les phases de vol, nous supposons qu’il n’y a pas 
de vent dans l’environnement. Done, lorsque l’avion vole verticalement, il n’existe 
pas de forces et de couples aerodynamiques car la vitesse laterale du vehicule est 
faible. Pourtant dans les phases de transition et de vol horizontal, la vitesse de hair 
est egale a la vitesse relative du drone, mais dans la direction opposee : V air = —V. 
Le vecteur de la vitesse de Fair cree avec Favion un angle d’attaque a (figure 3.15). 
Dans ces modes, les forces et les couples aerodynamiques sont donnes par : 



Figure 3.15 - Vol horizontal (mode d’avion) 


n = 


r 4 = 


- { T w + T ar + T al + T e + T r) 

P r 

P^- P S r + ps i + pb + ps 

\ ps ] 1 psi (pb I ps\ J _ pb] ps] ] 
L 1 ar h rn ' 1 arm y 1 e ' 1 e ) 1 1 r h r \ 


(3.127) 


Notons que pour l’equilibre dans la dynamique longitudinale, les distances l e et l w 
doivent satisfaire cette condition [45] : 


p b ] — p b ] 

1 e 1 W Lw 


^r~ 

lit 


K 

Pb e 


P.. 


S « 
5. 


(3.128) 


Done, on a : 



[-(T‘+T‘) 0 P‘+P‘] T 

1 ps; 4 _ ps ] ps] ps] 1 T 

L 1 apm ~ 1 a rm 1 e <, e 1 r b r J 


(3.129) 
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En ecrivant concretement, les forces et les couples agissant sur 1’ avion basculant 
tri-rotor sont : 

• Les forces F\ t = F\ + Fj, : 


(Ti + r 2 ) cos 5 


-C^ + T e & ) ' 

0 

+ 

0 

(Ti + T 2 ) sin 5 + T 3 


pb + pb 

L tu 1 e 


(3.130) 


FXt = 


(Ti + T 2 ) cos 6~(T b + T b ) 

0 

(T + To) sin 5 + To + P b + P b 


(3.131) 


F b — 

r A,T — 


b (sl + 4) cos 5 + \pV 2 [(S V CZ + S e cr ) sin a - (S W C£ + S e C T e ) cos a] 

0 

L b (4 + 4) sin5 + bs l + \p y2 [{ s w c w + s e°e) cos a + (S W C^ + S e Cj) sin a] J 

(3.132) 


Les couples V a,t = Ta + 


(T-2 - Ti) l m sin 5 + (Q 2 ~ Q i) cos 6 


' ps ; 4 - P s ] 

1 ar l rn ~ ± a p m 

T 3 / 2 -(T 1 + T 2 )/ 1 sin <5 

+ 

P S l 

1 e 

_ (Ti - T 2 ) l rn cos 5 + (Q 2 - Qi) sin 5 + Q 3 _ 


psi 
i r < r 


(3.133) 


P s e 




si 

2 




e 


P . = \ (& +*¥) s ' c p sis -<° 

" (xi, + ¥) S ‘ C F » <*- > 0 

JDS _ ( bs 2 I pv 2 i QS r<P X 

F11 \Am ' 2 ) i ~ , al'^al s °a,l 

— P s — f i_ pV 2 A Qs fiP r 

ar ~ “r 2 y °ar^ar s u ar 
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- On peut considerer que : 

r~<P _ r~<P /~iP _ /~iP 

0 ar s ~ ( ~ / r l5 ~ U e s ~ L 'S 

Os _ Cs Qs 

J al ~ J ar ~ °a 

8 ar V al 8 a 

- En mode de vol vertical (8 = 90°) : 



(3.133) —* Ta,t v — T t — 


b ( s 2 - s l) lm 
bs\l 2 -b(sj + sl) l\ 
k (sl — sf) + ks\ 


(3.134) 


(3.135) 


- En mode de transition, le flux d’air genere par les actionneurs n’agit pas 
les surfaces de controle (les ailerons, l’elevateur et le gouvernail) : 


(3.133) —> T a,t t 


b (sl — sf) l m sin 5 + k (sl — sf) cos 8 + pV 2 l m S^C[8 a 

bs\l 2 -b(sl + sl) h sin 5 + ^ l e S s e C[8 e 

b (sf — sl) l m cos 8 + k (sl — sf) sin 8 + ks\ + ^-l r S^.C[8 r 

(3.136) 


- En mode de vol horizontal (8 = 0°,S3 = 0 rad/s), le flux d’air genere par 
les actionneurs agit sur les surfaces de controle : 


(3.133) -> ? A)Th 


s\) + 



lmS S a C[8 



1 

2 



USiCfS 


e 


b (sl - sl) l m +(^ + ^) l r S s r C[8 r Si8 r < 0 

b (sl - sl) l m + (£j + P ^~) l r S s r Cf8 r si8 r > 0 

(3.137)' 


Les parametres du modele de l’avion basculant tri-rotor qui seront utilises dans les 
simulations en boucle fermee dans le chapitre 4 sont detailles dans le tableau 3.5. 
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Table 3.5 — Parametres du modele de l’avion basculant 


Parametres 

Description 

Valeur (USI) 

m 

masse de l’avion 

1.25kg 


Distance entre les 2 moteurs d’avant 

a l’axe x du rep ere E b 

0.16m 

h 

Distance entre les 2 moteurs d’avant 

a l’origine du repere E b 

0.23m 

h 

Distance entre le moteur arriere 

a l’origine du repere E h 

0.48 m 


Distance entre le centre de portance de l’aile 
a l’origine du repere E b 

0.08m 

le 

Distance entre le centre de portance de l’elevateur 
a l’origine du repere E h 

0 .68m 

lr 

Distance entre le centre de portance du gouvernail 
a l’origine du repere E b 

0.75m 

s w 

Superhcie de l’aile 

0.24 m 2 

Se 

Superhcie de l’elevateur 

0 .02m 2 


3.5 Conclusion 

Ce chapitre a presente en detail le modele mathematique du mouvement d’un corps 
rigide dans l’espace tridimensionnel. Ce modele est applique pour la modelisation 
des drones : le dirigeable, l’helicoptere a quatre rotors, l’avion convertible et l’avion 
basculant tri-rotor. Les forces et les couples agissant sur les modeles du prototype 
sont bien analyses pour chaque exemple. Ces analyses sont ensuite utilisees pour 
construire des lois de commande dans le chapitre prochain. 
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Chapitre 4 


Loi de commande 


L’objectif de ce chapitre est de presenter les lois de commande developpees. A 
partir de I’analyse dn modele mathematique decrit dans le chapitre precedent, on 
propose ici des lois de commande pour controler un corps rigide dans hespace a 
trois dimensions. Ces lois de commande sont ensuite appliquees pour les configu¬ 
rations de drones suivantes : l’helicoptere a quatre rotors, l’avion convertible et 
havion basculant tri-rotor. Les lois de commande developpees sont simples, ne ne- 
cessitent pas un cout eleve de calcul pour une boucle de commande et elles sont 
totalement appropriees pour les applications embarquees. Notons qu’on utilise un 
microcontroleur ( dsPICSSFJ256GP710 ) limite en ressources pour effectuer toutes 
les taches du systeme embarque. En dehors des lois de commande, le microcon¬ 
troleur doit effectuer plusieurs autres functions, par exemple la communication 
avec la Centrale Inertielle, avec la station sol, la commande de vitesse des moteurs 
etc. Done, si les lois sont compliquees, elles consomment des ressources et peuvent 
egalement produire des erreurs dans le programme du microcontroleur. Puisque 
les travaux pour debugger le programme du systeme embarque en temps reel sont 
difficiles (surtout impossible dans les phases de transition), des lois de commande 
simples ont ete developpees. Ces lois de commande simples doivent cependant ga- 
rantir la convergence asymptotique des erreurs vers zero, ce qui a ete prouve en 
utilisant une analyse basee sur la theorie de Lyapunov. Puis, dans ce chapitre, des 
methodes supplementaires comme la generation de la trajectoire pour le profil de 
vitesse, le changement des gains en utilisant les reseaux de neurones ameliorent les 
lois de commande proposees. 
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Les simulations developpees dans ce chapitre ont ete realisees en utilisant les pa- 
rametres des modeles presentes dans les tableaux 3.2 et 3.5 du chapitre 3. Les 
resultats de simulation nous montrent que les lois de commande sont performantes 
et peuvent etre appliquees dans la realite. Le programme de simulation est ecrit 
en C dans le bloc S-fuction Builder de Simulink, et peut ainsi etre importe direc- 
tement sur le programme embarque du drone. Le rnodele mathematique du corps 
rigide est realise par le logiciel VisualNastran fD , ce qui garantit la precision de 
rnodele complique du corps rigide. 


4.1 Contexte et etat de l’art 

Ces dernieres annees, des efforts significatifs de recherche ont ete diriges vers la 
commande, la modelisation et la conception de robots volants. Dans le cadre de 
cette these, nous nous interessons a des vehicules evoluant dans un espace a trois 
dimensions et pouvant etre modelises comme un corps rigide. Plus precisement, 
ce sont des vehicules sur lesquels agissent une force de poussee dans une direction 
privilegiee du vehicule et un vecteur couple permettant un controle complet de 
Lorientation. Cette classe de vehicules est parfois designee par le terme " thrust- 
propelled vehicles " dans la litterature anglo-saxonne. La commande de ces vehi¬ 
cules sous-actionnes (moins d’actionneurs que de degres de liberte) suscite depuis 
de nombreuses annees un grand interet pour des applications diverses et variees. 

L’un des avantages de ce type de vehicule est 1a, capacite de decoder et d’atter- 
rir verticalement. Ces vehicules aeriens sont done souvent appeles VTOL ( Vertical 
Take-Off and Landing). En fait, le fait que ces systemes evoluent dans un espace tri- 
dimensionnel (3D) pose de nouveaux problemes de recherche par rapport au cas 2D 
de la robotique mobile terrestre (vehicules a roues notamment) ou PVTOL ( Planar 
Vertical Take Off and Landing). En dehors des VTOLs (helicopteres, quadrotors, 
etc), e’est aussi le cas des dirigeables, fusees, missiles ou encore sous-marins. Evi- 
demment, chaque classe de ces systemes possede aussi des caracteristiques propres 
qu’il convient de prendre en compte. Celles-ci sont essentiellement liees a la nature 
des forces exterieures exercees sur le vehicule. Pourtant, cette similitude structu- 
relle peut etre exploitee dans un cadre general de commande. Cependant, ce travail 
n’a pas la pretention de couvrir 1a, commande de tous ces systemes, mais il essaie 
de donner une analyse commune et de proposer une approche de commande ex- 
ploitant la, structure d’actionnement commune a tous ces systemes. 

Le rnodele dynamique qui decrit le mouvement d’un vehicule sous-actionne dans 
l’espace tridimensionnel (six degres de liberte) peut etre decompose en trois sous- 
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systemes : la translation (£p), la rotation (Xa) et la dynamique des actionneurs 
(Em)- Le systeme en boucle ouverte du vehicule sous-actionne est represente a par- 
tir du diagramme bloc de la figure 4.1. Dans le cadre des VTOL, les actionneurs 
ont le but de fournir une force de poussee T le long d’une direction privilegiee 
du vehicule (c’est la force qui contre la force de gravite) pour generer le mouve- 
ment longitudinal et un vecteur T de couple permettant un controle complet de 
la dynamique de rotation. Pour les avions evoluant en mode horizontal, c’est la 
portance generee qui s’oppose a la gravite. Dans nos etudes, nous considerons que 
la dynamique de rotation est completement actionnee et la dynamique de trans¬ 
lation est sous-actionnee (une seule entree de commande T pour la dynamique de 
translation). Etant donne que la rotation conditionne le deplacement du vehicule 
dans l’espace, la commande (et/ou la stabilisation) d’attitude est privilegiee par 
rapport a la commande de position. En effet, une fois que l’attitude est stabilisee 
et controlee, le comportement du corps rigide est garanti et done le systeme est 
stable. 



T 


Figure 4.1 — Schema-bloc de la dynamique d’un vehicule sous-actionne (VTOL) 


Le probleme de la stabilisation d’attitude d’un corps rigide est un sujet qui a attire 
enormement d’interet au cours des 30 dernieres annees. Plusieurs approches ont 
ete appliquees telles que la commande linearisee par retour d’etat , [3], [4], la com¬ 
mande proportionnelle-derivee, [5], [6], la commande predictive dans [7] et [8], la 
commande integrate par backstepping dans [9] et la commande robustesse applique 
aux missiles tactiques, [10]. Parmi de ces approches mentionnees, la commande li¬ 
nearisee par retour d’etat et la commande proportionnelle-derivee sont largement 
appliquees. Cela est du au fait que ces lois de commande sont simples a mettre en 
oeuvre. La principale critique des deux dernieres approches est que pour une grande 
erreur d’attitude ou de vitesse angulaire, la loi de commande peut donner une en¬ 
tree de commande que les actionneurs sont incapables de produire. Ce fait devient 
evident en raison de leur limite physique. En pratique, les restrictions sur l’energie 
disponible imposent que les signaux d’entree soient bornes. En fait, lorsque l’erreur 
est assez grande, le signal de commande produit par ces lois est grand lui aussi. 
Ce dernier effet est dans plusieurs situations indesirables, voire dangereux car le 
signal de commande risque de provoquer une saturation voire rendonnnagement 
des actionneurs. Par ailleurs, il est souvent vrai qu’en raison de la limitation des 
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capteurs, la sortie du systeme est bornee. Pourtant, les approches de la commande 
d’attitude citees ci-dessus ne considerent pas le probleme qui prend en compte des 
contraintes sur les entrees. 

D’autres approches utilisant des lois de commande bornees sont alors proposees, 
par exemple dans [11] une loi de commande par saturation pour le systeme cine- 
matique d’un vaisseau spatial sous-actionne est etudiee. Cette loi de commande 
utilise une representation d’attitude non-standard, ce qui permet la decomposi¬ 
tion du mouvement general en deux rotations : une paire de coordonnees reelles 
et conjuguees. Dans [12], une loi de commande qui dirige un vehicule sous-marine 
d’une region initiate a une region finale de l’espace d’etat en satisfaisant les limites 
sur la commande et l’etat est proposee. Son approche exploite les idees de l’ana- 
lyse formelle des systemes continus et hybrides. Les auteurs dans [15] ont etuclie la 
robustesse de la stabilisation par mode glissant de la dynamique de l’attitude d’un 
vaisseau spatial en presence de la saturation de l’entree de commande base sur 
une approche VSC (Variable Structure Control ). L’application de la commande 
optimale au probleme de l’attitude d’un corps rigide a aussi interesse de nombreux 
chercheurs, (voir [16] et ses references). Malheureusement, les methodes citees ci- 
dessus sont couteuses en temps de calcul a cause du probleme de la resolution 
des equations de Hamilton-Jacobi. Cela explique la difficulte dans les applications 
pratiques. 

Dans le cadre de cette these, une loi de commande bornee pour la stabilisation 
d’attitude et de position d’un corps rigide dans l’espace tridimensionnel est propo¬ 
see. Avec l’approche proposee, les bornes de controle sur les six axes peuvent etre 
totalement librement choisis. La loi de commande bornee peut forcer separement 
la vitesse de chaque degre de liberte en boucle fermee vers une region determinee 
prioritaire de 1’espace d’etat et les maintenir dedans ensuite. Par ailleurs, en plus 
des contraintes de l’entree la loi de commande satisfait des contraintes de vitesse 
angulaire, c’est-a-dire des limites des gyrometres. La loi de commande est done 
robuste par rapport aux mesures bornees des vitesses angulaires, comme il existe 
avec des gyrometres classiques. La loi de commande est aussi robuste envers les 
erreurs de modelisation car elle ne depend pas de la connaissance des parametres 
d’inertie du systeme, ce qui est difficile a determiner dans la realite. La loi de com¬ 
mande proposee est inspiree de 1’approche des fonctions de saturation imbriquees 
proposee dans [17] pour le probleme de la stabilisation des systemes lineaires avec 
des connnandes bornees. 

A partir du modele preliminaire d’un degre de liberte, la loi de commande est 
developpee pour le modele d’un corps rigide a six degres de liberte dans l’espace 
a trois dimensions. Premierement, la loi de commande avec l’attitude representee 
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par les angles d’Euler est abordee. Les angles d’Euler sont simples a utiliser, pour- 
tant dans le systeme dont l’etat varie et peut approcher une region de singularite 
comme dans le cas de l’avion convertible, une nouvelle loi de commande d’attitude 
doit etre utilisee. C’est une loi de commande dans laquelle l’attitude est para- 
metree en utilisant le vecteur unitaire des quaternions, ce qui permet des grandes 
manoeuvres non lineaires des angles d’attitude. Cela fait aussi une difference avec le 
cas utilisant trois parametres de representation des angles d’Euler, ou une singula¬ 
rite geometrique apparait. La principale propriete de la loi de commande proposee 
est son extreme simplicity par rapport aux autres approches existantes, appro¬ 
priate pour l’implementation embarquee, ou des ressources de calcul du systeme 
embarque sont limitees. 


4.2 Stabilisation d’un corps rigide dans l’espace 3D 

4.2.1 Modele preliminaire 

Modele simple d’un degre de liberte 








x d 


X 



l«—> 



e 


Figure 4.2 — Modele preliminaire d’un degre de liberte 

On examine ici le modele mathematique simple d’un degre de liberte, ou le corps 
rigide se deplace suivant l’axe x sous une force externe F x sans frottement et des 
perturbations comme montre dans la Figure 4.2. C’est le modele principal d’un 
corps rigide quand il se deplace suivant un degre de liberte. Dans les lois de com¬ 
mande proposees pour commander un corps rigide dans l’espace tridimensionnel, 
la force de chaque degre de liberte agissant sur le corps rigide sera compensee pour 
qu’il revienne au modele preliminaire. Notons que c’est un modele pour la dyna- 
mique de translation, la dynamique de rotation est examinee de la meme maniere 
en changeant simplement la force par le couple et la masse par l’inertie de rotation. 

rnx = F x 


(4.1) 
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On dehnit des erreurs : 


e 

= x - x d 

(4.2) 

e 

= X - x d 

(4,3) 

e 

= X — Xd = x 

(4.4) 


On peut considerer que Xd = 0 car on n’a pas besoin de commander Pacceleration. 
En fait, si Xd ^ 0, on va le compenser en ajoutant une valeur x d a F x pour retour- 
ner au modele preliminaire. 

Theoreme 4.1 : La loi de commande non-saturee : 

F x = x d ~ k d e - k p e = -k d e - k p e (4.5) 

permet de converger x —> Xd quand t —* oo. kd et k p sont choisis par la methode 
de placement des poles pour que Pequation : 

e + k v e + k p e = 0 (4.6) 

ait deux racines dont la partie reelle est negative. 

La loi de commande ci-dessus n’est pas bornee, ce qui est delicate a mettre en 
oeuvre dans la realite. En effet, les actionneurs ont les limites physiques de leurs 
capacites, ce qui limite des valeurs des forces ou des couples generees. Alors, la loi 
de commande doit prendre en compte les contraintes physiques des actionneurs 
et par consequent, elle doit etre bornee. De plus, nous voulons ou devons parfois 
limiter les forces ou les couples generes pour des raisons de securite : si les forces 
ou les couples sont tres grands, le comportement du systeme n’est pas garanti et 
ces forces ou couples peuvent causer des accidents imprevus et depenser inutile- 
ment de l’energie. Dans le cadre de these, nous utilisons des lois de commandes 
bornees. Enhn, nous sommes interesses par la conception d’une loi de commande 
simple et continue, telle qu’elle puisse etre implantee en temps reel sur des systemes 
embarques ou la capacite de calcul est reduite. 


La loi de commande bornee 

Dans cette section, nous dehnissons la fonction de saturation, ensuite quelques lois 
de commande bornees de la litterature sont presentees. Elies sont applicables pour 
le controle du modele preliminaire, pourtant elles ne sont pas les plus simples a 
mettre en application. Nous proposons done une loi de commande bornee dont 
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la principale caracteristique est son extreme simplicity comparee aux autres ap- 
proches bornees existantes. La convergence asymptotique de la loi de commande 
est garantie par line analyse utilisant la theorie de Lyapunov. Dans le cadre de com¬ 
mande d’attitude, les lois de commande sont robustes par rapport a des erreurs 
dans la connaissance des parametres inertiels et la limite de mesure des donnees 
du gyrometre. 


Definition 4.1 (Fonction de saturation) : Soit M une constante positive. 
La fonction satu : M —>• R est une fonction continue et non decroissante telle que : 

_ / x M < M 

satM ~ \ Msign (x) si |x| > M [ } 

ou sign (x) est une fonction signe, definie par : 


sign(x) 


1 si x > 0 
—1 si x < 0 


(4.8) 


Quelques lois de commande bornees sont representees ci-dessus. Ces lois de com¬ 
mande sont originalement developpees pour commander un modele a n integrateurs 
et elles sont utilisees et adaptees ici pour controler le modele preliminaire ayant 2 
integrateurs. 

Theoreme 4.2 : La loi de commande bornee utilisant des fonctions de saturations 
emboitees [17] : 

F x = —satM 2 (oie + satMx («i«2e + a2e)) (4.9) 

permet de faire converger x vers Xd quand t —> oo, ou a\ et a 2 sont choisis par la 
methode de placement de poles suivant [18] et M\ < |M 2 . 

Dans [33], ces lois de commande sont modifiees en utilisant des fonctions de satu¬ 
rations separees. Chaque saturation joue sur une combinaison lineaire des etats du 
systeme. Elies peuvent aussi etre utilisees pour commander le modele preliminaire : 

Theoreme 4.3 : La loi de commande bornee utilisant des fonctions de saturations 
separees [33] : 

F x = —sat M2 (k 2 e) - sat Ml (h (k 2 e + e)) (4-10) 

permet de faire que x converge vers Xd quand t —> oo. 

Ces lois sont ensuite simplifiees dans [34] en utilisant des fonctions de satura¬ 
tions avec etats separes : 
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Theoreme 4.4 : La loi de commande bornee utilisant des fonctions de saturations 
avec etats separes [34] : 


F x = -sat Ml (k d e) - sat M2 (k p e) (4.11) 

permet de faire converger x vers x d quand t —> oo. 

Nous pouvons constater que nous devons utiliser deux fonctions de saturation 
dans les lois de commande bornees presentees ci-dessus. Ici. nous proposons une 
nonvelle loi de commande bornee qui utilise une seule fonction de saturation pour 
commander le modele preliminaire : 

Theoreme 4.5 : La loi de commande bornee : 

F x = —sat M (k d e + k p e) (4-12) 

permet de faire converger x vers x d quand t —>■ oo, ou M est choisi comme la force 
maximum identifiee par la limite physique de l’actionneur ou par notre besoin. 

Preuve : On a : 

rnx = F x = —satM ( k d e + k p e) (4.13) 

—> me — F x — —satM ( k d e + k p e) (4.14) 


Notons que : e —>• 0 x —> x d 
On definit la fonction de Lyapunov : 

V = ^me 2 + ^pe 2 (4.15) 

avec p > 0. On derive V : 

V = mee + pee = —esatM (. k d e + k p e ) + pee (4.16) 


On appelle : 

E = k d e + k p e (4.17) 

On suppose au debut de Pargument que la loi de commande F x est saturee, c’est- 
a-dire : 

\E\ > M > 0 —> me = F x = —Msign ( E ) (4.18) 

On definit la fonction de Lyapunov : 


14 = -mE 2 


(4.19) 
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-7Ve — mEE = rnE ( k d e + k p e) 

(4.20) 

-)• V E = E ( 

— k d Msign (. E) + mk p e ) 

(4.21) 

Notons que la condition pour que V E soit positif est : 


V E A 0 

e>m~ si E> 0 

TTlkp 

e< kd ™ si E < 0 

(4.22) 


— mk v 


Pour prouver qu’il existe un certain temps, tel que V E < 0 et done E diminue, on 
definit la fonction de Lyapunov suivante : 



—>■ Ve — mee = —eMsign (. E ) 

Si on suppose Ve > 0 de (4.22), on obtient : 

' V 6 = -eM < - < 0 si E > 0 

c — mkp 

Vi = eM < ^ k "]' < 0 si E < 0 

TTt kp 


(4.23) 

(4.24) 


(4.25) 


Alors, Ve < 0 avec la condition (4.22), e’est-a-dire e diminue et done il existe un 
temps t\ afin que e converge vers : 


k d M 

mk p 


< e < 


k d M 

mk p 


=>■ Ve A 0 


(4.26) 


Par consequent, il existe un temps tel que E converge vers : 

\E\ < M (4.27) 

et puis, la loi de commande n’est plus saturee : F x = — ( k d e + k p e). Done : 

—* V = — e (. k d e + k p e) + pee = —k d e 2 — k p ee + pee (4.28) 

En choisissant : p — k pi on a : 

V = -k v e 2 < 0 (4.29) 

Afin de completer la preuve, le principe d’invariance de LaSalle est invoque. Toutes 
les trajectoires convergent vers le plus grand ensemble invariant D dans 12 = 
j(e, e) : V = o| = {(e, e) : e — 0}. Dans cet ensemble invariant : 

me = 0 = F x = — (. k d e + k p e) = — k p e (4.30) 

Ainsi, 0 est reduit a l’origine : e = e = 0. Cela finit la preuve de la stability 

asymptotique du systeme en boucle fermee. En conclusion, la loi de commande 
bornee : F x = —satM ( k d e + k p e) est telle que : e —> 0 ou x —> x d . 
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Remarque 4.1 : 

Notons que si on peut representer le modele d’un corps rigide sous la forme de 
modeles preliminaires presente ci-dessus en compensant les forces et les couples 
exterieurs agissant sur le corps rigide, alors la loi de commande precedente peut 
etre utilisee pour commander separement chaque sous-systeme. L ’amplitude de la 
force et du couple qui s’exerce sur des axes principaux du corps rigide est bornee 
par M. La borne M peut etre choisie de fagon independante par chaque axe et elle 
depend des caracteristiques des actionneurs utilises. On peut constater que la loi de 
commande (4-12) est simple. Elle peut alors etre appliquee a la stabilisation d’un 
systeme d ’un corps rigide en temps reel avec des capacites de calcul contraignantes. 


Contrainte sur la vitesse angulaire : 

En pratique, la mesure de la vitesse angulaire dans chaque axe du corps est ac- 
complie par des gyrometres. Ces capteurs mesurent la vitesse angulaire dans une 
plage de mesure qui depend principalement de la technologie de fabrication. Par 
consequent, la capaeite a maintenir la rotation autour d’un axe sans depasser une 
certaine limite ou une borne dans la vitesse angulaire mesuree est essentielle. En 
fait, la loi de commande (4.12) permet de prendre en compte les contraintes sur la 
vitesse angulaire. 

Soit tc max la vitesse angulaire maximale 1 que le gyrometre peut mesurer lors d’une 
manoeuvre autour de l’axe de rotation. Considerons la contrainte suivante : 

|tc| ^ tC max 

Nous pouvons etablir la propriety suivante. 

Propriety 4.1 La loi de commande (402) stabilise globalement et asymptotique- 
ment Vattitude du corps rigide a Vorigine. De plus, si les parametres kd,k p sont 
choisis tels que la condition suivante est satisfaite : 

kd _ fWjnax 

k p M 

ou I est I’inertie de rotation du corps rigide sur l’axe examine, la contrainte 
M ^ Wmax reste satisfaite a partir d’un certain ti > t 0 si |cc(f 0 )| > w max et 
pour tout t si |ca(to)| < cc max . 

1. La Centrale Inertielle 3DM-GX3-25 utilisee dans le drone a une plage de mesure de vitesse 
angulaire de ±300°/s. 


(4.32) 
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Preuve : Supposons que les conditions initiates dn systeme sont telles que a l’ins- 
tant to, co appartient a l’ensemble $ = jcn : |cu| < j- A partir de la preuve 
du theoreme, nous pouvons observer que si co appartient a 4>, to va rester dans $ 
pour tout le temps t > to. Par consequent, si nous choisissons : jf = Jl ^ ax , nous 

pouvons conclure que |cu| < = c n max , Vt > t 0 . 

Par contre en realite, des perturbations externes s’exercent sur le corps rigide 
et dans certaines situations la vitesse angulaire du corps rigide depasse la plage de 
mesure du gyrometre : 

M > Wmax (4.33) 

Dans ce cas, les gyrometres seront satures et la loi de commande suivante doit 

faire converger la vitesse de rotation vers l’ensemble d 1 : 

T = -sat M {k d u max + k p e ) (4.34) 

oil e est l’erreur entre Tangle actuel et desire du l’axe examine. Rappelons le modele 
preliminaire de la dynarnique de rotation : 

Ico = T (4.35) 


On appelle : 


E k d co max "T kp€ 


On propose la fonction de Lyapunov : 




1 

2 


Ico 2 


(4.36) 


(4.37) 


V u — Icoco = —cosatM ( k d co max + k p e) = —cosatM (E) (4.38) 

On suppose pour commencer que la loi de commande est saturee : 

\E\ > M (4.39) 


La fonction de Lyapunov de E : 

V E = l -E 2 ^V E = EE = E (kpe) = E ( k p co ) 


(4.40) 


On a deux cas : 
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• Si Ve > 0 sign ( E) = sign (u) 

—>■ V u = —uisatM ( E ) = —c oMsign ( E ) = —c oMsign (u) < 0 (4.41) 

Done, il existe un temps t\ afin que la vitesse de rotation u converge vers 
l’ensemble <f> et reste dedans. En consequence le gyrometre n’est plus sature 
et le systeme evolue de la meme fagon que le systeme nominal 

• Si Ve < 0, il existe un temps t 2 afin que E converge vers zero et parce que 
|o;| > uj may , done e —>■ 0 


->V U = - UJSatM (feWmax) < 0 (4.42) 

car sign (uj) = sign (ta max ). Done, il existe un temps t 3 afin que la vitesse de 
rotation uj converge vers l’ensemble et reste dedans. 


4.2.2 Commande d’un corps rigide dans l’espace a trois di¬ 
mensions : 

Modele mathematique 

Dans le repere E b attache au corps rigide, le deplacement d’un corps rigide dans 
l’espace a trois dimensions est decrit par les equations dynamiques utilisant le 
formalisme de Newton-Euler (voir la section 3.3) : 

mv + muj x fa + mu x v + mu x (uj x re) = Fb (4.43) 

lu + mr G x v + mfc x (u x v) + u x (Iu) = Qb (4.44) 

oil [ ] T — [ F x F y F z T x Y y r z ] T sont les forces et les couples qui 

agissent sur le corps rigide. 

Quand le corps rigide se deplace dans 1’espace a trois dimensions, l’attitude definit 
l’orientation du mouvement. En general, la dynamique de rotation est complete- 
ment actionnee. C’est-a-dire que Y x , Y y et sont capables d’etre modifiees par 
les actionneurs. En revanche, la dynamique de translation peut etre actionnee ou 
sous-actionnee. Si elle est sous-actionnee, elle peut manquer d’actionneurs sur un 
axe (F y = 0) ou deux axes (F x — F y — 0). Dans ce chapitre, nous allons examiner 
tous les cas. Pourtant, nous devons d’abord garantir que l’attitude est totalement 
controlable et stabilise pour que le deplacement de la translation soit garantie. 
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Commande d’attitude : 

A partir des equations du mouvement de rotation : 

Iu + mr G x v + mr G x (u x v) + u x (Iu) = Q B (4.45) 

Iu — Q b — (mr G x ~v + mr G x (u x n) + cJ x (Iu)~j = Qc (4.46) 

En appliquant le theoreme (4.12), la loi de commande suivant permet la conver¬ 
gence de 1’ attitude : 0 = [ <f) 9 i/j ] T —» @ rf = [ cf) d 9 d ijj d ] T : 

Qc = sat Me (l< d Ee + (4.47) 

oil Eq, Eq sont respectivement le vecteur d’erreur de la vitesse de rotation et de 
l’attitude. 

Par ailleurs, on a les couples agissant sur le corps rigide (voir section 3.4) : 

Q b = r G x W b + r A x A 6 + f c (4.48) 

Oil r c = r A + r T est les couples de commande generes a partir des actionneurs. 

Done, en remplagant dans l’equation (4.46) : 

f G x W b + r A x A b + f c — ( mf G x v + mf G x (iu x v) + u x (Iu)) 

V ^ . \ ' (4.49) 

= -satM e [KfE q + KfE e J 

->• fc = —sat Me (^KfE e + KpEo'j + ^mr G x u + mf G x (cJ x u) + tJ x (Jen)j 

— ^'r*c x W b + C4 x A b j 

(4.50) 


Commande de position : 

A partir des equations du mouvement de translation : 

mv + mu x r G + mu x v + mu x (u x r*c) = -As (4-51) 

oil : = W b + + Fc est les forces agissant sur le corps rigide. Notons que 

F b G = F b A + Ft sont les forces de commande generees par les actionneurs. 
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On transforme le modele vers le repere inertiel Ej : 

R (mv + men x Eg + mu xi?+ mu x (en x Eg)] = R (w b + $ + Fq ) (4.52) 


—» mfi^ + xf^i men x n + men x (en x r G )^ = W f + A* + RF b , (4.53) 


—>■ rnv J 


E f = —R (mu x Eg + men x v + men x (en x r G ) j + A f + RF^ = F c 

(4.54) 

Si on choisit la loi de commande : 


F c =-sat M (K p d Ep + KPE P 


(4.55) 


Nous assurons la convergence de la position : p — [x y z ] T —> pd = [ Xd Vd Zd 
on Ep, Ep sont respectivement le vectenr d’errenr de la vitesse de translation et 
de la position. 

On doit identifier F b c , les forces agissent an corps rigide et generees par les ac- 
tionneurs : 


—> —R (men x Eg + men x v + men x (en x Eg) ) + W* + A* + RF ( 


c 


= -sat M „ ( K p d E P + K p E P 


(4.56) 


—> RF^ = —satM p (jF p Ep + KVEp^j + R (mu x Eg + men x v + men x (en x r G )j 

- (w> + -4') 

(457) 

On appelle qne la partie droite de Pequation (4.57) : 


(4.57) RF b , = 


h 

h 

fs 


(4.58) 


En fait, depend a la configuration des actionnenrs dn modele dn corps rigide : 
• Modele complet : le corps rigide a 3 actionneurs sur 3 axes x, y, z de E^ : 

Fc = in, n, n,] T 


(4.59) 
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• Modele sous-actionne suivant 1 degre de liberte : le corps rigide 
manque d’un actionneur sur Taxe y : 

F b c=[F* Cx 0 Fgf (4.60) 

• Modele sous-actionne suivant 2 degres de liberte : le corps rigide 
manque de deux actionneurs sur Taxe x, y : 

*c=[° 0 fg] T (4.61) 


On examine ensuite chaque cas : 


4.2.3 Modele complet : 


Le corps rigide a 3 actionneurs sur 3 axes x, y, z de E0 C’est le cas de l’octarotor 
(figure 1.8) par exemple. C’est le drone qui peut se deplacer dans le plan x,y sans 
necessairement incliner Tangle de tangage ou Tangle de ronlis et qui a le deconplage 
entre la dynamique de rotation et 1a. dynamique de translation. Dans ce cas, Fq 
est defini facilement : 


—>• Fq = R T satM p (^K%E P + RPEp'j + R (mu x vq + mu x v + mu x 

-R t (w f + 

^F b c = - R r sat Mp (^K^Ep + KPEpJ 

+ (mu x fa + mu x v + mu x (u x — (w b + A b j 


(u x f G ) 


(4.62) 

(4.63) 


4.2.4 Modele sous-actionne suivant 1 degre de liberte : 

Ce sont le cas du dirigeable qui est capable de s’auto-stabiliser, on ajoute gene- 
ralement un actionneur sur Taxe x pour commander la dynamique de translation. 

Nous avons en fait : F b , — [ F^ x 0 F^ z ] T . 

Car : 

' F b Cx 

0 

_ F c z 


RF b c = 


cos 9 cos if) — cos 0 sin 0 + sin 0 sin 9 cos 0 
cos 9 sin 0 cos 0 cos 0 + sin 0 sin 9 sin 0 
— sin 9 sin 0 cos 9 


sin 0 sin 0 + cos 0 sin 9 cos 0 
- sin 0 cos 0 + cos 0 sin 9 sin 0 
cos 0 cos 9 


(4.64) 



4.2 Stabilisation d’un corps rigide dans Pespace 3D 


108 


Puisqu’il y a un actionneur sur l’axe x, done il n’est pas necessaire d’incliner l’angle 
de tangage pour le deplacement lateral. Par ailleurs, la loi de commande d’attitude 
permet de converger Pangle de tangage vers zero, e’est-a-dire : 9 ~ 0. Done : 

cos ip — cos 0 sin ip sin <p sin ip i r f c* i 
RFq — sin0 cos <p cos ip — sin 0 cos ip 0 (4.65) 

0 sin (p cos 0 I L F C* \ 

Fp Jx cos ip + Fq z sin 0 sin ip f\ 

—>■ RF]p = Fp< x sin ip — Fp Jz sin 0 cos if) = f 2 (4.66) 

Fq z c°s 0 J U. 

Puisque la dynamique laterale sur Paxe y n’est pas actionnee, done au moins un 
autre degre de liberte de Pattitude doivent etre utilises pour realiser le deplacement 
sur cet axe. Les consignes desirees de ces degres de liberte sont instantanement 
fournies au sous-systeme de rotation a chaque boucle de commande. Ces consignes 
desirees sont identifiees comme suit. 

II y a 2 methodes : 

Utilisation de l’angle de lacet : Si on utilise Pangle de lacet, e’est-a-dire le 
cap du corps rigide se tourne vers la cible pendant le deplacement. Dans ce cas, 

Pangle de roulis est stabilise a zero degre : 0 « 0 et on a : 

T Fq x cos ip d 1 [ fi " 

RF b c = Fq x sin ip d = f 2 (4.67) 

F c z \ L /a _ 

Enfin, on deduit les forces des actionneurs et Pangle de lacet desire pendant le 
deplacement : 

( FL = VW+fl 

(4.67) —>■ < tan ip d = ^ —>■ ip d = arctan (4.68) 

{ 1 Fq z = / 3 

Utilisation de Pangle de roulis : Si on incline Pangle de roulis pour le depla¬ 
cement et la direction du cap est fixee a une direction predefinie : 

Fp Jx cos ip + Fg z sin (p d sin ip /i 

RFjp = Fp, x sin ip — Fp, z sin <p d cos ip = f 2 (4.69) 

Fp, z cos (p d J L/ s _ 
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{ tan (f) d = h sin 'P+h cos 'P _> 0 d = arctan 

F cx = A cos 0 + /2 sin 0 (4.70) 

F Cz = \J (/i sin0 + hcosijtf + f$ 

4.2.5 Modele sous-actionne suivant 2 degres de liberte : 

C’est le cas du VTOL (l’helicoptere a quatre rotors, helicopteres ...) ou il n’y a 
pas d’actionneurs places sur les axes x, y. Dans ce cas, on doit incliner Tangle 
de tangage et Tangle de roulis pour qu’il puisse se deplacer dans Tespace a trois 
dimensions. On a F^, — [ 0 0 F b , z ] T et : 

cos 9 cos 0 — cos 0 sin 0 + sin 0 sin 9 cos 0 sin 0 sin 0 + cos 0 sin 9 cos 0 0 

RF^i = cos 9 sin 0 cos 0 cos 0 + sin 0 sin 9 sin 0 — sin 0 cos 0 + cos 0 sin 9 sin 0 0 

— sin 9 sin 0 cos 9 cos 0 cos 9 Fq z 

(4.71) 

F Cz (sin 0 sin 0 + cos 0 sin d cos 0) 1 I" /i 1 

—>■ RFc = Fq z (— sin 0 cos 0 + cos 0sin 6 1 sin 0) = f 2 (4.72) 

[ F^ z cos 0 cos 9 J [/ 3 _ 

On doit identifier les entrees du sous-systenre de rotation. On a aussi deux cas : 


Utilisation de Tangle de lacet : Si on utilise Tangle de lacet, c’est-a-dire le 
cap du corps rigide se tourne vers la cible pendant le deplacement. Dans ce cas, 
Tangle de roulis est stabilise a zero degre : 0 ph 0 et on a : 

Fq z (sin 9 d cos 0 d ) 1 [" /i 

RF b c = Fq z (sin 9 d sin 0<0 = f 2 (4.73) 

F& z cos 9 d \ [f 3 _ 



tan0rf = F —y ip d = arctan 


tan 2 9 d = t/ 2 

J 3 


—> 9 d = arctan 


— V fi + fi + fi 


(4.74) 


Utilisation de Tangle de roulis et de tangage : Si on incline Tangle de roulis 

et Tangle de tangage pour le deplacement et si la direction du cap est fixee a une 
direction predefinie : 



F^ z (sin cj) d sin 0 + cos (j) d sin 9 d cos 0) 


' fi ' 

RF b c = 

F Cz ( — si n cos 0 + cos 0d si n @d sin 0) 

= 

/ 2 


Fq z cos 0 d cos 9 d 


h 


(4.75) 
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( tan 9 d = sinv 

-> \ tan tb d = /i^A'-/ 2 c 3 osV cog Qd (4.76) 

pb _ h - 

1 02 cos (f>d cos 9 d 

Car il doit prendre le temps pour que <f> —> (p d , 6 —> 6 d . Done, on utilise : 

( 9 d = arctan 

l fa = arctan cos 0^ (4.77) 

I FA = —^— 

V ^ z cos 4 > cos 6 

4.3 Application a un helicoptere a quatre rotors 

Parmi des robots volants, 1’helicoptere a quatre rotors a suscite un grand inte- 
ret en raison de sa manoeuvrabilite et la grande capacite de vol stationnaire. Ce 
vehicule VTOL a certains avantages par rapport a des helicopteres classiques : 
symetrie, simple a concevoir et a construire. Pourtant, le faible moment d’inertie 
de ce robot aerien le rend vulnerable aux accelerations angulaires importantes. 
Par consequent, pour de nombreuses missions potentielles une commande d’atti- 
tude eflicace est cruciale, car elle permet au vehicule de conserver une orientation 
desiree. Plusieurs techniques de commande ont ete utilisees pour la stabilisation 
d’attitude de l’helicoptere a quatre rotors, par exemple, le PID et le controle de 
LQ dans [19], la commande par mode glissant de [20] et l’approche backstepping 
appliquee dans [21]. 

Bien que les resultats presentes dans les travaux cites ci-dessus sont interessants, la 
bornitude de hentree de commande n’est pas consideree. Un aspect particulier de 
la dynamique d’helicoptere a quatre rotors est que le modele dynamique peut etre 
reduit a celui d’un corps rigide. Par consequent, l’approche proposee pour la stabi¬ 
lisation d’attitude et de position d’un corps rigide, est naturellement etendue a la 
stabilisation d’un mini helicoptere a quatre rotors. C’est pour cela que le quadriro- 
tor est un exemple choisi pour la simulation des lois de commande pour stabiliser 
un corps rigide dans l’espace a trois dimensions. D’une part, il est sous-actionne 
et il a deux degrees de liberte de translation non-controlables directement a cause 
du manque des actionneurs sur les axes x et y. D’autre part, le quadrirotor est 
le drone le plus utilise dans le domaine scientihque et civil grace a une structure 
mecanique simple et a une compromis de poussee/stabilisation tres eleve. Les tra¬ 
vaux sont aussi une contribution a la commande d’un helicoptere a quatre rotors. 
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Les equations dynamiques de Newton-Euler pour le deplacement d’un quadrirotor 
dans l’espace a trois dimensions sont les suivantes : 

f mv + mu x v = T B = W b + F £ , . 

\ Iu + u x (Ito) = Qb — Tc 

ou ; Fc = [ 0 0 Tc ] ;fc = [ Fci rC 2 Tc3 ] • 


Les lois de commande bornees proposees sont : 

f vnv s = W* + RF^ — mu x = — sat Mp {^E d E P + K p E^j 
| Jw = f c -wx(/a;) = -sat Me (k^Eq + K®E e ^ 

j RF£ = mu x v f - W f - sat-Mp (k%E p + K?E P ^ 
f c = ux (Iu) - sat Me (KfE e + K®E e ) 


(4.79) 


(4.80) 


Pour faciliter les lois de commande, nous considerons la partie composee de la force 
Coriolis (mu x yf) et le couple centrifuge (u x (Iu)) connne des perturbations au 
systeme. Enfin, les lois de commande suivantes sont utilisees : 


RF b c = -Wf - sat Mp (k£E p + K^Ep) = [ /, f 2 / 3 f 

Tc = —satM @ (KfEe + KpEe'j — [ Tea r C2 Tc3 ] 


Avec : 

• Wf — [ 0 0 —mg ] T 

• sat Mp (l<d Ep + Kp E^j = 

• satM e (yEfEe + KpE^j = 

Done : 


sat Mx (k d (x - x d ) + kp(x- x d )) 
sat My (k y d (y - y d ) + k%(y- y d )) 
sat Mz (k d (z - z d ) + k p (z — z d )) 

satM+ (k d (j> ~ <fid) + kf, (4> - (f) d )^ 

sat Mg [k° d (e - e d j + kp(9 — e d j) 

sat Mi (k% (ij> - ^ + k$ (^ - ij; d ) ) 


[h h h] T 


-sat Mx (k d (x - x d ) + k*(x- x d )) 
—sat My ( k\ (y - y d ) + k y (y - y d )) 
mg - sat Mz (k d (z - z d ) + k z (z - z d )) 


(4.82) 
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Suite aux resultats precedents, les lois de commande de position sont calculees 
partie de : 

( 9 d = arctan 

\ = arctan (^ /isin, y c ° s ^ cos (4-83) 

pb — _ h _ 

V ± Cz - cos <j) cos 9 


Finalement les lois de commande utilisees sont : 
• Pour la position : 


9 d = arctan 
c t> d = arctan 


-sat Mx (k d (x-x d )+k*(x-x d ))) cosip+(-sat M y (k d (y-yd)+kp(y-y d ))) sinV’ 
mg-sat Mz (k d {z-z d )+k z {z-z d )) 

-sat Mx (k% (x-x d )+k% (x-x d ))) sin ip-(-sat My (k y d (y-y d )+k y (y-y d ))) cos ip 

mg—sat Mz ( k%(z-z d )+k*(z-z d )) 

mg-sat Mz (k z d {z-z d )+k z (z-z d )) 
in — - 


cos 


cos (j) cos 6 


(4.84) 

Notons que les angles <p d , 9 d sont scalaires parce qu’ils sont calcules a partie 
des etats mesures a chaque boucle de commande, done <p d = 9 d = 0. 

• Pour l’attitude : 

Pci = - sat M + (k d <t> + k*{<f>- 0 d )) 

< r C 2 = -sat Mg (k d 9 + kp(9 - 9 d )^j (4.85) 

Tc 3 = - sat M * (ip - ip d ^ + kf (ip - V’d)) 


Remarque 4.2 : 

Notons que la loi de commande d’attitude ne depend pas de la connaissance des 
parametres du modele (I’inertie I , qui est difficile d’identifier). En consequence, elle 
est robuste aux incertitudes et erreurs du m.odele. Par contre, la loi de commande de 
position depend d’un parametre du modele, e’est la masse m. Meme si ce parametre 
est facile a identifier, parfois on estime incorrectement la masse et par consequent 
il nous induit une erreur statique dans la commande d’altitude. On suppose que la 
masse dans la loi de commande d’altitude est mesuree m! au lieu de m. Quand le 
drone est stabilise a une certaine altitude, on a : Tq = mg et <p> = 9 = 0°, i = z d = 
0 m/s . L’erreur statique est calculee a partir de Vequation : 

m'g — kp(z — zf) = mg ->• e z = (z - z d ) = — 9 f 0 (4.86) 

K~ 

Pour annuler cette erreur statique, on ajoute une partie integrate dans les lois de 
commande. Le terme integral accumule les erreurs a chaque iteration de la boucle 
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d’asservissement et reduit les petites erreurs statiques. Pourtant, si la sortie de la 
loi de commande est saturee pendant la periode de mise a jour precedente, Verreur 
actuelle n’est pas ajoutee a la valeur d’integration. C’est pour eviter une condition 
appelee "integrator-windup" qui se produit quand le terme integral continue a ac- 
cumuler des erreurs lorsque la sortie est saturee. 

Les forces et les couples sont calcules a partir des vitesses de rotation des action- 
neurs : 


/ r i \ 


/ 0 db 0 — db \ 


( s\\ 

r 2 


db 0 — db 0 


4 

r 3 


k —k k —k 


4 

V T > 


\ b b b b ) 


\4J 


Done la vitesse de rotation desiree pour chaque actionneur est obtenue au moyen 
de la transformation lineaire suivante : 
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(4.88) 


S d = N- 1 T 


(4.89) 


La figure 4.11 montre le schema-bloc de la strategic de commande. La loi de com¬ 
mande four nit un signal de commande Tc et Tc qui doit etre injecte au systeme 
propulseur. Ainsi, le systeme propulseur a l’objectif de generer le couple de com¬ 
mande qui doit s’exercer autour de chaque axe de helicoptere et la poussee totale. 
Par consequent, le couple de commande et la poussee totale de commande Tc 
doivent etre transformes en quatre signaux de vitesse de rotation desiree s d n avec 
n G {1,2, 3,4}, ou n represente le nombre d’actionneurs. Les commandes PI sont 
utilisees ponr l’objectif de donner les consignes PWM pour que le moteur n ait 
une vitesse de rotation s n egale a la vitesse desiree s^ donnee par l’equation (4.88). 
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Figure 4.3 — Schema de commande de l’helicoptere a quatre rotors 


4.4 Amelioration des lois de commande 


4.4.1 Autocorrection des gains de la loi de commande par 
les reseaux de neurones 

Les parametres de controle pour la commande PID doivent etre ajustes pour rea- 
liser un controle plus precis. II est evident que la methode de controle PID est 
limitee parce que les parametres k p , k d , ki sont constants. Cela signifie que la sor¬ 
tie de commande n’est pas adaptable et optimale dans tous les cas. Un controle 
intelligent adaptatif, qui a la capacite d’adaptation des parametres de controle 
pour minimiser l’erreur de position, est utilise. Avec la capacite d’apprentissage et 
d’adaptation des reseaux de neurones, le controleur peut resoudre ces problemes. 
Les parametres PID seront mis a jour de maniere adaptative et de fagon optimale 
afin de minimiser l’erreur par rapport aux perturbations externes et les erreurs 
internes. 

La figure 4.4 montre la structure de commande PID en utilisant les reseaux de 
neurones pour un degre de liberte. Le schema de reseau de neurones est illustre 
dans la figure 4.5. k p , kd, ki, e p , e* et e d sont le gain proportionnel, integral et derive, 
l’erreur du systeme entre la sortie desiree et reelle, l’integrale de l’erreur du sys- 
teme et la difference de l’erreur systeme, respectivement. Les reseaux de neurones 
sont entraines par l’algorithme de retro-propagation classique afin de minimiser 
l’erreur du systeme entre la sortie desiree et reelle. Dans la figure 4.5, l’entree de 
la function sigmoide Jm(x) dans la couche de sortie est : 


u{k ) = k p (k)e p (k ) + k^e^k) + k d (k)e d (k) 


( 4 . 90 ) 
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Figure 4.4 — Structure du controleur non lineaire PID en utilisant le reseau de neurones 



ou, 


' e p (k) = x(k) - x d (k) 

k 

< ajk) = y^e p (n)AT 

n= 1 


(4.91) 


e d (k) 


e p (k) - e p (k - 1) 

AT 


AT est le temps d’echantillonnage, k est la sequence discrete, x d (k) et x{k) sont 
des positions desiree et reelle. Pour regler les gains du regulateur PID, la methode 
de descente de gradient utilisant l’equation suivante a ete applic^uee : 


k p (k + 1) = kp(k) - 

ki(k + 1) = h(k) - Vi—7^ 
k d (k + 1) = k d (k) - 


(4.92) 


ou 7] p ,7]i,7] d sont les taux d’apprentissage determinant la vitesse de convergence et 
E(k) est une fonction de cout. Les parametres PID sont mis a jour a l’etape k afin 
de diminuer la fonction de cout. Puisque le but de la commande est de diminuer 
maximalement l’erreur e p {k ), par consequent nous choisissons la fonction de cout 
E(k) : 

E(k) = 5 Mk) - x(k)f 


(4.93) 



































4.4 Amelioration des lois de commande 


116 



Figure 4.6 - La function sigmoide fj\i{x) 


Selon [30], les parametres sont mis a jour : 

k p (k + 1) = k p (k) + r]p e p(k)e p (k) d ^^ 
ki(k + 1) = h{k) + r)ie p (k)ei{k 

ou 

k d {k + 1) = k d (k) + rj d e p (k)e d (k) d ^^ 

ou fnix) est une fonction sigmoide non-lineaire de saturation : 

M( 1-e-W) 

n _g g-2 x/M\ 


(4.94) 


(4.95) 


Notons que lim {/m (x)) = Msign(x), done M est la bornee de saturation. La 

IE—>-=b00 

figure 4.6 montre la forme de la fonction sigmoide /m(x) avec divers M. La fonc¬ 
tion sigmoide devient lineaire lorsque M tend vers l’ infi ni. 

Nous avons : 


df M {k) _ d (M( 1 - e ~ 2u ^/ M ) \ _ 4 e -2«(fc)/M 

du ~ du l (1 + e- 2u ^/ M ) J ~ (l + e - 2 u(fc)/M) 2 


(4.96) 


















117 


Loi de commande 


Done, les parametres PID sont regies a chaque boucle : 


k p (k + 1) 
h(k + 1) 
k d {k + 1) 


kp{kd) ~\~ f]p€p(k')6p(k') 


h(k ) + r)ie p (k)ei(k) 


4 e -2 u(k)/M 
(1 -f e -2 u(k)/My 
4 e -2 u(k)/M 


(1 _|_ e -2 u(k)/My 


k d (k) + rj d e p (k)e d (k) 


4 e -2 u(k)/M 
(1 + e -2 u{k)/My 


(4.97) 


Remarque 4.3 : 

On peut remarquer que, lorsque les taux d ’apprentissage r] sont mineures, la vitesse 
de convergence sera lente, ce qui diminuera Vadaptation du systeme. Mais si rj est 
grande, le reseau de neurones s’adapte plus rapidement mais la fonction de cout de- 
passe plus facilement la valeur minimum et les gains ne peuvent pas converger vers 
les valeurs necessaires. A fin d’ameliorer Valgorithme de descente de gradient, les 
taux d ’apprentissage rj sont egalement mis a jour de maniere adaptative suivant : 

r)(k + 1) = (1 + gsign ( E(k + 1) — E(k ))) rj(k) (4.98) 

oil g est un parametre positif. 


4.4.2 Generation de trajectoire 

Pour le cas d’une stabilisation a une certaine altitude, les entrees souhaitees de 
1’attitude sont les angles de roulis et de tangage a 0°, un angle de lacet desire 
constant et la position desiree. Pourtant, dans le cas ou le corps rigide doit realiser 
le deplacement pour un ou plusieurs degres de liberte, les controleurs demandent 
les prochaines entrees desirees pour ces degres de liberte. Un algorithme de gene¬ 
ration de trajectoire doit alors etre utilise pour le deplacement optimal. 

La generation de trajectoire peut ne pas etre utilisee en forgant simplement la 
position et l’attitude desiree finale. Cependant, quand la position et l’attitude de¬ 
siree finale est atteinte, la deceleration necessaire pour arreter le vehicule sera tres 
elevee. Puisque le drone ne peut pas realiser une deceleration si elevee car il vole 
dans l’environnement avec un frottement tres faible (l’air), il va depasser la posi¬ 
tion et Lattitude finale. Par consequent, le vehicule oscille autour de la position et 
Lattitude finale avant de Latteindre. Cette oscillation et deplacement ralentissent 
Lobjectif final et peut avoir des consequences dangereuses dans le cas d’une transi¬ 
tion ou d’un atterrissage autonome. De plus, si on donne tout de suite la position 
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desiree finale, l’erreur entre la position actuelle et la position finale sera tres elevee. 
La sortie de la loi de commande est done toujours saturee, ce qui peut causer des 
problemes pour les actionneurs et des accidents imprevus. 



Figure 4.7 — Trajectoire trapezoidal de la vitesse 



Figure 4.8 — Angle de tangage desire genere par l’algorithme de la generation de profil 
et son adaptation pendant la transition 


Afin d’avoir un deplacement doux vers la cible, nous avons besoin d’un algorithme 
pour generer la trajectoire pour controler la vitesse et l’acceleration. Line trajec¬ 
toire de vitesse lineaire par morceaux est mis en oeuvre : la vitesse est incrementee 
par une valeur d’acceleration constante jusqu’a ce que une vitesse maximale spe- 
cifiee soit atteinte. La vitesse maximale est maintenue pour un intervalle de temps 
necessaire et ensuite diminuee de la meme acceleration (deceleration) jusqu’a ce 
que la valeur zero de la vitesse soit atteinte (voir figure 4.7). La distance de depla¬ 
cement, la vitesse limite et la valeur d’acceleration sont necessaires pour executer le 
profil. A partir de ces donnees, le generateur de trajectoire va produire des courbes 
de vitesse en forme de trapeze pour un long mouvement et des courbes triangu- 
laires pour un court mouvement ou la vitesse maximale n’a pas ete atteinte. Le 
profil trapezoidal de la vitesse doit etre dans le profil de vitesse de z/ max ( la vitesse 
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physique maximale) et z> max (Tacceleration maximale). Les lois de commande de la 
position et de l’attitude vont forcer le corps rigide a suivre les nouvelles positions 
et attitudes generees par le generateur de trajectoire afin de realiser le deplacement 
vers la cible. 

Un exemple de la generation de trajectoire est montree sur la figure 4.8. Elle nous 
montre Tangle de tangage desire genere par le generateur de profil de mouvement 
et son adaptation dans la simulation lorsque l’avion convertible fait la transition 
du vol vertical au vol horizontal. L’angle de tangage ne depasse pas la position 
finale souhaitee (—2°), ce qui aurait pu provoquer des oscillations et au pire les 
dommages de l’avion apres la transition en atteignant Tangle de tangage final. 


4.4.3 Resultats de simulation 

Cette section presente les resultats de simulation des lois de commande proposees. 
Trois cas sont etudies : le premier est que le quadrirotor stabilise la position et 
Tattitude vers l’origine sans utiliser le generateur de trajectoire, le second cas avec 
le generateur de trajectoire et le troisieme avec une trajectoire simple d’helicoi’de. 
Dans les simulations, nous supposons qu’il n’y a pas de perturbation. Pour des 
raisons de securite, les couples et la poussee de commande sont bornes comme 
suit : 


G [0.5 mg, 1.5 mg] N 

r x ,r» G ('-0.4,0.4] N.m (4.99) 

r* e [—0.15, 0.15] N.m 


Sans la generation de trajectoire 

Le quadrirotor est stabilise a la position (6m, —4m, 3m) et a Tattitude (0°, 0°, 20°) 
vers Torigine. Sans utiliser la generation de trajectoire, le drone depasse la position 
finale (voir figure 4.9(a)) (plus de 2m sur Taxe x et lm sur Taxe y et z) et il perd 
done plus de 2 secondes pour stabiliser par rapport au cas oil on utilise la generation 
de trajectoire (figure 4.10(b)). 

La video de simulation dans le cas sans utiliser la generation de trajectoire est 
telechargee sur : 

http://www.youtube.com/watch?v=SlGsSaDQT2s 
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(a) Position du quadrirotor (b) Attitude du quadrirotor 

Figure 4.9 — Sans le generateur de trajectoire 


Avec la generation de trajectoire 

Le quadrirotor est stabilise a la position (6m, —4m, 3m) et a l’attitude (0°, 0°, 50°) 
vers l’origine. Les vitesses de la position ont une forme trapezoidale et diminuent 
vers zero quand le drone arrive a la cible (voir la figure 4.10(c)). Notons que les 
gains des lois de commande ne sont pas changes et ils sont les memes avec ceux 
dans le cas sans utiliser le generateur de trajectoire. L’adaptation est meilleure 
grace a la generation de trajectoire qui genere les positions desirees pour que le 
drone puisse les suivre vers la cible sans le depassement. 

La video de simulation dans le cas utilisant la generation de trajectoire est tele¬ 
chargee sur : 

http://www.youtube.com/watch?v=i3kWxD41ZrU 


Suivi d’une trajectoire simple d’helicoide 

Cette section presente les resultats de simulation ou le quadrirotor suit une tra¬ 
jectoire d’helicoide cree par les equations suivantes : 

Xd = 3cos (f t) [m] 

y d = 2 sin (ft) [m] (4.100) 

z d = O.lt [m] 

Les figures (4.12) et (4.13) nous montrent que le drone peut suivre la trajectoire 
desiree d’helicoide en inclinant bangle de tangage et bangle de roulis pendant le 
mouvement. 





























121 


Loi de commande 




(a) La position desiree generee par le genera- (b) La position actuelle du quadrirotor 
teur de trajectoire 




(c) La vitesse de translation 


(d) La vitesse de rotation 




(e) L’attitude desiree generee par la loi de com¬ 
mande de position (<j), 6 ) et le generateur de tra¬ 
jectoire (tj}) 


(f) L’attitude actuelle du quadrirotor 


Figure 4.10 - Le deplacement de quadrirotor en utilisant le generateur de la trajectoire 
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Figure 4.11 - Les couples et forces dans le cas ou on utilise le generateur de trajectoire 

La video de simulation d’un suivi de trajectoire d’helicoide est telechargee sur : 

http://www.youtube.com/watch?v=eY5Yze59dlM 


Simulation du decollage et atterrissage autonome 

La figure 4.14 montre que la performance de la loi de commande est amelioree 
en utilisant le reseau de neurones lorsque le vehicule realise 3 fois le decollage et 
Latterrissage de Om a 2m. La courbe (1) est la position desiree par le generateur 
de trajectoire et le second (2) est 1’adaptation de la commande d’altitude lorsque 
les gains ont ete choisis correctement. La courbe (3) est l’adaptation quand les 
gains ont ete choisis incorrectement. L’adaptation est presque acceptable grace au 
generateur de trajectoire, mais il y a des erreurs statiques au debut. Cependant, 
ces erreurs statiques sont rapidement eliminees par la partie integrate dans la loi 
de commande. 

Dans la courbe (4), le decalage de la loi de commande d’altitude mg a ete modifie 
a 0.6mg (c’est le fait d’estimer les parametres du systeme inexactement) et les 
gains restent incorrects comme dans le cas (3). Les parametres de controle dans 
le test (5) sont similaires a ceux de (4), mais ici, des bruits de capteur et des per- 
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- current position 

. desired position 



Figure 4.12 — La position du quadrirotor en suivant une trajectoire d’helicoide 




(a) La position desiree (b) L’attitude desiree 




(c) La position actuelle (d) L’attitude actuelle 


Figure 4.13 — La position et l’attitude de quadrirotor en suivant la trajectoire d’helicoide 

































































4.4 Amelioration des lois de commande 


124 



Figure 4.14 — Amelioration de la loi de commande par les reseaux de neurones 




Figure 4.15 — L’evolution des taux d’apprentissage r) 


turbations externes (moyenne nulle, des perturbations bornees avec des deviations 
standard de a z = 0.03m et a F b = AN) ont ete ajoutes a la mesure de l’altitude 
et respectivement. Les conditions de l’essai (6) sont les memes que (5) (des 
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gains et parametres incorrects de la loi de commande, des perturbations d’entree 
et de sortie ), mais les gains sont maintenant corriges de maniere adaptative par 
les reseaux de neurones. La figure 4.15 montre le changement adaptatif des taux 
d’apprentissage rj pendant le test. 

En conclusion, grace aux deux methodes ajoutes : les reseaux de neurone et le 
generateur de trajectoire, la performance des lois commande, qui sont simples, 
realisables en realite, est maintenant amelioree. Les lois de commande sont adap- 
tatives et robustes par rapport a des perturbations externes et des erreurs internes 
de la loi de commande. Particulierement, le generateur de trajectoire permet au 
corps rigide de se deplacer a la position finale sans depassement, ce qui est tres 
important dans le cas de transition et pour 1’atterrissage. 


4.5 Loi de commande pour l’attitude utilisant le 
vecteur quaternion 


Les lois de commande d’attitude proposees precedemment est utilisable seulement 
si l’attitude est limitee. Lorsque Tangle de tangage du corps rigide s’approche de 
±90°, ces lois de commande ne sont pas utilisables a cause du probleme de singula¬ 
rity de la presentation de Tattitude des angles d’Euler. Afin d’eviter la singularity 
lorsque l’avion convertible vole verticalement, une loi de commande bornee utili¬ 
sant l’erreur d’attitude estimee par l’erreur de quaternion est utilisee. 


D’abord, nous definissions une fonction scalaire H(qo) : 


Definition 4.2 : H(qo) est une fonction non-negative sur Tintervalle go G [—1,1] 
et il s’annule seulement a go = — 1 et/ou go = 1- -S’(go) satisfait aussi la condition 

de Lipschitz sur Tintervalle go £ [—1,1], ce qui signifie 


dH 


dqo 


< 


7- 


Theoreme 4.6 : Pour toutes les fonctions H(qo) definies ci-dessus, la loi de 
commande bornee suivante permet de faire converger globalement et asympto- 
tiquement le corps rigide vers Tattitude desiree : 


x, , x dH 

To; scitMx. ( AeCUe p — q ei 

oq eo 


T. , X dH 

Pj/ I XyUJy P~rpZ 9e5 


dq, 


eo 


(4.101) 

(4.102) 
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r* SO't'Mz ( P o Qes 


(4.103) 


ou A i, p > 0; Mi ^ 3/ry; i G {x,y,z}; q e represente l’erreur d’attitude entre l’atti¬ 
tude actuelle du corps et l’attitude desiree. 


Preuve : Nous allons definir la fonction de Lyapunov suivante : 

\ 

V = -co T Jco + 2 nH(q eo ) ^ 0 

ou k > 0. En differentiant V par rapport a temps, on a : 

• _, T dH 

V = co Jco + 2 k —— q en 

dq eo 

. -i dH f 1 _fp\ 

V = co 1 —co x Jco + T + 2k —— ( —qi co 
L J dq eo V 2 J 

T> -*T-p -T - 

V = co I — k—— q e co 
dq eo 


(4.104) 


(4.105) 


(4.106) 


(4.107) 


• ^ dH ^ dH ^ dH 

V (jJ x l x ^7^ Qei^x ~ y ^77 Qe 2 ^y~^~^z^z ^77 Qe^^z (4.108) 

dq e o %„ &Qeo 


V est la somme de trois tenues (Vi, V^, V 3 ). Nous analysons d’abord le terme : 


- r - dH 

V\ CJ X L X K, Q e 1 Cd x 

dq eo 

L’application de la loi de commande bornee (4.101), on a : 

■ L dH 1 dH 

\ 1 A x 00 x P~,n qei ^7 Qei^x 

dqe 0 J uq eo 

Nous supposons que : 


(4.109) 


(4.110) 


\Ku x \ > 2p7 ^ 2 p 


77 \C 0 T \ > 


(4.111) 


nous avons ensuite : \ x co x — p-§f~q ei ^ pq + £ (avec £ > 0) et [\ x co x — p-§^-q e ^j 
a le meme signe que celui de co x . Done : 


• l\ dH 1 dH 

m u x sat Mx a x lo x p— q ei q ei co x 

l dq eo J dq eo 

^ - \u x \ sat Mx [p 7 + e] + ^7 1^1 


(4.112) 
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En choisissant : 


k 7 < min (M x , py + e) 77 n < 


min (M x , py + e) 


7 


(4.113) 


nous pouvons assurer la decroissance de V\ parce que V\ < 0. En consequence, u x 


entre en $ 7 , = cm. : |u;J < 


\ 

At; J 


dans un temps fini t\ et reste dedans par la suite. 


Dans ce cas, [\ x u x — p^-q e ^J G [—3py,3py]. En choisissant M x qui satisfait 
la condition suivante : M x ^ 3py. Ensuite, pour le temps f 2 > ti, 1’argument 


satM x sera bornee comme suit : 
une condition pour k de (4.113) : 


P a n Qe 1 


dq, 


e 0 


P7 + £ £ 

K < - = p + - 

7 7 


^ 3py ^ M x et done nous avons 

(4.114) 


En consequence, satM x opere dans une region lineaire : 

r - \ dH 


dq, 


et : 


eo 


OH \ dH 

~ K dT qe ^ 

UL Le 0 / UL Le 0 

• dH dH 

Vi = ~\ x u x + p -7— q ei u x ~ k- — q ei u x 


Vl UJ X ( \ x oj x p 


dq, 


eo 


dq, 


(4.115) 

(4.116) 

(4.117) 


eo 


En choisissant k = p qui satisfait la condition (4.114), on a : V\ — —A x co x ^ 0. Si 
les memes arguments sont appliques a V 2 et V :i , V devient : 


V — Vi + V2 + V3 — — (A x tu^ + A y Uy + A 2 tu^) ^ 0 


(4.118) 


Ahn de completer la preuve, le principe d’invariance de LaSalle est invoque. Toutes 
les trajectoires convergent vers le plus grand ensemble invariant dans D = 

=, u) : V — 0 } = {(q e ,Lo) : u = 0}. Dans cet ensemble invariant : 


(4.119) 


Ju = -u x Juj + r = p-— [ q ei q e2 q e3 ] = 0 

OQeo 


Ainsi, D est reduit a l’origine : q e — uj — 0. Cela finit la preuve de la stability 
asymptotique du systeme en boucle fermee. 


Le tableau 4.1 presente quelques fonctions H(q 0 ) avec leurs points d’equilibre 
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Table 4.1 — Exemples de fonction H(qo) 


H{qo) 

Points d’equilibre asymptotiquement stable 

Points d’equilibre instable 

1 - 90 

90 = 1 

90 = -1 

1 + 90 

9o = -1 

9o = 1 

1 - \qo\ 

90 = ±1 


i - M” +1 

9o = ±1 


oos(T) 

9o = ±1 


cos n (4^) 

90 = ±1 



asymptotiquement stable et les points d’equilibre instable [24], Afin d’optimiser et 
d’ameliorer les performances de la loi de commande utilisant le vecteur de quater¬ 
nion, les parametres positifs A* et pi ont ete calcules selon la methode de placement 
de poles de controle bornee proposee dans [18]. 

Propriety 4.2 (Contrainte sur la vitesse angulaire) La loi de commande (4-101) - 
(4-103) stabilise globalement et asymptotiquement Vattitude du corps rigide a I’ori- 
gine. De plus, si les parametres p , A, sont choisis tel que la condition suivant est 
satisfait : 

P ^ma Xj 

A* 2 q 

oil ta maxi la vitesse angulaire maximale 2 que le gyrometre peut mesurer lors d’une 
manoeuvre autour de Vaxe de rotation. La contrainte \uji\ ^ ta maXi reste satisfaite a 
partir d’un certain t x > t 0 si |cuj(f 0 )| > w m axi et pour tout t si \u>i(t 0 )\ < Wmar 

Preuve : Supposons que les conditions initiates du systeme sont telles que a l’ins- 
tant to, u>i appartient a 1’ensemble = |c Oi : |a^*| < ^ j. A partir de la preuve 
du theoreme, nous pouvons observer que si uii appartient a <3>j, ta, va rester dans 
pour tout le temps t > t 0 . Par consequent, si nous choisissons : + = , nous 

pouvons conclure que |ta,| < ^ = ca maXi , Vf > t 0 . 

Par contre en realite, des perturbations externes exercent sur le corps rigide et 

2. La Centrale Inertielle 3DM-GX3-25 utilisee dans le drone a une plage de mesure de vitesse 
angulaire de ±300°/s. 
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dans certaines situations la vitesse angulaire du corps rigide depasse la plage de 
mesure du gyrometre : 

M > Wma Xi (4.121) 

Dans ce cas, les gyrometres seront satures et la loi de commande suivante doit 
converger la vitesse de rotation c <7 a l’ensemble <t>j. : 


r i SCltMi ( ^ i 


■'ma Xj 


d H 

p W qe ' 


En analysant V) de l’equation 4.110, nous avons : 
V) (jJiSdljMi 


dH 

P ~7\ Qei 


dq, 


eo 


eo 


dH 

- k- — q ei Ui 


dq, 


eo 


(4.122) 


(4.123) 


Puis que \u maXi = 2py, 


dH 


dqo 

A iOJ, 


< 7 , il suit : 


dH 


max 7 ; P o Qe 
UQeo 


> pi = 


Aitu, 




(4.124) 


Du au fait que le gyrometre est sature nous avons : sign(uji ) = sign (cu maXi ). En 
suite, V{ devient : 

t V , II , AjO; max , , , 

Vi < - \uJi\sat Mi —-— + «7 Kl 


Si nous prenons : 


K'y < min ( Mi , 


AjCu n 


min (Mi 


<3- k < 


AjCOn 


L ii 2 


7 


(4.125) 


(4.126) 


il est possible d’assurer que V) < 0. Done, il existe un temps t\ afin que la vitesse 
de rotation u t converge vers l’ensemble et reste dedans. En consequence le 
gyrometre n’est plus sature et le systeme evolue de la meme fagon que le systeme 
nominal. 


4.6 Avion Convertible 

4.6.1 Vol longitudinal 

Dans la dynamique longitudinale (voir figure 4.16), l’avion convertible effectue 
les phases de vol suivantes et la transition entre les controleurs de chaque phase 
depend de bangle de tangage : 
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Figure 4.16 - Vol longitudinal de l’avion convertible 


• Le vol vertical : L’avion est dans la position verticale, c’est-a-dire Tangle 
de tangage est —90°, et done il a la capacite de decoder, voler et atterrir 
verticalement coniine un helicoptere. Les parties utilisees pour commander 
l’avion sont les deflecteurs de la surface de control des ailes immerges dans le 
flux d’air cree par l’helice. Ils genereront les couples pour controler l’avion. 

• Le vol horizontal : dans cette phase, l’avion convertible vole coniine un 
avion conventionnel. Pourtant, en mode avion certains parametres comme 
Tangle d’attaque, la vitesse de Tavion ne peuvent pas etre identifiee precise- 
ment dans la realite. Alors, on rencontre une difficulty pour developper une 
loi de commande pour stabiliser l’avion en mode horizontal parce que les 
parametres du systeme dans cette phase sont inconnus. En plus, les travaux 
pour debugger les codes source dans cette phase sont vraiment difficiles a 
cause de la vitesse elevee de l’avion. Dans la simulation, les lois de com- 
mande d’attitude dans le vol vertical sont utilisees pour stabiliser Tattitude 
pour le vol horizontal, mais en realite l’avion est pilote manuellement. 

• La transition : La phase la plus importante de l’avion convertible est la 
transition entre le mode vertical et le mode horizontal. Sachant que la tran¬ 
sition de l’horizontale a la verticale est inverse a la transition de la verticale 
a Thorizontale, on va construire une loi de commande unique pour ces deux 
transitions, mais avec deux generateurs de trajectoire de Tangle de tangage 
desiree differentes. Pour que le vehicule puisse faire la transition doucement, 
une trajectoire trapezoidale de la vitesse de Tangle de tangage est construite : 
au debut, la vitesse accelere de zero a une certaine valeur constante et puis 
diminue vers zero a la fin de cette phase de transition. Notons que pour la 
transition du mode vertical a au mode horizontal, Tangle de tangage varie 
progressivement de —90° a 0°, tandis que pour la transition de Thorizontale 
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a la verticale, il varie de 0° a —90°. 

• Etablissement de la vitesse de l’avion : Pendant la phase de transition 
de la verticale a l’horizontale, le vehicule doit accelerer pour atteindre line 
certaine vitesse pour que le vehicule puisse avoir une portance generee par 
l’effet aerodynamique plus elevee ou egale que son poids. 

• La deceleration : l’inverse de la phase d’acceleration precedente est la phase 
de deceleration de l’avion pendant la transition de Fhorizontale a la verticale. 
En fait, pour realiser cette phase les vitesses de rotation des deux moteurs 
sont reduites avant la transition a la verticale. Elies sont ensuite augmentees 
apres la transition aux valeurs auxquelles elles contribuent une poussee totale 
egale du poids de drone. La phase qui suit cette phase est le vol vertical pour 
commander l’avion comme un helicoptere. 


4.6.2 Les lois de commande 

Suite au modele longitudinal decrit precedemment, le drone n’est qu’autonome 
(ou semi-autonome) dans la phase de vol vertical et de la transition. Dans le vol 
horizontal, le drone est pilote manuellement et les lois de commande ne sont pas 
developpees a cause de la difficulte de debugger le programme embarque avec les 
lois de commande car havion fonctionne avec une vitesse elevee. En fait, dans la 
simulation les lois de commande developpees pour le vol vertical sont encore utili- 
sees pour commander l’avion pour le vol horizontal mats nous avons decide de ne 
pas les developper dans la realite a cause de la securite (le drone peut causer des 
degats importants s’il fonctionne mal en mode horizontal). 

Les lois de commande sont developpees pour commander quatre degres de liberte. II 
s’agit de l’altitude et l’attitude (l’angle de roulis, de tangage et de lacet). Puisqu’en 
vol vertical, bangle de tangage s’approche de 90°, la loi de commande d’attitude 
presentee par le vecteur quaternion est utilisee (voir section 4.5). En fait, les lois 
de commande sont les suivantes : 

• Pour l’altitude : La loi de commande de baltitude est utilisee lorsque 
l’avion vole a la verticale, mais elle n’est pas utilisee dans la phase de tran¬ 
sition et le vol horizontal. En extrayant les equations de la dynamique de 
translation sur l’axe z de Ef et en considerant que les forces de Coriolis et 
centrifuges sont des perturbations au systeme, nous avons : 

rriz = T cos — d'j cos <f> — mg 


(4.127) 
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oil T est la force de poussee totale des deux moteurs. Une loi de commande 
bornee simple PID garantit la convergence de l’altitude a F altitude desiree : 


T C = T 1 + T 2 


—sat T ( k p (z - z d ) + k d z + k { f (z - z d )) + mg 
cos (l — 6) cos (f) 


(4.128) 


• Pour l’attitude : Dans les lois de commande d’attitude en utilisant le 
vecteur de quaternion de la section 4.5, puisque |g 0 1 ^ 1 et pour simplifier 
la commande embarquee en temps reel, la function H(q 0 ) = 1 — |g 0 | dans le 
tableau 4.1 a ete choisi. Sa derivee est = —sign(q 0 ). Notons que cette 
fonction n’a pas de points d’equilibre instables, et la loi de commande de 
l’attitude est donnee par (coniine dans [22]) : 


r c x = 

—sat Mx (Kuj x + psign (q eo ) q ei ) 

(4.129) 

V Cy = 

-sat My (Ayta^ + psign (q eo ) q e2 ) 

(4.130) 


-sat Mz (A z uj z + psign (q eo ) q e3 ) 

(4.131) 


La function sign(q eo ) assure qu’entre les deux rotations des angles (3 e et 
27 t — /3 e , le plus petit angle sera choisi. 

A partir des equations des forces (3.118) et des couples en vol vertical (3.119), 
nous devons retrouver des valeurs des actionneurs. II s’agit des vitesses de rotation 
des deux moteurs (^ 1 , 52 ) et les angles de deflexion des ailerons et du gouvernail 
(^ai, fi ar , dr)■ Nous avons les equations suivantes : 

t c = b(4 + 4) 

To, = bl ?f^(s 2 A l - S lS„)+k(sl-st) 

Am 

_ bl s S°C[ 

1 C y ~ -1+ S 2 0 ar ) 

Am 

\ S^CfsjSr + bl m (sl - s?) si S r > 0 
L £ S r C F4Sr + bl m (sl - sj) Si5 r < 0 

Nous avons 5 racines mais il n’y a que 4 equations. Pour resoudre ces equations, 
nous pouvons considerer que l’angle de deflexion du gouvernail 5 r est proportionnel 
avec le couple de lacet T Cz , c’est a dire : 

■ ^ S* r Cfsl5 r si r Cz > 0 
_ £s°Cfsl5 r St T Cz < 0 


(4.132) 

(4.133) 

(4.134) 

(4.135) 


k r T Cz — 


(4.136) 
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et 

k t T Cz = U m (si - s?) (4.137) 

ou k r et k t sont des taux proportionnels : k r + k t — 1. Par consequent, les vitesses 
de rotation des deux moteurs sont calcules par : 


kt r Cz = bl 


m °1) 


So = 




T c = b (sj + s 2 2 ) 
et 1’angle du gouvernail d r est done donne par : 


st = 


1 { k t T Cz , Tg_ 

2 l blm ' b 

1 ( Tc _ ktTcr. 

2 \ b bl m 


(4.138) 


3 r = 


k r r 


V 1 c z 


r r <5 V b blm ) 


k r T r 


bl s Qsf^P ( k t F C z ■ T C 
2 Am r $ \ blm b 


) 


si T Cz > 0 


si T Cz < 0 


(4.139) 


Les angles de deflexion des ailerons droit d ar et gauche S a i sont calcules a partir 
des equations : 


■p _ bl m S S a C 6 


(sldal - sld ar ) + k (s| - sf) 


r c „ = (° 2 A, + 4 U 


(4.140) 


Am 


2 bSgC^sl \ Is 


r c y r Cx -fc(s|-s?)^ 


TZ 2 l /-1-i- — °al 
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(4.141) 


2 bSiCfsl \ i 


r Cy _ r C X -k(s%sj) \ _ 
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k t T C z 
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^ C y , ^ C x b i 

I l 




A — _4m_ 

“ r «Cf(«a + T 0 ) 


In 


kk+Ts~i 

T Cy 


(4.142) 


Le schema-bloc de la strategie de commande est montre dans la Figure 4.17. La 
loi de commande fournit un signal de commande Tc et Tc qui doit etre injecte 
au systeme. Par consequent, le couple de commande Tc et la poussee totale de 
commande Tc sont transformes en deux signaux de vitesse de rotation desiree 
des deux moteurs et les angles de deflexion des ailerons et du gouvernail. Les 
connnandes PI sont utilisees pour l’objectif de donner les consignes PWM pour 
que le moteur n ait une vitesse de rotation s n egale a la vitesse desiree . 
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Figure 4.17 — Schema de commande de l’avion convertible 

4.6.3 Simulation de vol longitudinal 



Figure 4.18 - Programme Simulink de simulation de l’avion convertible 

Dans cette section, quelqnes resultats de simulation a partir des simulateurs Matlab- 
Simulink / MSC. visualNastran J^D (Figure 4.18) sont presentes pour montrer la 
performance de la strategic de commande. Dans la simulation, nous supposons 
qu’il n’y a pas de perturbation. Selon les logiques de la transition (1 = transition 
vers le vol horizontal, 0 = transition vers le vol vertical) et du vol vertical (1 = 































































































































































135 


Loi de commande 



Figure 4.19 - Chemin du vol longitudinal 


decollage, 0 = atterrissage), nous pouvons diviser les phases de vol comme suit : 

(1) . Decollage (depuis le debut jusqu’a, la seconde j) : le drone decode automa- 
tiquement jusqu’a 2m (figure 4.19), la vitesse d’altitude a la forme trapezoidale, il 
accelere et diminue la vitesse quand le vehicule atteint 2m (figure 4.20). 

(2) . Transition vers le vol horizontal (de la seconde j a 5,5) : l’avion fait la tran¬ 
sition a la seconde 4, Tangle de tangage varie de —90° jusqu’a —2°. Pendant la 
transition, le generateur de trajectoire cree Tangle de tangage desire qui varie pro- 
gressivement de —90° a —2° (transition vers le vol horizontal), de —2° jusqu’a —90° 
(transition vers le vol vertical). Les lois de commande d’attitude force le drone a 
suivre la trajectoire generee pour avoir une transition stable. 

(3) . Vol horizontal (entre les temps 5,5s et 11s) : Apres la transition, le drone 
accelere pour atteindre une certaine vitesse. Sa vitesse v x augmente tres vite et 
converge vers une valeur constante a cause de la trainee de frottement. A ce mo¬ 
ment, l’acceleration longitudinale du drone est egale zero. Sur la figure 4.23, on 
peut constater que Tangle d’attaque converge aussi vers une valeur constante. Mal- 
gre cette valeur petite (environ 1°), elle est suffisante pour creer la portance grace 
a la vitesse tres elevee de Tavion. II vole horizontalement environ 5s et puis fait la 
transition vers la verticale a la seconde 11 (la periode (4)) . Apres la transition, le 
drone stabilise son altitude dans la periode (5) et ensuite il descend et atterrit a 
la seconde 21 (la periode (6)). 

Notons que la loi de commande de l’altitude est deactivee pendant le vol horizontal 
et pendant la phase de transition. Le profil de la vitesse de l’altitude v z a la forme 
trapezoidale dans les periodes de decollage (1) et d’atterrissage (6) tant que le profil 
de la vitesse de rotation de tangage c o y a la forme d’un triangle dans les phases 
de transition (les periodes (2) et (4))(figure 4.20). Il n’y a pas de depassement 
dans les deux cas. Sur le chemin du vol longitudinal (figure 4.19), le drone perd 
son altitude dans deux cas : le premier apres la transition a l’horizontale parce 
que la vitesse mineure du drone ne permet pas de creer la portance necessaire; la 
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Figure 4.20 — Les variables d’etat de l’avion pendant la periode de vol 



Figure 4.21 — Forces aerodynamiques dans le repere Ef, 
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Figure 4.22 - Vitesse de chaque rotor 
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Figure 4.23 - Angle d’attaque 



Figure 4.24 — Transition vers le vol horizontal dependant des vitesses de moteur (si 

«2 = s) 
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deuxieme avant la transition a la verticale a cause de la deceleration. Le programme 
de commande de la simulation a ete developpe par le code C embarque (dans le 
bloc S-function Builder de Simulink). Ainsi, ces codes sont importes directement 
dans le programme de controle du microcontroleur. 


Transition vers le vol horizontal dependant des vitesses des moteurs 

La figure 4.24 montre la phase de transition dependant des vitesses des moteurs. 
On peut constater que si les vitesses de deux moteurs sont de Si = s 2 = 370 rad/s 
l’avion va tomber, si les vitesses de deux moteurs sont de Si = S 2 = 380 rad/s 
le drone perd une altitude de plus d’un metre pendant la transition. La vitesse 
a laquelle le vehicule fait la transition avec succes sans perdre d’altitude est de 
plus de 385 rad/s. En conclusion, la transition vers le vol horizontal est echouee 
si la poussee initiate avant la transition n’est pas suffisante, dans ce cas l’avion 
ne peut pas accelerer immediatement vers une vitesse a laquelle elle permet de 
creer la portance necessaire. En realite, cette vitesse est choisie par Toperateur en 
changeant le gaz et en regardant le comportement du drone pendant la transition. 
Notons dans tous les cas de simulation, Tangle de tangage varie progressivement 
de -90° a -2°. 


4.7 Avion basculant a trois rotors 

4.7.1 Strategie de commande 

Le principe du vol longitudinal de l’avion basculant a trois rotors est similaire a 
1 ’avion convertible mais ici les deux moteurs d’avant basculent de la verticale a 
Thorizontale et inversement au lieu du fuselage pendant les phases de transition. 
Le vol longitudinal est montre dans la figure 4.25. Notons que dans la phase de 
vol verticale, le drone a le comportement de Thelicoptere tri-rotor. 

Dans la cadre de la loi de commande, les angles d’Euler sont utilisees pour la 
representation d’attitude grace a leur simplicity. En fait, on utilise les lois de 
commande suivantes : 

• Pour l’altitude : Identiquement au cas de Tavion convertible, la loi de 
commande de l’altitude n’est utilisee que lorsque le drone vole a la verticale, 
mais elle n’est pas utilisee dans la phase de transition et le vol horizontal. En 
extrayant les equations de la dynamique de translation sur l’axe z de Ef et 
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Figure 4.25 — Vol longitudinal de l’avion basculant a trois rotors 


en considerant que les forces de Coriolis et centrifuges sont des perturbations 
au systeme, nous avons : 

mz — Tcosdcoscj) — mg (4.143) 

ou T est la force de poussee totale des trois moteurs. Une loi de commande 
bornee simple PID garantit la convergence de l’altitude vers l’altitude desi- 
ree : 


T = 


sat T ( kp{z d - z) + k d (z d - z) + k z f (z d - z)) + mg 
cos 9 cos 0 


(4.144) 


T 


sat T (k z e z + kfef + k z d e z d ) + mg 
cos 9 cos 0 


(4.145) 


• Pour Vattitude : La commande de l’attitude est utilisee pour stabiliser 
le drone quand il decode et atterrit verticalement. Pourtant, dans la phase 
de transition les deux moteurs avant doivent tourner avec la meme vitesse. 
Done seulement Tangle de tangage est stabilise dans cette phase de transi¬ 
tion, les mouvements de roulis et lacet sont commandes manuellement en 
changeant des ailerons et le gouvernail. La commande d’attitude PID bornee 
est utilisee : 

Tj = sat Fj (k 3 p e 3 v + k\e\ + k 3 d e 3 d ) (4.146) 

avec j = 0, 9,0 et Tj sont la saturation de la commande de couples autour 
des axes x, y, z. 


L’angle de basculement desire est genere par le generateur de trajectoire pour qu’il 
puisse varier progressivement de 90° a 0° et inversement. En fait, dans le cas de la 
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transition vers le vol en croisiere, si on donne tout de suite l’angle de basculement 
final (5 = 0°), le drone change immediatement vers le vol horizontal sans passer 
la phase d’acceleration. Le vehicule tombera certainement parce que sa vitesse ini- 
tiale de zero ne suffit pas pour creer la force de portance necessaire. Done, Tangle 
de basculement doit etre diminue lentement vers 0°. Pourtant dans la transition 
vers le vol stationnaire, Tangle de basculement final (6 = 90°) peut etre donne 
immediatement pour que le drone change rapidement de l’etat instable a l’etat 
stable (stabilise par la loi de commande en vol vertical). 

Pendant la transition vers le vol horizontal, quand les deux rotors d’avant com- 
mencent a basculer progressivement, la force des moteurs sur Taxe x apparait. Par 
consequent, la vitesse de Tavion augmente. La portance aerodynamique apparait, 
pourtant on ne peut pas mesurer ou estimer cette force. Done, la loi de commande 
de Tangle de tangage est utilisee. L’equation du couple de tangage est : 


^c e — T 3 I 2 — ( T'i + T 2 ) l\ sin 5 


(4.147) 


ou Tc g est le couple de tangage deduite de la loi de commande d’attitude pour 
stabiliser Tavion, T\ = T 2 est la force de deux moteurs d’avant donnes par l’ope- 
rateur via le recepteur R/C . Par cette methode, l’operateur peut varier la vitesse 
du drone pendant la transition. Sachant que Tangle de basculement 5 diminue de 
90° a 0°, la vitesse de rotation de troisieme moteur doit varier d’apres Tequation : 


rji ^c 6 + (^1 + ^2) h shi 8 

- - 7 - 


(4.148) 


Quand la force des moteurs sur Taxe z diminue (T 3 et (7\ + T 2 ) l\ sin 5), la force sur 
Taxe x augmente. Par consequent, la portance aerodynamique F\ va maintenant 
augmenter pour compenser la reduction de la force F j, grace a vitesse elevee du 
drone. En conclusion, grace a la generation de trajectoire de Tangle de bascule¬ 
ment la phase d’acceleration est realisee pour que le drone puisse avoir une vitesse 
initiate avant la diminution totale des forces des moteurs sur Taxe z. 


Le schema-bloc de la strategic de commande est montre dans la figure 4.26. La loi 
de commande four nit un signal de commande Tc et Tc qui doit etre injecte au sys- 
teme. Le couple de commande Tc et la poussee totale de commande Tc sont done 
transformes en deux signaux de vitesse de rotation desiree des deux moteurs et les 
angles de deflexion des ailerons, de Televateur et du gouvernail. Les commandes 
PI sont utilisees pour Tobjectif de donner les consignes PWM pour que le moteur 
n ait une vitesse de rotation s n egale a la vitesse desiree s^. 
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Figure 4.26 — Schema de commande de l’avion basculant tri-rotor 

4.7.2 Simulation de vol longitudinal 

Dans cette section, les resultats de simulation a partir des simulateurs Matlab- 
Simulink/MSC. visualNastran 4D (figure 4.27) sont presentes pour demontrer la 
performance de la strategic de controle. Dans cette simulation, on suppose qu’il 
n’y a pas de perturbation et que l’avion est stable dans la dynamique laterale et 
la reaction du servomoteur de Tangle de basculement est beaucoup plus rapide 
que la dynamique du vehicule. Cela signifie que seule la dynamique longitudinale 
est examinee et le servomoteur de Tangle de basculement s’adapte immediatement 
par rapport a Tangle de basculement desire. Selon les logiques de la transition (1 
= transition vers le vol en croisiere, 0 = transition vers le vol stationnaire) et du 
vol vertical (1 = decollage, 0 = atterrissage), nous pouvons diviser les phases de 
vol comme suit : 

(1) . Decollage (depuis le debut a la seconde 3) : Initialement, le drone decolle 
automatiquement jusqu’a 2m, la vitesse d’altitude a une forme trapezoidale, il 
accelere et reduit la vitesse lorsque le vehicule atteint 2m. 

(2) . Transition vers le vol en croisiere (de la seconde 3.5 a 8.5) : Tavion fait la 
transition a la seconde 3.5, Tangle de basculement varie de 90° a 0° et Tangle de 
tangage incline 2° afin d’avoir un angle d’attaque eleve. Cet angle de tangage peut 
etre augmente plus , mais 2° est suffisant pour prouver la performance de la stra¬ 
tegic de commande. Le drone perd Taltitude au debut de la transition parce que 
sa vitesse mineure ne peut pas permettre de creer la force de portance suffisante. 
En fait, on peut augmenter la vitesse de rotation des deux moteurs d’avant pour 
que Tavion puisse eviter de perdre de Taltitude pendant la transition mais ici ils 
sont encore tenus a la valeur (349 rad/s) a laquelle ils contribuent a egaler le poids 
de Tavion dans le vol vertical (figure 4.31). Pendant la transition, la vitesse de 
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Figure 4.27 — Programme Simulink de simulation de l’avion basculant tri-rotor 


rotation du moteur arriere se reduit a 0 rad/s d’apres l’equation (4.148) et l’avion 
accelere grace a l’augmentation de la poussee (voir figure 4.30). La transition 
se termine a la seconde 8.3 lorsque l’angle de basculement atteint 0°. 

(3) . Vol horizontal (de la seconde 8.5 a 12.5) : Apres la transition, la vitesse de 
l’avion augmente encore a une certaine vitesse, malgre que la vitesse de rotation 
des moteurs d’avant ait ete reduite a 270 rad/s (voir figure 4.31). La vitesse de 
drone ensuite converge vers une valeur constante du fait de la trainee de frotte- 
ment. Cela signifie qu’a ce moment l’acceleration longitudinale de l’avion est zero. 
Dans la figure 4.32, on peut remarquer que l’angle d’attaque converge aussi vers 
une valeur constante. Malgre que cette valeur est petite (environ 1.5°), elle suffit 
pour creer la force de portance grace a la vitesse elevee de l’avion (environ 15m/s). 

(4) . Transition vers le vol stationnaire (de la seconde 12.5 a 14-5) : Pendant la 
transition, le moteur arriere augmente la vitesse de zero a la valeur necessaire afin 
de maintenir l’angle de tangage a —2° pendant que la vitesse des deux rotors avant 
reste encore a 270 rad/s. La vitesse de l’avion diminue parce que la force Ff, se re¬ 
duit tres vite (figure 4.30). Par consequent, la force aerodynamique de portance F\ 
diminue egalement tandis que la force totale de poussee des moteurs F augmente. 
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Figure 4.28 - Chemin du vol longitudinal 
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Figure 4.29 — Les variables d’etat de l’avion pendant la periode de vol 
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Figure 4.30 — Forces generees dans le repere 


(5) et (6). Vol vertical et atterrissage (de la seconde 14-5 a 43) : Sans la force 
Ff, l’avion vole encore en croisiere parce qu’il a encore de l’elan dans le vol ho¬ 
rizontal. Pourtant, sa vitesse se reduit a cause de la force de trainee et l’angle 
de tangage —2° (qui cree la force inverse a la direction de mouvement). L’avion 
a encore de la portance aerodynamique F\ (Figure 4.30) grace a sa vitesse mais 
cette portance se reduit aussi a zero d’apres la vitesse du drone. Quand la vitesse 
atteint zero, l’angle desire de tangage est zero pour stabiliser l’avion et ensuite, le 
drone atterrit a terre a la seconde 34. 

Notons que pendant toute la periode de vol, la somme de la force aerodynamique 
de portance F\ et la force totale de poussee F est toujours approximative au 
poids du vehicule (mg ~ 12.26A) (Figure 4.30). Cela garantit que l’avion ne peut 
pas perdre l’altitude. La loi de commande d’altitude est inactivee pendant le vol 
horizontal et la transition, le prohl de la vitesse d’altitude est trapezoidal pen¬ 
dant les periodes de decollage et d’atterrissage et il n’y a pas de depassement. Le 
programme de commande de la simulation a ete developpe par le code C embar- 
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Loi de commande 



Figure 4.31 - Vitesse de chaque rotor 



Figure 4.32 - Angle d’attaque 


que (dans le bloc S-function Builder de Simulink). Ainsi, ces codes sont importes 
directement dans le programme de controle du microcontroleur. 
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4.8 Conclusion 

Le travail presente dans ce chapitre porte snr la stabilisation globale de la position 
et d’attitnde d’un corps rigide. Une nouvelle loi de commande bornee a ete proposee 
en utilisant une approche non lineaire a base de fonctions saturees. L’objectif de 
la loi de commande est de maximiser l’efficacite des actionneurs sans prendre de 
risque dans la stability du systeme. L’approche est simple et elle ne necessite pas 
la connaissance des parametres inertiels pour garantir la stabilisation d’attitude. 
La loi de commande proposee pour stabiliser un corps rigide dans l’espace a trois 
dimensions est ensuite etendue a la stabilisation de la position et l’attitude d’un 
quadrirotor et a la stabilisation d’attitude et d’altitude d’un avion convertible et 
d’un avion basculant tri-rotor. L’implantation en temps reel a ete effectuee et des 
experiences en temps reel seront presentees dans le chapitre prochain. 



Chapitre 5 


Result at s experimentaux 


Ce chapitre presente les resultats experimentaux realises au laboratoire. Les lois 
de commande presentees dans le chapitre 4 sont programmees et executees par 
le microcontroleur du systeme embarque presente dans le chapitre 2. Les donnees 
des etats mesurables sont enregistrees par la station sol dans un fichier, qui sont 
ensuite affichees par le logiciel Matlab. Les experiences sont realisees en plusieurs 
essais. Etant donne qu’apres la transition vers le vol horizontal, l’avion conver¬ 
tible volera avec une vitesse elevee, les tests sont divises en deux categories : la 
premiere concerne le vol vertical et la deuxieme la transition et le vol horizontal. 
Les videos, qui peuvent etre telechargees, illustrent les vols reels. Par ailleurs, les 
experiences du vol vertical de l’avion convertible ont egalement ete presentees lors 
de nombreuses demonstrations au laboratoire. Ces vols en pratique montrent la 
performance et Lefficacite des travaux realises pendant la these. 


5.1 Avion convertible 

Dans cette section, les experiences de l’avion convertible seront presentees. Les 
lois de commande d’attitude utilisent les quaternions (4.129 - 4.131) mais pour 
une mcilleure lisibilite des resultats, nous allons montrer l’attitude avec une repre¬ 
sentation des angles d’Euler. Cependant, 1’algorithme de commande est base sur 
un formalisme des quaternions. Les gains de la loi de commande ont ete choisis ahn 
d’obtenir une bonne stability du systeme pour le vol stationnaire et une reaction 
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rapide lorsque l’operateur souhaite deplacer le vehicule. 


5.1.1 Vol vertical 

Les experiences du vol vertical suivantes ont ete effectuees : 

- Stabilisation 

- Suivi de consigne de l’attitude pour le vol semi-autonome 

- Robustesse par rapport aux perturbations externes 

- Decollage et atterrissage autonome 


Stabilisation 

Cette premiere experience verifie l’efficacite de la loi de connnande pour stabili¬ 
ser l’avion convertible a l’origine (0 = 0°, 6 = 90°, 0 = 0^, z — z ( i)- Dans la figure 
5.1 nous constatons que l’attitude et I’altitude sont stabilises autour de l’origine. 
L’amplitude de l’erreur de l’angle de roulis et de tangage est de moins de 4°, ce 
qui garantit la stabilisation. Notons que les ailerons de l’avion sont des actionneurs 
utilises pour stabiliser l’angle de tangage et Tangle de lacet en meme temps. En 
effet dans la stabilisation d’attitude, la stabilisation de Tangle de tangage est plus 
importante que Tangle de lacet. Ainsi, c’est Tangle de tangage qui a la priorite et 
il est done stabilise d’abord. C’est pour cela que Tadaptation de Tangle de lacet 
est lente par rapport a celle de Tangle de tangage. 

La video du vol autonome ("Hands-off") a Tinterieur est telechargee sur le site : 

http://www.youtube.com/watch?v=lPL8GaaZHGk 


Suivi de Tangle de roulis et de tangage 

Dans ces essais (figure 5.2), nous analysons le comportement du systeme boucle 
lorsque l’attitude desiree est variable dans le temps. Les entrees du systeme sont 
des angles desires de roulis et tangage donnes par l’operateur via l’emetteur R/C. 
L’attitude de reference est donnee en termes d’angles d’Euler. Ainsi, les angles 
d’Euler desires sont transformes en un quaternion desire q r i qui est utilise ensuite 
par la loi de commande. En suivant l’attitude desiree, le drone peut se deplacer 
d’un point a un autre point dans Tespace a trois dimensions. 
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20 30 

time [s] 


(a) Angle de roulis 


(b) Angle de tangage 




(c) Angle de lacet 


(d) Altitude 


Figure 5.1 — Stabilisation (vol stationnaire) de l’avion convertible en vertical 


Suivi de l’angle de lacet 

La figure 5.3 nous montre l’efficacite de la loi de connnande pour maintenir en vol 
stationnaire l’avion convertible et pour suivre la trajectoire angulaire autour de 
l’axe Zb- La loi de connnande maintient le drone dans un voisinage de l’attitude 
desiree tout au long de la trajectoire (c’est-a-dire un tour 360° de l’angle de lacet). 
Les consignes de Tangle de roulis et de tangage sont donnees pour retourner le 
drone a la position initiate. Ces deux dernieres experiences montrent que l’avion 
est capable de faire un suivi de trajectoire angulaire, ce qui est une etape impor- 
tante dans la connnande d’un corps rigide dans Tespace tridimensionnel. 
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(a) Angle de roulis (b) Angle de tangage 




(c) Angle de lacet (d) Altitude 

Figure 5.2 — Vol semi-autonome : l’operateur deplace le drone en donnant des consignes 
d’attitude 


La video du vol de suivi de Tangle de lacet est telechargee sur le site : 

http://www.youtube.com/watch?v=70Agvkj_aNk 


Robustesse d’attitude par rapport aux perturbations externes 

Dans cette experience, la performance de la loi de commande pour rejeter les per¬ 
turbations externes a ete testee. Des couples perturbateurs sont introduits autour 
des deux directions de chaque axe de rotation de tangage et de roulis. Ces per- 
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(a) Angle de roulis (b) Angle de tangage 




(c) Angle de lacet (d) Altitude 


Figure 5.3 — Suivi de trajectoire de Tangle de lacet : evolution des angles de roulis, 
tangage et lacet lorsque le drone suive un trajectoire autour de l’axe Zb 


turbations provoquent des grandes erreurs dans les angles de roulis et de tangage 
comme le montrent les figures 5.4(a) et 5.4(b). Nous pouvons constater que les 
angles de roulis et de tangage sont re-stabilises apres environ 1 seconde, ce qui est 
un result at tres satisfaisant. La robustesse d’attitude est aussi garantie lorsque des 
perturbations consecutives sont introduites au systeme. 

La loi de commande a rejete ainsi de fagon satisfaisante les perturbations autour 
des deux directions des axes Xb et z/b, alors que les contraintes dans la commande 
restent satisfaites. Ces resultats montrent que la loi de commande bornee propo- 
see est robuste par rapport aux perturbations externes. Notons que la propriety 













































5.1 Avion convertible 


152 


de robustesse d’attitude est necessaire dans des missions reelles, elle garantit la 
stabilisation du drone ou le suivi d’attitude pour le deplacement. 

La video du vol ou le drone subi des perturbations est telechargee sur le site : 

http://www.youtube.com/watch?v=r5fTaXGqvEY 




(a) Angle de roulis (b) Angle de tangage 

Figure 5.4 — Robustesse par rapport aux perturbations externes 


Decollage et atterrissage autonome 

Dans cet essai, Toperateur donne les angles desires de roulis et tangage tandis 
que Tangle desire de lacet est fixe. L’altitude desiree est donnee par le generateur 
de trajectoire et la poussee totale est generee par la loi de commande d’altitude 
(4.128). La figure 5.5 represente les resultats de l’algorithme de controle d’attitude 
et d’altitude dans la phase de decollage et atterrissage. Le graphe de Taltitude 
5.5(d) montre son evolution par rapport a celle desiree. Le decollage et Tatterris- 
sage autonome ont ete testes avec succes. Notons que 19cm est Taltitude initiate 
du capteur ultrason par rapport au sol car le capteur ultrason est place a cette 
distance de l’extremite de l’avion. 

Les videos de decollage et atterrissage autonome de l’avion convertible sont te¬ 
lechargees sur le site : 

http://www.youtube.com/watch?v=iaONBY8PeOE 

et 

http://www.youtube.com/watch?v=RLoC4HvkGMs 


Le vol stationnaire a l’exterieur de l’avion convertible est montre dans la figure 5.6, 
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time (s) 


(a) Angle de roulis 



time (s) 

(b) Angle de tangage 



(c) Angle de lacet 



(d) Altitude 


Figure 5.5 - Vertical Take Off and Landing Flight Test 
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ou le drone est stabilise a une hauteur de plus de 4m . En conclusion, le vol vertical 
de l’avion convertible a ete teste avec succes. Les experiences de la transition et 
du vol horizontal seront presentes dans la prochaine section. 

La video du vol stationnaire a l’exterieur est telechargee sur le site : 

http://www.youtube.com/watch?v=UwQs56II16A 



Figure 5.6 — Vol stationnaire de l’avion convertible a l’exterieur 


5.1.2 Transition et vol horizontal 

Avant d’effectuer la transition reelle, ce qui est delicat et difficile a realiser avec suc¬ 
ces la premiere fois, des essais au sol ont ete effectues a&i de debugger le programme 
du microcontroleur. Malgre que ce code source a ete teste dans le programme de 
simulation, les essais sont obligatoires pour modiher des parametres necessaires, 
pour ajouter ou modiher le code afoi d’adapter au mieux la realite. Le mecanisme 
des tests de transition est montre dans la hgure 5.7. Ce mecanisme bloque cinq 
degres de liberte de mouvement et ne permet que la rotation de tangage autour 
de l’axe y^. La transition est realisee apres un signal de l’operateur via un levier 
dans la radio commande. Ce levier change le programme en mode de transition. 
Puis le programme est change automatiquement en mode de vol horizontal apres 
la terminaison de la transition. En mode de vol horizontal, l’avion convertible est 
pilote manuellement. Ensuite, l’operateur change le programme en mode de tran¬ 
sition vers le vol vertical via le levier de commande. Dans les transitions vers le 








155 


Resultats experimentaux 


vol horizontal ou le vol vertical, la loi de commande d’attitude genere le couple de 
tangage pour que Tangle de tangage actuel puisse suivre Tangle de tangage desire 
cree par le generateur de trajectoire (section 4.4.2). 



Figure 5.7 — Le mecanisme des tests de transition 


Les angles de tangage desires et leurs adaptations dans les essais de transition 
au sol sont montres dans la figure 5.8. Les vitesses angulaires de tangage co y ont 
une forme triangulaire pendant deux transitions. La vitesse s’annule lorsque la 
transition est terminee, ce qui evite le depassement. Notons que les periodes de 
transition sont differentes dans les deux transitions : elle est d’environ 3 secondes 
dans la transition vers le vol horizontal et d’environ 1,2 secondes dans la transition 
vers le vol vertical. Cela permet au vehicule d’accelerer pendant la premiere tran¬ 
sition et d’arriver tout de suite au mode stable (stabilise par les lois de commande 
en vol vertical) dans la deuxieme transition. Les vitesses maximum sont done aussi 
differentes : pendant la transition vers le vol stationnaire, la vitesse atteint jusqu’a 
180°/s, une valeur double par rapport a celle de la transition vers le vol en croisiere. 

Les essais a l’exterieur du vol horizontal et de la transition ont aussi ete reali¬ 
ses. La video du vol horizontal est telechargee sur le site : 
http://www.youtube.com/watch?v=vlXBheVax8A 
La video des tests de la transition est telechargee sur le site : 

http://www.youtube.com/watch?v=igw7XvhAms8 

Les essais reels de la transition et du mode horizontal ont ete effectues manuelle- 
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ment pour montrer les capacites du prototype a voler dans ces modes. Un pilote 
a perrnis de valider ces deux modes a l’exterieur. La phase de transition du mode 
horizontal au mode vertical a pourtant ete effectuee en autonome (sous la super¬ 
vision du pilote), mais les resultats n’ont pas ete transmis au sol. Seule la video 
montre cette phase. 




(a) Transition vers le vol horizontal (b) Transition vers le vol vertical 




(c) Vitesse angulaire u y (transition vers le vol (d) Vitesse angulaire u> y (transition vers le vol 
horizontal) vertical) 

Figure 5.8 — Transition entre le vol vertical et le vol horizontal 


5.2 Conclusion 

Ce chapitre presente des resultats experimentaux de l’avion convertible. Nous 
avons pu montrer que les lois de commande sont efficaces et robustes grace a 
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plusieurs tests. Notons que dans les tests a l’exterieur, le drone subit des pertur¬ 
bations dn vent. Pourtant, les videos montrent que des lois de connnande sont 
capables de stabiliser le drone. Les tests de la transition et du vol horizontal ont 
aussi ete realises avec succes. 
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Chapitre 6 


Conclusions et perspectives 


6.1 Conclusions 

D’un point de vue general, l’objectif de cette these concerne la commande d’un 
mini-drone avion a decollage et atterrissage vertical dans ses phases de transitions 
dedicates du vol stationnaire au vol horizontal et inversement. Ce travail de these 
a permis d’aborder differents domaines allant de l’automatique, de l’electronique, 
de la mecanique, de 1’aero dynamique et de l’informatique embarquee, rendant le 
travail multidisciplinaire, avec comnie application des resultats a la commande 
d’un drone convertible. Ce domaine applicatif est particulierement attrayant et les 
resultats obtenus ont ete appliques sur les prototypes developpes. 

Deux prototypes experimentaux ont ete developpes au laboratoire : l’avion conver¬ 
tible et l’avion basculant tri-rotor. Nous avons aussi fait la description de leurs 
modes de vol, le modele dynamique, congu et implements une loi de commande 
pour 1’asservissement en vitesse de chaque actionneur const it ue de l’ensemble "dri¬ 
ver - moteur - helice". L’identihcation experimentale de certains parametres et une 
recherche bibliographique sur la dynamique des drones ont permis de developper 
un simulateur numerique des prototypes congus. L’implantation en temps reel a 
egalement ete effectuee. Un grand nombre de tests ont ete realises, montrant la 
performance du schema propose par rapport au temps de reponse, rejet de pertur¬ 
bations et suivi de trajectoires angulaires. Les tests sont effectues a l’interieur ainsi 
qu’a l’exterieur. Nous avons pu voir que l’utilisation de techniques de commandes 
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relativement simples et implementables dans un microcontroleur embarque permet 
d’effectuer le vol stationnaire, le decollage, l’atterrissage, ainsi que les transitions. 

De plus, nous avons presente les lois de commande pour la stabilisation en at¬ 
titude et en position d’un corps rigide. Nous nous sommes plus specihquement 
interesses aux approches de commande bornee existantes et nous avons propose 
une commande bornee qui garantit la convergence du corps rigide vers une atti¬ 
tude et position desirees, meme dans le cas d’une grande erreur initiale d’attitude 
et de position. La loi de commande proposee est basee sur des fonctions de satu¬ 
ration et elle est tres proche par sa structure de celle utilisee pour la stabilisation 
de systemes lineaires avec une commande bornee. La preuve de stability asymp- 
totique traditionnelle utilisant la fonction de Lyapunov a montre la convergence 
de Lerreur a zero. Cette approche de commande bornee est interessante parce que 
chaque saturation ne joue que sur un seul degre de liberte, ce qui facilite le reglage 
du systeme reel. L’objectif de cette approche est de maximiser la performance des 
actionneurs sans mettre en peril la stability du systeme. De plus, le schema de 
controle permet egalement de prendre en compte des contraintes sur la vitesse 
angulaire. La principale caracteristique de notre approche est qu’elle est extreme- 
ment simple a implanter. 

Nous avons aussi constate que la commande en attitude et en position d’un corps 
rigide peut etre aussi appliquee en general a la commande des vehicules VTOL 
(Vertical Take-Off and Landing). II est remarquable que dans le cas general d’une 
plateforme VTOL, le modele dynamique peut s’assimiler a une structure en cas¬ 
cade : dynamique des actionneurs - dynamique de rotation - dynamique de position. 
La dynamique et les lois de commande utilisees de chaque etape dans la structure 
en cascade sont analysees et presentees dans les chapitres 3 et 4. Nous avons choisi 
le quadrirotor, un exemple tres populaire des vehicules VTOL, pour demontrer la 
performance de la loi de commande. 

Du cote experimental, l’electronique embarquee a ete recongue car les solutions 
employees auparavant par l’equipe etaient limitees. Le microcontroleur utilise est 
maintenant plus puissant et un systeme de commande de la vitesse de rotation 
des moteurs a ete developpe. Une station sol a aussi ete developpee, ce qui permet 
d’envoyer les parametres de vol vers le sol pour observer revolution des etats pen¬ 
dant les tests et pour debugger le programme embarque. Les problemes de pertes 
de donnees et de conflits de transmission de la communication sans hi ont aussi 
ete resolus en utilisant parallelement la radiocommande classique et la communi¬ 
cation sans hi et en realisant un protocole hable entre la station sol et le drone. La 
stability en vol stationnaire de la plateforme a ete montree lors de la Journee des 
demonstrateurs en Automatique 2010 a Angers, pour la visite du Comite AERES 
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au laboratoire, dans les journaux scientifiques et lors de visites. 

Afin d’eviter les degats dans la periode du developpement, surtout dans la phase 
de transition autonome, une programme de simulation a ete developpe, ce qui per- 
met de tester au sol sans danger tous les programmes et algorithmes sur un drone 
virtue! Ce simulateur permet aussi d’eliminer les erreurs de programmation avant 
d’effectuer le vol reel. Les programmes sont alors plus garantis et peuvent etre plus 
complexes puisqu’ils sont testes facilement sans avoir recours au drone. Les temps 
de developpement sont done reduits car tout se fait sur un seul ordinateur, sans 
avoir a reprogrammer un microcontroleur, ce qui est une operation relativement 
longue, ni avoir a attendre la recharge des batteries. 


6.2 Perspectives 

Les perspectives de ce travail sont multiples et plusieurs points dans cette these 
pourraient etre ameliores. En effet, nous allons en enoncer quelques uns ci-apres : 

Tout d’abord, les plateformes peuvent etre ameliorees. En fait, les prototypes sont 
congus et developpes par nous-meme au laboratoire. Pourtant, le domaine aeronau- 
tique n’est pas notre atout. Malgre que l’avion convertible fonctionne tres bien en 
vol vertical, le vol horizontal n’est pas toujours garanti. En effet, le fait d’effectuer 
un prototype convertible performant reste encore un defi pour nous, la solution 
proposee de la structure de 1’avion convertible ou l’avion basculant tri-rotor etant 
encore trop limitee. Une tres bonne comprehension du domaine aerodynamique 
et de nombreuses optimisations sont necessaires afin d’obtenir un prototype de 
l’avion hable. Le travail de l’optimisation du prototype est done necessaire parce 
qu’il n’existe actuellement aucune plateforme de l’avion convertible "pret a uti¬ 
lise!'" avec un systeme embarque pour le faire fonctionner. Quelques suggestions 
pour optimiser la structure mecanique actuelle de l’avion convertible sont les sui- 
vantes : le poids des avions peut etre diminue jusqu’a environ 100 grammes en 
changeant les moteurs moins puissants (et moins lourds) et en changeant les mate- 
riaux; le prohl de l’aile peut etre change par un prohl symetrique pour augmenter 
le coefficient de portance et la taille de l’aile peut etre reduite pour diminuer des 
perturbations du vent lateral en vol vertical. 

La centrale inertielle utilisee est tres precise et performante, pourtant elle est tres 
couteuse. Lin travail interessant mais ambitieux est d’effectuer nous-meme la fusion 
des capteurs inertiels (accelerometre, gyrometre et magnetometre) afin d’arriver 
a concevoir une centrale inertielle moins chere. En fait, le probleme de l’estima- 
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tion de l’attitude a ete aborde par la communaute scientifique depuis longtemps 
et beaucoup des applications embarquees ont ete realisees. Pourtant, la centrale 
inertielle faite maison doit satisfaire les contraintes de precision et les resultats 
d’attitude doivent etre fournis a une haute frequence car la commande d’attitude 
decide de la stabilisation et du comportement du drone. 

Dans cette these, nous avons propose une loi de commande pour controler globale- 
ment un corps rigide dans l’espace tridimensionnel. Les simulations sont presentees 
pour la stabilisation d’un quadrirotor. Au niveau experimental, le fait de savoir 
precisement la position de drone est encore un probleme. Le GPS nous donne la 
position avec une erreur jusqu’a 2m. Quelques methodes utilisant le systeme de 
vision pour estimer la position relative du drone par rapport a des marquages sont 
realisees au laboratoire. Les resultats obtenus sont prometteurs. Pourtant ces me¬ 
thodes ne sont pas vraiment performantes a cause du calcul lourd de l’algorithme 
du traitement d’images, la frequence de l’estimation de la position etant done tres 
basse. Recemment nous avons effectue le traitement d’images embarque sur un 
DSP, ce qui permet d’augmenter la frequence de fournir la position relative. Par 
ailleurs, l’idee est d’avoir un environnement experimental assez performant par 
exemple le systeme de capture de mouvement par cameras ( VICON ) qui permet 
de fournir la position absolue et done de realiser des algorithmes de commande 
de position. Ce systeme est cependant limite a un environnement maitrise, ou les 
conditions sont ideales pour tester des lois de commande. 
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